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1.     INTRODUCTION 

This  report  represents  the  second  phase  of  a  long-range  project 
to  study  self-adaptive  control  principles  and  apply  them  to  the 
control  of  aircraft.     The  primary  purpose  of  this  phase  is  to  develop 
a  suitable  mathematical  model  for  the  study  of  aircraft  dynamics 
including  all  degrees  of  freedom,   and  the  secondary  purpose  is  to 
test  the  feasibility  of  using  active  compensators  in  the  autopilot 
as  a  basis  for  self-adaptive  schemes  of  control. 

It  was  decided  to  restrict  the  model  to  that  of  a  rigid  airframe, 
and  it  was  also  decided  to  use  a  block  diagram  model  because  of 
well  established  correlations  between  the  block  diagram  and  a 
determinantal  array  on  one  hand,   and  between  the  block  diagram 
and  analog  simulation  on  the  other  hand.     The  initial  block  diagram 
model  is  to  a  linear  model  for  two  reasons:  firsts   it  is  necessary 
that  the  model  be  reducible  to  simpler,   well-known  linear  block 
diagram  models  when  the  proper  restrictions  are  imposed;   second, 
analysis  of  a  linear  model  is  readily  obtained  and  provides  a  basis 
for  estimating  the  importance  of  nonlinear  effects  as  they  are 
introduced.     It  is  anticipated  that  nonlinear  effects  can  be  introduced 
by  the  injection  of  properly  derived  signals,   and  it  is  expected  that 
this  modification  of  the  model  will  lead  to  simulation  studies  since 
the  complications  caused  by  the  nonlinearities  may  make  an  analytic 
treatment  impossible. 


Another  reason  for  choosing  a  block  diagram  model  for  the 
airframe  is  the  ease  with  which  an  autopilot  configuration  can  be 
added  to  the  model.     It  is  expected  that  the  self-adaptive  compen- 
sation devices  will  be  incorporated  in  the  autopilot,   though  it  is 
believed  that  the  block  diagram  model  may  lead  to  profitable 
analysis  of  the  aerodynamic  characteristics  of  the  airframe , 
with  resultant  indications  of  ways  to  improve  its  design. 


2,     DEVELOPMENT  OF  THE  STANDARD  BLOCK  DIAGRAM 
AND  THE  STANDARD  DETERMINANT. 

2.  1     Purpose 

It  is  the  purpose  et  this  seeuen  to  develop  the  mathematical 

background  leading  to  the  derivation  of  the  block  diagram  model. 

It  will  be  shown  that  there  is  a  unique  relationship  between  a 

properly  arranged  block  diagram  and  a  determinantal  array  such 

that  either  may  be  written  down  by  inspection  of  others.     The 

theoretical  development  starts  with  the  block  diagram  and  leads 

to  the  mathematical  representation;  when  applied  to  the  airframe 

one  starts  with, the  simultaneous  equations  and  arrives  at  the 

block  diagram. 


2.  2    The  Standard  Block  Diagram 

The  theory  of  simple,    single  loop  control  systems  is  reason- 
ably well  developed,,  and  the  block  diagrams  of  such  systems  are 
net  §©mgJie&t§d,     On  the  ©the?  hand  the  theei-y  ef  multiple- loop, 
multiple -input,   multiple -output  and  multicoupled  control  systems 
is  not  well  developed  because  of  the  complexity  of  the  problems 
involved.     The  block  diagrams  of  such  systems  are  naturally 
complicated  and  sometimes  difficult  to  interpret.     The  obvious 
approach  to  the  analysis  of  such  a  system  is  the  formulation  and 
solution  of  a  set  of  simultaneous  differential  equations  for  the 
system.     Perhaps  the  easiest  way  to  formulate  the  equations  is 
from  a  block  diagram  of  the  system.     In  using  this  approach 
certain  advantages  accrue  from  standardizing  the  block  diagram, 
and  certain  other  advantages  accrue  from  generalizing  the  diagram 
By  standardizing  is  meant  specifying  the  arrangement  of  the  sum- 
ming points,   blocks  and  pick- off  points,   and  specifying  a  notation 
for  the  transfer  functions.     By  generalizing  is  meant  extending 
the  standardization  rules  so  that  the  maximum  number  of  blocks 
is  clearly  specified  and  organizing  the  mathematical  approach  so 
that  the  presence  or  absence  of  any  block  is  readily  evident  in  the 
equations.     It  is  clear  that  few  systems  will  contain. the  maximum 
number  of  blocks,    so  that  the  ability  to  recognize  the  presence  or 
absence  of  blocks  from  the  equation  is  in  itself  an  advantage.     It 


is  also  clear  that  block  diagrams  set  up  from  physical  analysis 
must  be  manipulated  to  get"  them  into  the  standard  form. 

Any  system  contains  a  finite  number  of  summation  points, 
a,n<2  eeme  e£  these  are  in  a  mpin  transi-niseisn  path  in  the  sense 
that  the  signal  output  from  the  summer  is  fed  into  a  block  which 
produces  a  signal  output  closer  to  a  system  output.     If  all  other 

summing  points  are  ignored  these  main  summing  points  are 

* 

lettered  in  sequence,   a,   b,    c,   d   n,   and  the  standard  block 

diagram  is  drawn  by  arranging  these  summing  points  in  a  straight 

line  from  left  to  right,  with  intermediate  blocks -as  required..     It 

is  obvious  that  the  numbering  sequence  is  not  completely  arbitrary 

In  general  summer  "a"  i.-eds  a  block  whoce  output  feeds  summer 

"b",    etc.  ,   though  there  may  be  breaks  in  the  signal  flow  as  will 

be  shown  later.     Also  the  output  of  the  block  following  the  nth 

summer,   is  normally  a  system  output  though  additional  system 

outputs  may  be  derived  at  preceding  points. 

Each  summer  may  receive  a  number  of  signals,   but  delivers 

only  one  output.     At  each  summer  there  is  the  possibility  that 

one  of  the  signals  received  is  an  input  command  from  an  external 

source  or  a  disturbance  signal  such  as  noise.     No  distinction  is 

i 
made  here.     It  is  convenient  to  designate  such  inputs  (at  their 

sources)  by  numbers  in  sequence  1,   2,    3,   etc.     It  is  also  conven- 
ient to  relate  the  input  sequence  with  the  summer  sequence  sp 


that  input  1  feeds  summer  "a",   input  2  feeds  summer  "V\   etc. 
If  any  summer  does  not  have  an  external  input  signal,   the 
corresponding  number  for  the  input  signal  is  not  used.     In  writing 
tht  equation*  the  signal*  from  the  input  gtationa  are  designated 
by  Roman  numerals- 

The  .transfer  functions  designated  in  the  blocks  shall,   in 
general,   be  designated  by  double  subscript  notation.     To  imple- 

* 

ment  this  the  signal  at  the  output  of  a  summer  is  designated  by 
the  same  lower  case  letter  which  designates  the  summer,  and 
the  signal  output  of  the  block  which  receives  this  signal  is  desig- 
nated by  the  same  letter,   but  using  the  upper  case  form.     This 
has  the  advantage  of  clearly  indicating  all  important  signal  sources, 
and  the  disadvantage  of  permitting  two  symbols  which  represent 
the  same  transfer  function,   as  will  be  shown.     The  double  sub- 
script notation  shall  consist  of  two  letters  or  a  combination  of  a 
number  and  a  letter.     The  first  subscript  designates  the  input  to 
the  block  and  the  second  subscript  designates  the.  summer  into 
which  the  output  of  the  block  is  fed. 
Thus, 

G,    =  transfer  function  of  a  block  between  input  1  and  summer  "a" 
la 

G,    =  transfer  function  of  a  block  between  input  1  and  summer  "b" 

•        lb 

G      =  G  ,    =  transfer  function  of  a  block  between  summer  "a" 
aA        ab 

and  summer  "b" 


G        =  transfer  function  of  a  feedback  block  between  signal 
Aa,  ',,'••' 

A  and  summer  "a" 


G      =  transfer  function  of  a  feedback  block  between  signal 
B  and  summer  "a" 


G      =  transfer  function  of  a  feedback  block  between  signal 
C  and  summer  "a" 


Many  other  interconnections  are  obviously  possible  and  are 
readily  handled  by  the  double  subscript  notation. 

It  is  common  terminology  to  speak  of  systems  in  terms  of  the 
numbers  of  loops  involved.     From  the  viewpoint  of  mathematical 
analysis  the  number  of  loops  is  not  a  convenient  description  of 
the  complexity  of  the  system.     A  better -criterion  is  the  number 
of  summing  points  or  "nodes,"  and  this  terminology  will.be  used 
in  this  text. 

The  development  or  manipulation  of  the  standard  block  dia- 
gram is  most  readily  shown  by  a"  series  of  illustrations.     Figure  2.1a 
shows  the  simplest  case  of  a  single  node  system.     For  many 

practical  systems  G,      =  GA      =1.0.     All  other  variations  of  the 

1  la  Aa 

single  node  system  may  be  reduced  to  this  form  if  desired 
Figure  2.  lb  shows  a  single  node  system  with  additional  blocks  and 
summers  representing  the  insertion  of  tachometer  feedback 
and  input  velocity  feedforward.     Note  that  the  simple  rules  for 
block  diagram  manipulation  permit  reduction  of  Figure  2.  lb  to 


a)     The  Standard  Diagram 
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b)     A  block  diagram  which  readily  manipulates  to  the  c.^idard  form 


Fig.   2-1     STANDARD  BLOCK  DIAGRAM  FOR  SINGLE  NODE  SYSTEMS 
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Fig.   2-2    A  TWO  NODE  SYSTEM  DERIVED  FROM  FIG.  2- lb, 


the  form  of  Figure  Z.  la.     Of  course,   it  is  not  always  desirable 

to  carry  outv  such  a  reduction.     If  the  particular  problem  involves 

determination  of  suitable  parameters  for  G       ,  a  slight  manipula- 

Bb 

tion  results  in  the  block  diagram  of  Figure  2.2,  which  is  a  two 

node  system.     In  this  .figure  G.      is  a  new  transfer  function  obtained 
•  la 

by  manipulating  Figure  2.1b,   and  G     =  1.0. 

The  block  diagram  of  Figure  2.-2  is  in  the  standard  form,   but 

* 

is  a  special  case  since  it  obviously  does  not  contain  all  the  signals 
and  blocks  possible  for  a  two  node  system.     A  generalized  block 
diagram  for  a  two  node  system  is  shown  in  Figure  2.3.     Note  that 
each  summer  receives  a  number  of  signals,  and  it  is  necessary 
to  designate  whether  the  signal  adds  {  +  )  or  subtracts  (-).     The 
configuration' of  Figure  2.  3  is  easily  expanded  to  the  block  diagram 
of  a  n-node  system  as  shown  in   Figure  2,4. 

The  preceding  figures  have  illustrated  the  multiple  node  block 
diagram  for  cases  with  only  a  single  output  but  many  inputs.     In 
practice  any  number  of  outputs  may  exist  in  combination  with  one 
or  many  inputs.     The  basic  procedure  for  putting  the  block  diagram 
is  standard  form  is  unchanged,   but  the  results  of  the  manipulations 
sometimes  take  a  peculiar  configuration.     A  number  of  block  dia- 
grams are  therefore  presented  to  illustrate  some  of  the  possible 
configurations.  ti^' 

'  Figure  2.5  shows  a  system  with  one  input  and  two  outputs. 
The  summers  are  lettered  and  arranged  in  sequence  as  usual. 
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Fig.   2.3    GENERALIZED  BLOCK  DIAGRAM  FOR  A  TWO  NODE  SYSTEM 
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There  is  a  break,   or  cut,  in  the  main  path  since  the  output  of 
G    .    is  not  fed  into  node  b.~  Figure  2.6  shows  a  system  with 
two  inputs,   two  outputs,  and   a  reduced, number  of  paths  coupling 


2.  3  Mathematical  Analysis  Based  on  the  Standard  Block  Diagram 

The  obvious  approach  to  system  analysis  is  to  derive  a  set 
of  simultaneous  equations  for  the  system.     This  is  readily  accom- 
plished by  writing  an  equation  for  the  signal  summation  at  each 
summer  or  node.     If  this  is  done  for  the  n-node  standard  block 
diagram  of  Figure  2.4  (with  proper  attention  to  subscripts  in 
arranging  the  terms)  a  determinant  can  be  obtained  which  has  a 
sort  of  symmetry.     In  addition,   as  will  be  shown,   each  term  in 
the  determinant  has  a  specific  relationship  to  a  definite  path  in 
the  system  so  that  a  reasonable  amount  of  physical  correlation 
exists. 

Referring  to  the  n-node  block  diagram,   note  that  at  each 
node  there  is  only  one  output  signal,  but  many  signals  being  fed 
in.     Of  the  signals  being  fed  in  some  are  from  input  stations, 
some  are  feedforward  signals,   and  some  are  feedback  signals. 
From  the  nature  of  a  summing  device  it  follows  that  at  each  node 
the  output  signal  must  be  the  algebraic  sum  of  all  input  signals. 


12 


Fig.   2-6    TWO  INPUTS;  TWO  OUTPUTS 
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In  feedback vcontrol  systems  inverse  or  negative  feedback  is 
normally  used,   so  it  is  convenient  to  treat  feedback  signals  as 

subtraetive  •    Input  and  feedforward  lignalf  may  either   ad  or 

subtract,   so  it  is  convenient  to  consider  them  additive  in  the 

i  i 

generalized  treatment. 

At  node  "a"  the  output  of  the  node  is: 


a  =  -  (G       GA   )  a  -  (G      G     V  b  -  (G    ,  G.  )  c 

ab     Aa  be      Ba  cd      Ca 

-  (G       GTN)d-...-(G       G      )  n  +  G,     I  +  G„     II 
de      Da  no      Na  la  2a 

+  G.     III+  G.     IV  +   .  .  .   G 


(2-1). 


.      Ill  +   UA      IV  +    .  .  .    U        \A 
3a  4a  na  y  | 


This  rearranges  to 

(  1  +  G       G      )  a  +(G_      G      )  b  +  (G       G      )  c  + 
ab      Aa  be      Ba  cd      Ca 

m^_n  (U\) 

+  (GJ     G_  )  d  +  .  .  .  (G       GL.  J'na  ^T- 

v    de      Da'  no      Na'  <- — ,,x        (G  )m       (2-2) 

m  =  1  (I)        ma 


At  node  "b"  the  arrangment  of  signals  is  similar  to  that  at 
node  "a",  but  the  first  feedforward  signal  appears.     It  is  the 
signal  fed  forward  from  node  "a"  through  G  ,  .     The  equation  at 
node  "b"  may  be  written  as: 

b=(G     )a-(G       G      )  b  -  ( G        -(G       G      )'d 
v    ab'  v    be      Bb;  v    Cb      v    de      Db' 

-   .  .  .  (G       G__)  n  +  G1U  I  +  G.,    II  +  G_,    III  +   .  .  .  +  G  .    J/l   *** 

no      Nb  lb  2b  3b  nb  V\ 


■'■>* 
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Keeping  the  signals  in  normal  sequence  this  rearranges  to: 


-(G  ,  )a+  (1+G,    G_,  )  b  +  (G       G      )  c  +  G       G  ■    )     d  + 
ab  be  Bb  cd      Cb  de      Db 


(2-4) 


•••+<0noGNb,,1-€Z_  (Gmb,m 

.     m  =  1  (I) 


•  In' like  manner  the  equilibrium  equation  at  node  "c"  is: 

c  =  (G     )  a  +  (G.    )  b  -  (G       G_  )  c*-  (G,     G_    )  d 
ac  be  cd      Cc  dc      Dc 

(2-5) 

-   ...   -  (G       GXT  )  n  +  G,     I  +  Gn     II  G,     III  +  -.  .  .  +  G      M 

no      Nc  lc  2c  3c  nc  ¥ 


This  rearranges  to  • 

^y-a"nc)b+(1  +  Gcd  ^   C  +  (°de  GDc>'i+    '••' 

m  =  n  (IA)  (2-6) 

+(G       GM  )  n=  >>  (G       )  m 

no      Nc  r— ,_  .  me 

m  =  1  (I  ) 


By  inspecting  equations  2-2,    4,,  6,   the  symmetry  of  arrange- 
ment of  terms  is  easily  determined,   and  the  equation  at  any  node 
may  be  written  without  the  need  to  refer  to  the  block  diagram.     It 
should  also  be  noted  from  these  equations  that  all  of  the  input  signals 
or  forcing  functions  are  located  on  the  right  hand  side  of  each 
equation  leaving  on  the  left  hand  side  the  terms  which  refer  to 
components  and  signals  within  the  closed  loops.     Thus  the  left  hand 
sides  may  be  arranged  to  form  the  characteristic  determinant  of 


15 


the  system.     This  determinant  may  be  called  ZA  and  is: 


A  ■ 


1  +  G   G       +  G   ~G     .G   G      G   G   G   G 

aA  Aa       bB  Ba     cC  Ca     dD  Da    .  nN  Na 


-  G 


ab 


-G 


ac 


-G 


ad 


-G 


an 


1+G      G  G      G  G      G      G      G 

bB    Bb  cCMCb  dD   ab  nN   Nb 


-G 


be 


-G 


bd 


-G 


bn 


1+G  ^.G^,         G'^G^ G  „TG„T 

cC    Cc  dD    Dc  nN    Nc 


-G 


cd 


1+G      G      G      G 

dD    Dd  nN    Nd 


-G 


en 


1+G  __G._  ■ 
nN    Nn 


(2-7) 
The  arrangement  of  the  characteristic  determinant  has  certain 


features  which  can  be  useful: 

a)  Every  term  in  the  main  diagonal  is  of  the  form  1  +  G.   'G 
and  each  of  these  terms  is  itself  the  characteristic 
equation  of  the  basic  loop  at  the  designated  node. 

b)  Every  term  above  the  main  diagonal  represents  a  feedback 
path  consisting  of  forward  transmission  from  a  node  through 
one  block,   then  transmission  from  the  output  of  that  block 
through  one  feedback  block  to  a  prior  node.     The  node  at 
which  this  path  starts  is  designated  by  the  column,   and  the 
node  to  which  the  signal  is  fed  is  designated  by  the  row. 

c)  Every  term  below  the  main  diagonal  represents  a  feedfor- 
ward path  through  one  block  only.     The  starting  point  of 
each  path  is  the  output  of  the  node  designated  by  the 
column  and  the  node  to  which  the  signal  is  fed  is  designated 
by  the  row. 
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In  practical  systems  not  all  paths  will  exist.     When  a  path 
docs  not  exist  the  transfer  function  of  the  block  which  is  omitted 
is  zero.     On  the  other  hand  if  the  path  exists  but  there  is  a  direct 
connection  rather  than  transmission  through  a  block,   then  that 
transfer  function  is  unity.     Thus  the  characteristic  determinant 
of  any  specific  system  may  be  obtained  easily  by  arranging  its 
block  diagram  in  the  standard  form,   writing  the  determinant 

4 

of  the  general  system  for  the  given  number  of  nodes,  and  inserting 
the  transfer  functions  of  the  specific  system  at  the  proper  location 
in  the  determinant  structure. 

The  characteristic  determinant  of  any  system  is  simply  that 
portion  of  the  N-node  determinant  which  is  specified  by  the  number 
of  nodes  in  the  given  system.     For  example,   for  a  one  node  system: 
^  (one  node)  =  1  4-  G      G 


aA    Aa 


(2-8) 


for  a  two  node  system: 


J\^  (two  nodes)  = 


aA    Aa 


-G 


ab 


for  a  three  node  system: 


A  (three  nodes) 


G      G 
bB    Ba 


1  +  G^G^ 
bB    Bb 


(2-9) 


(2-10) 


1+GG  GG  GG 

aA    Aa  bB    Bb  cC    Ca 


-  G 


ab 


1+GbBGBb      GcCGCb 


-  G 


ac 


-  G 


bB 


1  +G  rGr 
cC    Cc 
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The  solution  of  the  system  equations  for  system  perform- 
ance is  obtained  with  the  usual  manipulations  of  the  determinant. 
That  is,  the  column  representing  the  input  functions  replaces 'the 
c&;vmin  far  the  desired  signal  in  thp  *ftr»ptpnqtie  aetefflUflltftt. 
and  the  resulting  determinant  is  the  numerator  determinant  Z-^ 
is  the  denominator  determinant,  and  expansion  of  these  deter- 
minants provides  the  solution  for  the  desired  signals.     If  all 
inputs  exist  and  are  activated,  then  the  column  for  the  input 
functions  is: 

m=  n(H) 


(G       )  (m) 

ma 


m=  1  (I) 


(G       )  (m) 
mb 


(2-11) 


m  =  n(V\) 


(G       )  (m) 

mc 


m  =  1  (I) 
m  =  n(|/|) 


(G       )  (m) 

mn 


m  =  1  (I) 


•   " 


This  column  replaces  the  nth  column  in  the  general  determinant 

of  equation  2-7  to  form  the  numerator  determinant.     When  any,  other 
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_  (I 


signal  is  to  be  computed  Equation  2-11  replaces  the  corresponding 
column  in  the  numerator  determinant. 

In  most  practical  cases  only  one  input  is  activated  at  a  time, 
so  that  the  terms  in  Equation  2-11  are  relatively  simple.     In  analysis 
and  design  computations  the  test  functions  likely  to  be  used  as 
inputs  are  the  step  and  the  ramp,  and  for  these  the  terms  in 
Equation  2-11  are  extremely  simple.     It  should  be  noted  that  this 
theory  has  been  developed  for  linear  systems,    so  the  principle  of 
superposition  is  valid.    Thus  even  in  the  case  x>f  multiple  inputs 
acting  simultaneously  the  system  may  be  analyzed  for  response 
to  one  input  at  a  time  and  the  results  superimposed. 
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3.     AN  ILLUSTRATION  OF  THE  USE  OF  THE  DETERMINANTAL 
APPROACH. 

3. 1    Introduction 

A  high  gain  uncompensated  positioning  servo  is  shown  in 
block  diagram  form  on  Fig.   3.1.    It  is  the  purpose  of  this  section 
to  study  the  effect  of  amplidyne  voltage  feedback,   either  positive  or 
negative,  tp  the  mag: amp' input. 

Values  of  K  and  K     have  been  determined  as 

K=  57.5  andK     =   .315 


3.2  The  Characteristic  Determinant  and  Steady-State  Error 

The  characteristic  determinant  may  be  written  down  by  inspec- 
tion of  the  block  diagram  Fig.   3.  1.     This  determinant  is: 

3 


A- 


1 


-50 


315  x  10 
s(s+100)2(s2+10s+100) 


1+    ".500        G 
s  +  200         c 


-  11,500 


s  +  200 


(3-1) 


Expansion  yields: 
£  11,500 


s  +  200 


G    +- 
c       s 


18.1  x  10 


10 


(s  +200)  (s+100)^  (s2  +  10s  +  100) 


(3-2) 


The  steady  state  error  from  a  step  input  is: 


ss 


=  Lim 
step        s  *  0 


1  G    A     A 

1  -      cd       ac 


A 


J 


(3-3) 
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Evaluating  each  term  independently  before  combining: 
From  Equation  3-2: 

Lim  &  =  1  +  57, 5G    +  905 

-  8 

8    »   0 


Lim      A 

s  ■>  0 


ac 


-  50 


1  +  57. 5G 


-57.5 


=  2876 


Lim      G    ,  =  .315 
cd     


s  >  0 


Substituting  into  Equation  3-3  yields 


&\ 


,.-:- 


ss 


step 


•0- 


.315 


(2876) 


1  +  57.  5G, 


+    905     J 


s  >  0 


ft- 225 1 

L        s  +  57,  5G   s  +  905  J 


57,  5G   s  +  9L- 

c  s  f  0 


L        57.5G     s  +  905J'     A 
c  s  >  0 


!•- 


905 


9ft".W:5s   K     (z  )  (z   )  ---  (z   ) 
-  c       1         L  n 


(pj)  (p2>  —  (pn) 


s  *  0 


This  error  will  be  zero  if  there  be  no  p  =  0. 
The  steady  state  error  for  a  ramp  input  is: 


ss 


ramp 


r  E- 


Jcd      ^^.c j 


A 


s  >  0 


.'iS" 


(3-4) 
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Substituting  values: 


E 


ss 


ramp 


s      u  1  + 


905/s 


) 


57.5  G    +  905/s 

c  s  *  0 


1  +57   Gc +  905/s  -  905/s 

1  +  57.5  G    +  905/s  ~ 

c 


s  *  0 


1  +  57. 5G 


s  +  57.5  G     s  +  905 
c 


s  ^  0 


1+57.5 


K     (z.)  (z.)  ---  (z   ) 
c        1         2 n 

(PI)    (Pt)  (P3)   ---  (Pn) 


ss'  57.5  G   s  +  905 

ramp  c 


s  ->  0 


Compensators  containing  terms  of  the  form are  not  allowed. 


Therefore, 


1  +  57.5  K 


E 


ss 


<Z1>  <*2)   —  <»n> 
c  (Pl)  (p2)  —  (pn) 

905 


ramp 
In  order  that  the  steady  state  error  be  as  small  as  possible,  the 


term 


57.5Kc(z1)(z2)-(zn) 
<Pl)  <P2)  —  (Pn)  ^ 


(3-5) 


This  may  be  best  obtained  by  having  one  or  more  of  the  z's  identi- 
cally zero.     Then  the  steady  state  error  is  the  minimum  obtainable. 
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3.  3    Summary  of  Compensation  Schemes  Which  Appear 

Feasible  and  Restrictions  to  be  Observed  When  Attempting 
Compensation. 

'   '  —  ■     '  ■■■»■■■  ■■■II     f    ■  ■ 

Compensators  containing  a  term  1/s  are  prohibited  because 

$u§  esienti&Uy  **sauita  \n  figuring  t&e  se*ve  ta  a  typo  aoi»a 

system  and  producing  a  K     =0  and  an  infinite  error  to  a  ramp 
input, 

In  all  other  compensators, 

57.5Kc(Zl)(Z2)  —  (zn)  . 

(  Pj)  (P2)  —  <Pn>     ' 

must  be  much  less  than  1.  This  is  needed  so  that  the  steady  state 
error  to  a  ramp  input  will  be  as  small  as  possible.  The  optimum 
case  will  be  for  one  of  z's  =  0. 


3.4    Compensation  Schemes  Attempted  in  the  Feedback  Path 
Around  the  Magamp. 

A  compensator  having  the  transfer  function 

K    (s2+2§        U)       s+    LOZ  ) 
—  .    .        _c z        nz  nz  . 

(S)  ~    (s  +  Pl)  (s  +  p2)  (3-7) 

was  investigated.     It  was  determined  that  if    s  =  -  10  +  j   10  were 
points  through  which  it  was  desired  to  pass  a  180  degree  locus 
and  if  p     =  p     =  230  were  selected  that  using  a  pair  of  complex 
conjugate  zeros  at  s  =  -  16  +  j  10  would  force  a  180     locus  to  pass 
through  the  desired  points. 
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Since: 


57.5  K       U>    2 

£^-£-«i 

P1P2 


57.5  K      (356) 

K     <<     2.585  (3-S) 

c 

Letting  K  =   1  does  not  satisfy  specifications,   but  it  is  a 

convenient  place  to  start.     Using  Equation  3-7  and  Equation  3-2 

yields: 

K     (s+   16  +  j  10)  (s+  16-j  10) 

(s  +  230)  (s  +  230) 

(3-9) 

11.5  x  103  s  (s  +  100)  2  (s2  +  10s  +   100) 

(s  +  197)  (s+135)  (s+60+j43)  (s+60-j43)  (s  +21+J28.5)  (s+21-j28.5) 

An  Esiac  investigation  performed  on  the  root  locus  of  Equation 
3-9  to  find  the  roots  of  the  characteristic  equation  with  this  com- 
pensator resulted  in  the  root  locas  of  Fig.   3-2.     For    K     =  2  the 

3 
system  was  unstable.     By  making  K     =6.17x10     the  system 

c 

could  be  made  stable,  but  this  violates  Equation   3-8.     This  com- 
pensator was  therefore  adjudged  unsatisfactory. 
A  compensator  having  the  transfer  function 
G(s)  =  Kc  (s  +  z) 


s2+10s+100  (S  +  P>  <3"10) 


was  investigated. 
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Substituting  Equation  3-10  into  Equation  3-2 /yields: 


K     {s  +  z) 
c 


(s     +  10s  +  100)  (s+p) 


,2       2 

llf500(s)  (s  +100)     (s   -HOs+100) 


(s+197)  (s+135)  (s+60+j43)  (s+60-j43)  (s+21+j28.5)  (s+21-j23„5) 


(3-11) 


K     (s  +  z) 
c 


(s  +  p) 


11,500  s  (s  +  100) 


(s+197)  (s+135)  (s+60+j43)  (s+60-j4.3)  (s+21+j28.5)  (s+21-j28.5) 


=  -  1 


By  examining  the  pole  -  zero  array  of  Equation  3-11  it  was 
determined  that  z  =  20,   p  =   183  would  cause  a  180     locus  to  pass 
through  s  =  -   10  +  j   10  „     Substituting  these  values  for  z  and  p 


57.3  K 


100 


(20) 
"183) 


«1 


K    <<    16.1 
c 


(3-12) 


Again  K     =1  does  not  satisfy  this  criterion  adequately,  but  it 
make  a  convenient  starting  point.     An  Esiac  investigation  (Fig.   3-3) 
was  performed  on  the  root  locus  of  Equation  3-11  and  was  found  to 
be  unstable  for  K     =  1,0.     Stability  could  be  achieved  by  setting 
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K  ^285,  but  this  violates  Equation  3-8,.   Therefore  this  com- 
c 


i 


pensator  was  adjudged  unsatisfactory. 

A  compensator  having  the  transfer  function 

/  \  2 


K 


G(s)  = 


I     8  t    Z 
{     S+   P 


s   n+  10s  +  100 
was  investigated. 

Substituting  Equation  3-13  into  Equation  3-2  yields: 

K     (s  +  z) 
c 


2  2 

(s   +10s+100)  (s+i) 


11.500  (s) (s2+  10s  +  100)  (s  +100) 


(3-13) 


(s+135)(s+197)  (s+60+j43)  <  s+60-j43)(  s+2i+j2&  .5)(s+21-j28,  5) 


K  '  (s  +z) 

c 

(s+p) 


(3-14) 


"1 


(11.500)  (s  +  100)   s 


(s+135)  (s+197)  (s+60+j43)  (s+60-j43)  (s+2i+j2,8,  5)  (s+21-j28.5) 


=  -  1 


By  examining  the  pole  zero  array  of  Equation  3-14,   it  was  determined 
that  z  =43,  p  =250  would  cause  a  180     locus  to' pass  through  the  points  ^>v 


,.«#" 


s  =  10+jlO.     Applying  the  criterion  of  Equation  3-8  yields        ^^' 

.2 


■>" 


57,5  K 


100 


43 
250 


^<1 


(3-15) 


K         </<         58.7 
c 


29 


Assume  K     =  1  satisfies  this  criterion.     Performing  an  Esiac 
c 

investigation  (Fig.   3.4)  on-the  root  locus  of  Equation  3-14  showed 
that  the  system  was  unstable  for  K    =  1 .     Stability  could  be 
achieved  for  K     *T  452,  but  thin  violates  Equation  3*15.     Therefore 
this  compensator  was  adjudged  unacceptable. 


3.  5    Ross-Warren  Method  Used  for  Compensation 

* 

From  the  root  locus  plot,  it  was  desired  to  bring  the  unstable 
roots  (s  =  21.66  +  28.  6j)  over  into  the  left  hand  plane  at  a  desirable 
location  of  s  =  -  10  +  10  j.     This  location  should  give  a  damping  ratio, 

j     ,  within  the  specifications.     This  location  should  also  give  a 
good  transient  response.     In  order  to  obtain  the  original  K     of  905, 
a  4  section  compensator  was  tried.     From  the  root  locus  plot,  the 
following  were  obtained: 

<V46° 


?  =  56. 


35 


4 

"\/K/G=  .908 


cot  0A  =   .966 
4 

esc  0.  =  1.39 

4 


A   =  cot  "     (cot  0     -        \/K/G         esc  0  ) 
=  cot-1  (.966  -  1.26)  =  cot   -1  (-.294)  =  106° 
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<V  TT  +^4  +  ?-X4 


=  180  +46-56.  35  -  106  =  63.  65 


and  since- 0    ^0    ,  this  compensator  is  allowable.     Hence: 


z  =  h 


r*  A 

"nX4  1 

sin  0„  J 
4 

=  18.908 

Vk/g 

sin  1.06 

• . 

sin  63. 65 

=  17,4 


P  = 


h  sin 


X 


sin  (04-  <Z>4) 


=   18  x  .956 
.304 


=  56.2 


and  the  compensated  system  becomes 

4 


-  1 


Is  +   17.4  \ 
Is  +  56.2  / 


11,500  (s+100)      s  (s2  +  105  +  100) '    " 

(s+200)  (s+100)2  s  (s2  +   10s  +100)  +18.  JxlO10) 


This  was  tested  for  stability  on  the  Esiac  and  the  results  are 
shown  on  Fig.    3.5.     From  this,  it  can  be    seen  that  the  system  was 
unstable  (having  two  roots  in  the  right  hand  plane) . 

The  next  compensator  investif  ated  by  this  method  was  a  3  section 
filter.     From  root  locus  plot, 


03=  61.3° 


N^/K/G=   .88 
h  =  18 
?=  56.35 


cot  0     -  .  ^47 


esc  03  =  1.14 
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X  3  =  cot"1  (cot  03-       \/K/G  esc  03) 

a  cot"1  (.547-1,012, 
=  cot"1  (-.463)  =  114.8° 
_      =  7T+  0      -  ?  -   X  3 '=  180  +  61.3  -  56.35  -  114.8 
=  70.15° 
and  ©.,>  0  ,.  which  shows  that  3  sections  is  allowable. 


->■'■ 


Hence    z=  h    \/  K/G  sin  X   =  18  x  .88  x  .907     =  15.4 

■     s  -940 

sin  © 

p  =  h  sin   X.  =  18jc  ^907  =106 


** 


<fiP 


sin  (G3  -03) 
The  resulted  in  the  following: 


.154 


5.4         3      11,500  (s+100)2  s'(s2  +  10s  +   100) 


■) 


s+106  (s+200)  (s+100)2  s  (s2+10s+100)+18.1  x  1010 


Thin  was  tested  for  stability  on  the  Esiac,   and  the  results 
shown  on  Fig.    3.6  show  that  the  system  is  still  unstable  (having 
two  roots  in  the  right  hand  plane  s=26 '+   36j) . 

3.6  Compensator  with  Zero  at  the  Origin 

A  compensator  having  a  transfer  function 

G(s)  =    Kc    S  (3-16) 

•      8  +    p 
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was  investigated.     Substituting  Equation  3- 1 6  into  Equation  3-2 

yields 

2      2  (3"17) 

11,100  s  (I  t  100)     (g  +  10s  +  100)  K   s 

, * , =  -i 

(s  +197)  (s+135)  (s+60+j43)  (s+60-j43)  (s-21+j28.5)  (s-21-j28.5)  (s+p)      , 

and  it  follows  that: 

57.5  K    (0) 

^ 4X'l 


c 

This  satisfies  the  criterion  for  any  finite  value  of  K    .     The  effects 

c 

of  varying  p  were  studied  by  selecting  several  values  of  p  and  plotting 

+** 
the  root  loci  on  the  Esiac.     It  was  determined  that  for  values  of  /p/  — , 

4,   root  loci 'from  the  two  right  half  plane  poles  of  Equation  3-17 

came  into  cancel   the  zeros  at  the  origin;  giving  either  roots  in  the 

right  half  plane  or  roots  in  the  left  half  plane  with  real  parts  —  1;  ' 

neither  of  which  is  acceptable. 

The  Esiac  investigation  further  revealed  that  an  optimum  value 

of  p  was  p  =  -  3.5.     With  this  pole  location,   there  could  be  produced 

two  pairs  of  complex  conjugate  roots,   the  real  part  of  each  being 

=  ^3.5.     K     required  to  obtain  these  root  was  K     =  8.96.     The  Esiac 
c        n  c 

plot  of  this  is  Fig.    3.7. 

However,  the  Esiac  is  not  sensitive  enough  to  determine  the 

•  * 

actual  value  of  K     and  root  locations.     Therefore  a  digital  computer 
c 

program  was  utilized  to  determine  the  exact  root  location  for  various 
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values  of  K    .    It  was  found  that  optimum  results  were  obtained 
c 

with  K    a  4.35.     This  gave  roots  as  follows: 

8  a  -  1.48  +  j  4.46  i 

s  =  -  2.58  +  j  5.925 

s  =  -  101.904 

s  =  -  97.978 

s  =  -  5.02  x  104 
These  results  satisfy  all  steady- state  and  transient  specifications 
except  settling  time  and  the  response  to  unit  load  torque.     Checking   , 
the  response  to  a  unit  load  torque,   the  block  diagram  of  the  system 
is  given  in  Fig.    3.8.     The  characterisitic  determinant  may  be 
written  by  inspection.     (G       =  0). 


A- 


50        1+ 


5x10    s 


(s+200)(s+3.5) 


0 

-11,500 

1 

s  +  200 

!06 

0 

0 

(8+ 100) 

0 

0 

0 

0 

0 

0 

0 
0 
0 


0 


1 


300i 
0      ,  0 

0         0 
9.43       A 


s+,246 


10.4 
s+10 


0     - 


9 .  09 
S+.246 
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For  the  unit. load  torque  input: 


ec 


..-  [ 


^W 


f^C 


■  to 


Obtaining  each  term  separately: 


A  ■  =  (  1  +  7.5s)  41  + 

s  ^  o  s 


36,800 


(3-18) 


fee 


i 


s  >0        300s 


A  e  =  36.9  +  2640s 

ef   s>0 

Substituting  these  into  Equation  3-18 

1 


Gc  =    (36.9  +  2640s)  300s 

88 


(1  +  71.5s)  41  + 


36.9 
(300)(36,800) 


36,800 
s 


1 


300,000 


Radian 


And  this  is 


s  >0 


ec       =   .0115  minutes  of  arc 
ss 


3>7  Combining  of  Amplidyne  Path  with  Tach  Path 

Prior  to  attempting  any  compensation  schemes  an  investigation 
was  made  to  determine  if  there  were  any  restrictions  on  the  types 
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of  compensators  that  could  be  used  in  these  paths.     Assuming 
an  unspecified  compensator  in  each  path,  the  block  diagram  of 
Fig.   3.9  is  obtained. 

From  the  block  diagram  of  Fig.    3.9,  the  characteristic 
determinant  may  be  written  down  by  inspection. 


(3-19) 


A= 


315  x  10 


-  50 


11,500 

s  +  200  cl 


11,500 
s  +  200 


s(s+100)2(s2+10s+100) 


.315  x  10   Kt  Ge? 

(s+   100)2(s2+10s+100) 


Now  for  a  step  input  at  O 


ss 


step 


R 
ac         c^c 

A 


(3-20) 


s-0 


Evaluating  each  term  separately: 


A  = 

s>0 


1 

50  1+5 

0 


.315 


cl 
-57.5 


.315K  G   _ 
t     c3 


905 
A  =  1  +  57.5  G   ,      +  18.1    K    G   0  +  T 

VA  cl  t      c3 

s>0 


>'" 


,{&> 


^ 


.«'. 


41 


0 

Q 


i  9 


,c7 


So 

3 


-J 


4- 
s5^ 


O 


o 

4. 

<^ 

-i. 


M 


& 


rJ'" 


4£ 


A„ 


=  2877 


bM) 


Gc0c 


.315 


s>  0  s 

Substituting  into  Equation  3-20 


ss 


step 


1  - 


,315 
2877  (      s       ) 


905 


1  +57. 5G   ,  +  18.1  KG,      +     s 
c  1  t    c3 


s-»  0 


E 


\  step 


905 


57.5  G   ,   s  +  18.1  K  G-      s  +  905 
cl  t    c3 


8-+-  0 


By  inspection,   it  is  seen  that  if  it  is  desired  that  E  be  zero, 

ss 
step 


then  G   .   and  G       can  contain  no  terms  of  the  form  1 
cl  c3  — 

s 
N=l  and  the  feedbacks  be  of  opposite  sign  such  that 


If  however, 


57.5  Cfc.    =-  18.11  K    G 

cl  t     c3 


-3.18  G  .,  =  K    G 

cl  t      c3 


Then  the  desired  zero  error  will  be  attained. 


Now, 


E 


ss 


RAMP 


1  - 


A       G 

ac       c0c 


s  ^  0 


(3-21) 
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Substituting  the  already  obtained  values  into  this  equation 


E 


S8 


RAMP 


,  909/8 , 

1  "l  +  47.5G   ;  +  16.1  KG     +905/s 
cl  t    03 


s>0 


• 


1  +  57.5  G   ,+  18.1  KG 
cl t    c3 

1+57.5  G   ,+  18.1  K    G-    +  905/ s 
cl  t     c  3 


s  >0 


1  +  57.5  G      +18.1K    G, 
cl t      c3 

57.5  G   ,s  +  18.1  KG  Os+905 
c  1  t    c3 


s  ^0 


From  this  equation,  it  can  be  seen  that  G    ,  and  G   _  cannot  contain 

c  1  c3 

terms  of  the  form  1       where  N  —  1    because  this  would  cause  the 

~N 
s 

numerator  to  go  to    CO  .     This  then  negates  this  possibility  of  a 

I 

solution  using  positive  feedback  in  one  path  and  negative  in  the  other 
while  using  a  term     1       as  suggested  from  the  equations  for  E 


ss 


step 


It  is  further  seen  from  this  equation  that  since  the  required 
E  = 


ss 


RAMP 


1 
905 


That  the  Lim.      G       =  Lim  G.    50 

s>0        c  1  s  >  0        c3 


Considering  now  the  possibility  that  the  response  to  a  unit  load 

torque  may  place  additional  restrictions  on  the  type  of  compensator 

adopted,  use  the  block  diagram  of  Fig.   3.8  with    4.  35  s        being 

s  +  3.5 


replaced  by  a  general  compensator  G 


cl 


44 


A 


ec 


ef  fee 


89 


u.l.t. 


s  >  0 


(3-22) 


Writing  down    ^\         by  inspection  of  Fig.   3.8  with  4.35  s    A    -, 

s-°  TTTs  =  Gci 


A    ■ 

s-0 


1 


300s 


50      1  +  57. 5G  0 


0        K.       G:_ 

t         c3 

300 


0       -57.?        1 


100      1        38.3       0 


1.04      1 


0  0        -36.9       1 


905 


=  1  +  57.5    G   ,  +  18.1K    G      +. 

cl  t      c3  s 


W 


=  36.9  +  2064  G 


s  >0 


cl 


fee 


i 

30UT 


s-0 


Substituting  these  equations  into  Equation  3-22, 


\ffi 


ec 


l/300s       ^36.9 +2064  Gc3 


88 


U.l.t. 


1  +  57.5  C.  +  18.1  KjGz-  + 

cl  .     t    c3        s 


905 


s>0 


-"      36.9  +  2064  Cr 


300s  (  1+  57.5    G  ,  +  18.1  K    G-     )  +  271,500 

cl  t      c3 
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Now  it  can  be  seen  that  if  the  previously  derived  restrictions        .,,<• 


are  observed;  i.e.   Lini  ^Q     CQl     G^  =  Urna^0    Q        =  ( 

c3 


•jf^ 


i 


then     ©c  =    •TTrz — zrrrr     -      ^-.^       radian     =  .468  minutes, 

ss   ,x  271,500  7530 

ult 

Now  that  restrictions  on  the  type  of  compensation  have  been 

determined,  it  is  useful  to  put  the  equations  of   l-±  into  a  convenient 

form  for  working  with  G       and  Q-  ,    simultaneously.     It  is  impossible 

to  obtain  a  form  of  the  characteristic  equation  in  which  G   ,  and 

cl 

G       appear  simultaneously  as  cascade  compensators.     It  is  possible 
however  to  obtain  an  equation  of  a  root  locus  which  will  be  of  the 
form 

k  fr  <.♦.;>  cc3 

-  i  =    — l_,=  i 


n 


i,k2  .■;    (.+  Zj)Gcl 


m 


n 


TV    (s  +  pJ 
k=i  K 


or  conversely: 


n 

'  I         (s  +  z.)n 


,  K3i.'l      <9+*j>"GcI 

n  ,;. 

1  +       KA  y[   (s  +  z.)  G   , 
4,j.l  }'      c3 

n 

8    m    IT      (8+   P.) 

k.l  J 
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Now  clearing  the  denominator: 


18,lxi010(8+100)a(8a+10s+100)  G   .  +  36.2xl08K  sG  , 

A '  -  ■  cl  "*    c3 


=  1  + 


s  +4105+5.41xl04s4+2. 54x10    s3+2. 5xl07s2+2xl08s+18. 1x10*° 


Letting  this  6th  order  equation    =  D  a 


18.1xl010(s+100)2(s2+10s+100) 

A=  i  +  : 


cl     cl 


DA 


36.2x10   K  s  N   6 
t        c3 


DA  D   „ 

A  c3 


where 


G.=K.N  ;        G   .  =    N  ' 

cl  cl         cl  c3  c3 


D   .  D   . 

cl  c3 

i 
Now    Z^  may  be  put  into  one  of  the  two  previously  mentioned  forms 

The  first  form  is  : 

36.2x10       KsN   „ 
t      c3 


A=l  +  - *- 


s(s+100)(s2+10s+100)N   .  K' 

cl      cl 

1    +      "H- 


DA  Dcl 
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Clearing  the  denominator 

36.2xl08  Kb  N-.D  , 
y\     .  ,  . t      el   el ;         

DAD   ,  +  8(6+100)   (s  +108+100)  K   ,N  D  ,  }*-**9 

L*-»    cA  cl    cl  c3 


The  bracketed  term  in  the  denominator  may  now  be  used  to  form 
the  root  locus  equation 


s(s+100)2(s  +10s+100)  K   ,N   . 
A    =  i  +  cl    cl 


^ 


Thus  it  is  seen  that  the  effect  of  the  compensator  G   ,  is  to  supply 

cl 

zeros  in  Equation  3-23  for  each  pole  of  G       and  also  to  provide 

cl 

the  poles  of  Equation  3-23  in  the  form  of  the  roots  of  Equation  3-24. 

The  alternate  form  of  the  equation,  by  the  principle  of  similarity 
is: 


A                18.1xl010  (s+100)2(s2+10s+100)  D'      N  '  K   , 
A    =  i  + c3      cl     cl 

A*     D   *    +  36.2xl08  K    a  N  J  D   _ 
D        c3  t         c3       cl 

Of  which  the  denominator  bracketed  term  may  be  taken  to  form 


the  root  locus  equation 

36.2xl08  K  s  N   , 
A      -  *        c3 

D  c3 
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A  4.35    g 

Assuming  a  compensator      G 


cl       s  +  3.5     ' 

This  function  for  G   .  is  the  compensator  which  almost 

cl 

met  specifications  previously.     It  is  therefore  taken  as  a 

starting  point  to  find  aQ  .  which  will  actually  meet  specifications. 

c-j 

Substituting  this  equation  for  G       into  Equation  3»23: 

36.2  x  102  K    s  (s+3.5)     G   „ 

At  c3 

=  1  + 


Roots  obtained  in  Section  VII  (e) 


It  is  seen  that  the  pole  zero  array  of  Fig.    3.10  will  result.  ,,t^" 

:  "  ■?-'  '  ' 

Now  G    _  must  be  a  compensator  containing  a  term  (s)  in  the  ;^'" 
■  c3 

numerator.     This  will  result  in  a  double  zero  at  the  origin  of 

Fig.   3.10.     Obviously,  a  root  locus  will  always  appear  in  the  region, 

between  one  pair  of  the  complex  poles  at  s  =-2.5  and  the  origin. 

Since  this  will  always  give  a  set  of  roots  with  their  real  parts  ^"-2.5  , 

this  compensator  will  not  produce  satisfactory  results  in  terms  of 

transerit  response  settling  time  specifications. 

It  can  be    seen  that  to  achieve  satisfactory  compensation  with 

the  assumed  G      ,   it  will  probably  be  necessary, for  G       to  cancej 
cl  c3 

all  four  complex  poles  at  s  =  -  2.5.  If  the  compensator  G  -  must 
be  this  complicated,  it  would  be  just  as  well  to  use  a  complicated 
compensator  for  G  and  achieve  compensation  by  use  of  the  mag 
amp  path-  alone . 
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It  may  be  seen  further  from  the  sketch  of  Fig.    3.11,   the 
root  locus  of 

1     2 

A     .  i ' x  UtS008(B+i0.0)   (s  ^108+100) 

**A!H  u^J..,uu^.u^y.f..  > '-!--  '-■■-'■■'     Gcl  ,     (3.26) 

D 

where  it  is  required  that  G       contain  a  term  (s)  in  the  numerator 

cl 

that  there  will  always  be  at  least  two  roots  of  Equation  3-26  in  the 

i  4 

area  between  s=20  and  s  =  -  5;  since  the  root  locus  from  the  right 
half  plane  poles  will  always  go  to  either  the  zeros  at  the  origin  or 
at  s  =  -  5+  j  8.6.     These  roots  would  then  become  poles  of  the 
equation 

A  36.2  x  108  K    s  D    , 

A=     t  cl     G 

B  n  c3 


sl'      (s+pi) 

and  since  as  mentioned  previously  G    „  must  contain  a  term  s 

c3 

in  the  numerator  there  would  in  all  probability  be  root  locus  between 

the  poles  under  discussion  and  the  double  zero  at  the  origin;  unless 

of  course,   G   _  produced  zeros  to  cancel  those  poles.     Since  it  is 

necessary  to  retain  the  zeros  at  the  origin,  the  obvious  alternative 

is  to  cancel  the  zeros  at  s  =  -  5  +  j8.6  with  an  active  compensator 

for  G   .  •     If  this  be  done,   all  specifications  can  be  met  using  the 
cl  . 

mag  amp  path  only  and  no  additional  advantages  accrue  from  also 
compensating  with  the  tach  path. 
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The  denominator  of  each  of  these  root  locus  equations  is 
another  root  locus.     The  denominator  root  locus  is  a  function  of 
only  ons  compensator  thus  one  compensator  may  be  used  to  deter; 
mine  the  poles  which  the  other  compensator  will  work  with.     In 
many  cases,  this  may  be  a  useful  tool.     Obtaining  these  forms 
for  the  two  compensation  paths  now  under  consideration: 

Expanding  the  determinant  of  Equation  3-19t 


A      .  18.1  x  1010 11,500 

^=  1  + _ +     G 

s(s+100)    (s+200)(s  +10s+100  s+  200 


(s+100)2(s+200)(s2+10s+100) 


=  1  + 


11,500 
s+200        cl 


18. lx  101 


s(s+100)2(s  +200)(s2+10s+100) 


36.2xl08K  G  ,' 
t  c3 


(s+100)2  (s+200)(s2+10s+100) 


18.1  x  1010 
1  +  


s  (s+100)2(s+200)(s2+10s+100) 


53 


+ 


8  g0*p 

36.2x10  KG  _ 
t  c3 


:.:-■' 


When  using  an  active  compensator  in  the  mag  amp  path  to 
cancel  the  zeros  at  s  =  -  5  +  J8.6,  as  a  first  consideration  this 

compensator  must  be  cascaded  with  a  compensator  of  the  form 

s 

An  active  compensator  may  be  used  to  produce  poles 


s  +  p 

at  s  =  -  5  +  j8.6.   'There  will  be  two  real  zeros  which  come  as 
an  integral  part  of  this  compensator,  which  may  be  placed  as 
desired.     The  farther  out  on  the  real  axis  these  zeros  are  placed, 
the  lower  the  root  locus  gain  will  become  to  place  roots  in  the   . 
desired  locations.     This  presents  no  real  problem  since  a  voltage 
divider  network  may  be  made  part  of  the  compensator  and  only  a 
part  of  the  mag  amp  output  feedback  as  compensation. 

Therefore,  for  simplicity,  let  the  two  real  zeros  cancel 
the  poles  at  s  =  -  197  and  s  -  -  135.     This  gives  the  equation' 


ll,500s(s+100)2(s2+10s+100)K   ,s    (s+135)(s+197) 

A- ?1 -2 +i 

&D  (s+p)    (s  flOs+100) 


11.500     Kcl  .'*    <s+100>2 


(s-21+j28.5)(8-211j28.5)(s+60+j43)(8+60-j43)(s+p) 


+  1 
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It  was  determined  from  Esiac  investigation  that  p  =  -  3  would 
produce  a  root  locus  through  the  desired    0  range.      (See  Fig.   3.12) 
The  root  locus  gain  to  produce  roots  with  real  part    =7,5  was  102. 
The  compensator  then  becomes 


102' 

(s  +  135)  (s+197)  s 


r  11,500 

Gci=7^~" 


(s     +  10s  +  100)  (s+3) 

,  4 

Substituting  this  into  Equation  3-26  and  setting     l\  =  0,  the 

characterisitic  equation  polynomial  becomes: 

0  =  s5  +  181.68s4+22,035.85s3 +  1 ,095, 320. 9752+8, 506,777 .62s+ 

25,070,121 
Solving  this  equation  on  the  digital  computer  gave  the  following 
roots: 

s  =  -  4.47 +  2. 88j 

s  =  -  53.28+  10 3j 

s  =  -  66.18 
Investigation  of  these  roots  showed  that  the  required  specifica- 
tions were  not  met.     Several  values  of  gain  were  used  until  satisfactory 
roots  were  obtained  to  meet  specifications.     The  value  of  K  used 
was  39.0  giving  the  following  roots: 

s  =  -  42 . 6  +  60j 

s  =  -  1 1 . 2  +  20j 

s  =  -  16.0 
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,    These  roots  were  obtained  from  the  root  locus  plot  shown  in 
Fig.   3, .13.     Since  the  compensator  satisfies  the  steady  state, 
criteria  and  these  roots  satisfy  the  transient  specifications,  this 
compensator  was  selected.     From  Fig.   3.  13,  it  can  be  seen  that 
the  dominant  (s  =  16.0)  root  is  far  enough  away  from  the  origin 
to  give  a  settling  time  within  the  required  specifications.     A 
damping  coefficient  of    o   =   .485  resulted  from  the  above  roots. 
Since  the  inclusion  of  the  tach  path  offered  no  advantages  in 
simplification  of  the  compensator  required,  it  was  not  investigated 
further.     Analyzing  the  transient  response  gives: 


c        10.  1  x  1010s  (s+8)  (s2+  105  +  100) 


°R        (s+16)  (s+11.2  +  20j)  (s  +42.6  +  60j) 


From  this,   it  can  be  seen  that  the  settling  time  is  4  x 


11.2 
.357  seconds  which  is  well  within  the  required  specifications. 

This  detailed  illustration  is  intended  to  shew   some  of  the  uses 
Of  the  determinant- block  diagram  relation.     The  solution  obtained 
is  not  a  complete  one,  nor  does  the  problem  itself  have  any  bearing 
on  the  self-adaptive  problem  other  than  as  an  illustration  of  manipu- 
lative methods. 
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4.      DEFINITION  OF  SYMBOLS  USED  IN  FOLLOWING  SECTIONS 

A 
U       =      x  component  of  velocity  (scalar) 

V       =      y  component  of  velocity  (scalar) 

W      m  m  i  ecmifrohont  Hi  VtiOtiUy  (ie*l*f) 

A  2 

A      =      Moment  of  inertia  about  x  axis,   slug-ft    (whole  mass) 

A  '  2 

B       =      moment  of  inertia  about  y  axis,   siug-ft     (whole  m.^ as) 

C       =       moment  of  inertia  about  z  axis,   siug-tt    (whole  mass) 

D       =       product  of  inertia,    j      yz  dm,   a  lug- ft 
A  ,  „.  '    /•  ,  ,        ..2 


E 


=  product  of  inertia,    J     xz  dm,   slug-ft 

A  r  2 

F  ,==  product  of  inertia,    J     xy  dm,   siug-ft 

P       =  angular  component  of  rotation  about  x-axis  (scalar)  radiao./sec 

Q       =  angular  component  of  rotation  about  y-axis  (scalar)  radian/ sec 

R       =  angular  component  of  rotation  about  z-axis  (scalar)  radian/ sec 

A 

K       =  relative  angular  velocity  of  rotor  term  on.  x-axis  (scalar) 

K_     =  relative  angular  velocity  of  rotor  term  on  y-axis  (acalar) 

A 

K_     =  relative  angular  velpcity  of  rotor  term  on  z-axis  (scalar) 

A  \ 

a       =  moment  of  inertia  of1  rotor  along  x-axis 

b        =  moment  of  inertia  of  rotor  along  y-axis 

c        =  moment  of  inertia  of  rotor  along  z-axis 

A  2 

a        =  acceleration  of  airframe  at  its  mass  center  ft/ sec 
c 

m      =  mass  of  control  surface 

e        as  eccentricity  of  control  surface  between  hinge  and  the  control 
mass  center 

A 

0       =  angular  displacement  about  the  reference  (x  )  axis 
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5.     DERIVATION  OF  THE  STANDARD  BLOCK  DIAGRAM  FOR  THE 

AIRFRAME. 
5.1    Introduction 

The  types  of  motion  that  result  from  a  disturbance  in  some 
equilibrium  flight  condition  and  the  transient  motion  of  the  aircraft 
in  response  to  control  movements  are  analyzed  in  the  study  of  the 

dynamic  characteristics  of  an  airplane.     Over  the  course  of  years  the 

* 
effects  of  dynamic  stability  in  aircraft  have  been  known  to  designers. 

Tj'he  mathematical  computations  were  discussed  by  early  aeronautics 
pioneers  which  included  Lanchester ,  Bryant  and  Glauert.     In  spite  of  the 
fact  that  dynamic  stability  effects  were  known,  most  people  paid  little 
concern  to  this  field.     As  design  groups  solved  other  problems  of  air- 
craft design,   the  dynamic  stability  effects  were  automatically  solved  as 
a  secondary  result.     However,    since  World  War  II,   aircraft  development 
has  pushed  into  the  field  of  very  high  performance  aircraft.     In  many 
cases  the  dynamic  stability  problem  doesn't  solve  itself  as  happened 
previously.     Furthermore,  the  fire  coiitrol  problem  has  entered  the  scene 
which  has  generated  a  control1  problem  of  no  small  magnitude. 

This  work  will  consider  the  aircraft  in  six  degrees  of  freedom  and 
explore  paths  of  compensation  and  solution.     It  is  the  prime  intention  to 
explpre  the  work  of  Chu  (1)  in  multi  loop  servo  systems  and  determine  the 
value  of  it  as  applied  to  the  aircraft  stability. 

The  first  section  will  set  down  a  general  definition  of  terms  and 


«**•" 


th^e,  basic  equations  of  motion  as  related  to  aircraft.     It  will  further  set  ;•-' 

1 1 
down  the  basic  assumptions  upon  which  the  equations  are  formed. 
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The   equations  motion  for  the  airframe  which  are  completely 
derived  in  Appendix  I  are: 

* 

F    =  (U  +  QW  -  RV)  m  =  X  -  mg  sin  e  (5-1) 

» 
F    t  (V  +  RU  -  PW)  m  =  Y  +  mg  cos  0  sin  V  (5-2)    ,.^ 

7  ^'" 

F    =  (,W  +  PV  -  QU)  m  =  Z  +  mg  cos  0  cos  V  (5-3) 

z  ... 

L,  =  AP  +  aK}  +  (C  -  B)  QR  +  D  (Q2  -  R2)  -  E  (PO .+  R) 


+  FPR  +  cK3Q  -  bK  R  (5-4) 

4 

M  =  BQ  +  bK2  4  (A  -    C)  PQ  +  E  (P2  -  R2)  -  FQR  +  DPQ 


+  aKjR  -  cK'  R  (5-5) 

i  ' 

N  =  CR  +  cK3  +■  (B  -  A)  PQ  +  F  (Q2  -  P2)  +  E  (QR  -  P) 


-  DPQ  +  bK2P  -  aKxQ  (5-6) 


H+Fslri^mla      =P       (RP  -  Q)  (5-7) 

e         e        e  • c  e   e   cz        ex 


H    +  F     »  I  >     s  m   1   a       -  (R  +  PQ)  P       -  (RQ  -  P)  rz  (5-8) 

r  r        r»  r   r   cy  rx 

2H    +  F     =1    X     +2P       (RQ  +  P  )  (5-9) 

a         a        a  A  aY 

The  angular  and  linear  velocity  of  a  free  body  in  space  are: 

P  =   (J  -     N*    sine  (5-10) 
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Q  =  0  cos  0    +     V    sin  0  cos  ©  (5-11) 

R  =    V    cos  0  cos  9  -  0  sin  0  (5-12) 


-™  *  U  cos  0  cos    ty   +  V  (sin  0  sin  0  cos  ty  -  cos  0  sin  ty   ) 
at 

•     +  W  (cos  0  sin  0  cos  V  +  sin  0  sin  V   )  (5-13) 


-r^U   *  U  cos  9  sin  ty   +  V  ( sin  0  sin  0  sin  ^  +  cos  0  cos  V  ) 
at 

+  W  (cos  0  sin  0  sin  ty    -   sin  0  cos  ty  )  (5-14) 


dz* 

-< —    5  -  U  sin  ©  +  V  sin  0  cos  0  +  W  cos  0  cos  ©  (5-15) 

at 


Consideration  of  the  aerodynamic  forces   must  be  resolved  into  the 
effects  they  produce  on  various  parts  of  the  body  in  some  flight  condi- 
tions and  also  the  resultant  moment  effects  they  produce. 

The  first  or  necessary  condition  is  to  establish  a  reference  condi- 
tion and  proceed  from  there.     A  reference  condition  is  some  flight 

■  '  i. 

condition  where  equilibrium  exists  in  all  equations  and  the  summation 

of  forces  is  zero.     This  is  a  steady  state  condition. 

The  set  of  axis  upon  which  this  rigid  body  is  oriented  is  commonly 

called  the  Eulerian  Axis  system.     General  aerodynamics  usage  modifies 

the  z  axis  so  the  positive  direction  points  downward. 
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The  first  rotation  moves  about  the  Z  axis  in  a  direction  4)  and  is  known 
as  yaw.     The  next  is  the  movement  about  the  Y  axis  which  is  term  0 
and  known  as  pitch.     The  third  motion  about  an  axis  is  roll  about  the  X 
axis  and  ia  maneuvered  in  terms  ef  0.     Caution  must  be  used  in  deter- 
mining these  positions  because  the  order  in  using  this  system  is  not 
commutative . 

The  derivation  for  the  yaw  rate  R,  the  roll  rate  P  and  the  pitch 
rate  Q  are  shown  in  detail  in  Appendix  Two.     They  are: 

P=0-Vsin0  (5-16) 

Q  =  0  cos  0  +  lf>  sin  0  cos  0  (5-17) 

R  =  ^  cos  0  cos  0  -  0  sin  0  (5-18) 

and  where  small  angles  are  involved 

P~0  'I  (5-19) 

Q~  0  (5-20) 

R^V  i'  (5-21) 

Similarly,  Etkin  shows:  1 

0  =  Q  cos  0  -  R  sin  0  (5-22) 

0  x  P  +  Q  sin  0  tan  0  +  R  cos  0  tan  0  ($-23) 

Y  =  (Q  sin  0  +  R  cos  0)  sec  0  (5-24) 

The  general  equations  of  motion  have  to  be  applied  to  the  aircraft  motion. 

t 

The  effects  of  axis  movement  upon  the  weight  force  of  the  aircraft  will  .-^ 
be  shown  as  well  as  the  reference  state  condition.    It  must  be  remembered 
that  the  gravity  force  is  always  oriented  toward  the  center  of  the 'earth. 
The  reference  body  axis  may  or  may  not  be  oriented  in  this  direction. 

68 


Origin  at 

center  of       q 
gravity 


t 


NORMAL  AIRFRAME  AXIS  ORIENTATION 


a 


<K* 


From  the  discussion  of  the  Eulerian  angle  concept,  it  is  deter- 
mined that  a  reasonable  reference  condition  would  be 
F    =  O  =  X    -  mg  sin  ©0 


ry  =  V     -  W  eee  &     sin  0 
o  o  o 

F  s  0  =  Z     cos  e     cos  0 
z  o  o  o 


(5-25) 

(5-27)"' 
This  allows  for  a  steady  state  condition  where  the  aircraft  is  not 

accelerating  but  neither  is  it  necessarily  in  straight  and  level  flight. 

« 

From  the  previous  illustration  of  Eulerian  angles  coupled  with  the  deri- 
vations in  Appendix  II,  the  gravity  forces  can  be  determined  as  shown 
below. 


cos  0  cos  lp 
cos  y  £  ain  0 
-sin  ty    cos  0 
cos   V   sin  ©  cos  0 
-  sin  W    sin  0 


cos  ©  sinty 

cos    Ipcos  0 

+  sir>  ^P  sin  ©  sin  0 

sin  lpsin  ©  ctis  0 

-cos  V   sin  0 


-sin  0 

cos  ©  sin  0 

cos  ©  cos  0 


—          — 

my 

X 

xmg= 

my 

y 

■ 

my 
z 

'Elscaeaipn  and  derivation  of  this  general  cosine  matrix  are 
found  in  several  sources  (2,3).     The  overall  result  of  the  equations 
of  motion  is  shown  below. 
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m  (U  +  QW  -  RV)  =  X  -  (mg  sin  ©  )  (cos  e  cos  V  )  +  mg  (cos  ©     sin  0  ) 

o  .  •  o  o 

(cps  6  sin  V  )  -  mg  (cos  ©     cos  0  )  (cos  ©  sin  0  )  (5-28) 

o  o 

t  t 

m  (V+  RU  -  PW)  ■  Y  -  mg  (tin  ©fi)  (coi  V  nn  e  lin  0  •  lin  V  cos  0) 

+rng  (cos  0     sin  0    )  (cos   ^    cos  0  +  sin   ty     sin  ©  sin  0) 
o  .  o 

+mg  (cos  ©    cos  0  )  (cos  ©  sin  0)  (5-29) 

o  o 

m  (W  +  PV  -  QU)  =  Z  -  mg  (  sin  ©  )  (cos  W  sin  ©  cos  0  -  sin  V  sin  #) 

+  mg  (cos  ©     sin  0    )  (sin   V  sin  ©  cos  0  -  cos   ty  sin  0) 
o  o 

+  mg  (cos  ©     cos  0  )  (cos  ©  cos  0)  (5-30) 

o  o 

,1 

L  =  AP  +  aK    +  QR  (C-B)  -E  (PQ  +  R)  +  cK_Q  -  bK  R  (5-31) 

1  5C 


M  =  BQ+  bK,  +  PQ.(A-C)  +  E  (P     -  R  )  +  aKR  -  cK.P  (5-32) 


_L] 

z 


N  =  CR  +  cK3  +  PQ  (B-A)  +  E  (QR  =  P)  +  bK2P  -  aK}Q  (5-33) 


The  above  equations  are  complete  except  for  the  external  forces 
X,   Y,  and  Z  and  the  external  moments  L,  M,  and  N.     These  forces  and 
moments  are  composed  of  aerodynamic  and  propulsion  effects  and  the 

moments  due  to  control  surfaces.     Thefse  equations  are  very  general  in 

i 

nature  and  contain  the  following  assumptions  for  simplification  from  the 
most  general  and  complex  case.     These  assumptions  are: 

1.  The  earth  is  assumed  to  be  fixed  in  space  and  the  Earth's 
atmosphere  is  assumed  to  be  fixed  with  respect  to  the  earth. 

2.  The  airframe  is  assumed  to  be  'a  rigid  body. 
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3.  The  mass  of  the  body  is  assumed  constant  during  the  time 
of  the  problem. 

4.  For  an  aircraft  or  missile,  the  XZ  plane  is  assumed  to  have 
mirr-e?  symmetry. 

The  above  equations  then  cover  the  six  degrees  of  freedom  found  in  an 
aircraft  or  missile.    Coupled  with  the  other  auxiliary  equations  covering 
roll,  pitch  and  yaw  effects  along  with  control  motions,    solution  is  pro- 
bably  possible.     To  date,  this  has  been  accomplished  either  by  a  machine 
using  specific  coefficients  for  various  terms  or  several  simplifications. 
Discussion  of  this  will  be  conducted  in  succeeding  sections. 
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5.2         Consideration  of  External  Forces  and  Moments  on  a  Body  in 

Motion  | 

In  the  previous  section  the  equations  of  motion  were  discussed. 
The  aerodynamic  thrust  and  control  effects  were  simply  grouped  in  X, 
Y,   Z,   L,  M  and  N.     It  is  the  aim  of  this  chapter  to  discuss  these  in 
detail  with  respect  to  the  effect  on  the  dynamic  control  stability  problem. 

Consider  the  term  X  when  the  aircraft  is  in  a  disturbed  condition. 

* 

The  reference  condition  is  X    and  thus  X  =  X    +  Ax  for  some  disturbed 

o  o 

condition.     The  question  is  then  to  determine  the  ^  X  components. 

The  force  and  moment  variables  can  be  expressed  in  a  Taylor's 
Series    .  I 

F  =  Fo  +  ,_££_, *♦  <^§-  >  -£-  -  ♦  4j->/>&   ^      M 

For  practical  reasons  the  effects  of  second  and  higher  order  terms  are 

i  -i  1 1  ;■:       -li 

omitted  from  further  considerations  as  will  be  discussed  later.     Con- 

I:  :  i  . ' 

sidering  the  motion  of  the  aircraft  itself  the  X  forces  would  be  .  -    >_. 

v      v         ./^X.  ,  0  X  v  *       .  &X  .  .  #  X.   •       .  0  X .  .  ^X .   « 

x  -  Xo  +  <-fv» u  +  *tt>  u  +  'tv'  w  +  (t** w  +  ('7T)  q  +  {TT)  q 

/  ^x  .            ^x  .   .         #  X .         .  0  X  »  •      .    <^  x  .         ,  <ff  x  .    • 
+  (— )  v  +   (  )  v  +  ( )  p  +  ( j-)  p  +  (— )  r  +  (-*£-*-  )    r 

Similarly,  I  '  (g_36) 

Y .  Yo  +  (4i )  u + (z^ )  u  +  '\£t )  w + (4i)  w(+  h^X)  q + (4?, • 

^u  ^u    '  Vw^wV<l^ciq^ 

i.i  '  ,;::■■■' 

For  convenience  all  terms  will  be  expressed  in  letter  subscript  form*.* 

It  can  readily  be  seen  then  that  considering  X,   Y,   Z,   L,  M  and  N  and 
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the  effects  pf  U,   V,  W,  P,  Q,  and  R  plus  the  acceleration  terms  U, 
V,  W,   P,  Q,     and  R  there  are  72  derivative  terms  required.     This  does 
not  include  the  control  effects  which  can  add  another  36  derivatives. 
Needless  to  say,  considerable  effort  has  been  expended  in  justifying 
Simplifying  assumptions.     In  addition  to  this,   obtaining  reliable  co- 
efficients  for  various  derivative  terms  is  not  always  possible. 

It  is  not  the  intention  of  this  work  to  cover  the  computational 
methods  for  these  various  derivatives.     The  scope  is  too  large.     A 
siearch  of  all  authoritative  books      3,   4,   5,   6,   on  the  subject  yields  a 
table  of  terms  which  can  be  reasonably  determined.     This  appears  as 
Table  1.     These  considerations  are  mentioned  to  bring  out  the  awareness 
of  the  different  causes  and  effects  of  external  forces.     There  are  con- 
siderable gaps  in  the  table.     Determination  of  these  to  complete  the 
table  would  be  a  topic  of  a  complete  study  in  itself.     Much  work  remains 
to  be  done  in  this  field. 

i  Furthermore,   some  of  the  methods  of  computing  the  known 

stability  derivatives  are  questionable.     Some  computations  from  theory 
prove  quite  accurate.     Others  are  very  poor  and  empirical  methods  or 
model  testing  are  used.     Even  these  fail  at  times  and  "educated  guesses" 
from  previous  airplanes  are  used  as  a  basis  of  determining  stability 
derivatives  for  a  new  design.     Tail  effects  are  the  most  notable  in  this 
last  category. 
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In  all  aerodynamics  problems,   considerable  effort  is  made  to 
non-dimensionalize  the  equations.     The  advantages  of  doing  this  are 
of  considerable  merit.     For  the  most  part,  it  is  rather  simply  done. 
It  noii  been  iquna  In  &etn  tfeeery  and  aetual  praetie*  that 

F  =  1/2  /?V2SC_  (5-37) 

Where  S  is  a  reference  area,  usually  the  v  ing  area,   V  is  the  absolute 

reference  velocity  and  C      is  some  non.dimensional  coefficient.     C 

varies  with  Mach  Number  and  Reynolds  Number  depending  on  the  force 

considered,     (lift,   drag,    etc.)    P    is  the  density  of  the  fluid  through 

which  the  body  is  passing.     Thus  a  term  such  as  Zw  has  the  following 

relationship: 

(5-38) 

Zw'  =  #*-    =   1/2  />U   2C  (-^-?)    =   1/2    /?U'  SC         =   l/2/7SU    C7 
^  w  /        o  S  v  3  w  /        o        Z  f         o    Z 

Here  again  V  and  S  are  the  velocity  and  wing  area  terms  respectively. 
The  same  concept  applies  to  the  moment  terms  except  that  a  moment 
arm  term  must  also  appear  in(  the  equation.     A  generalized  example  is: 

M  =  /?U2    SIC  (5-39) 

I        Q  m 

\  ■  ■  .  i 

In  this  example  M  is  some  moment  and  C      is  a  needed  dimensionless 

M 

coefficient  to  satisfy  the  equation.    P  ,  U     and  S  are  the  same  as 
mentioned  above  while  Jo  \b  the  moment  arm. 

Likewise  in  a  stability  moment  the  derivative  can  be  of  the  following 

.  '  ■       s  '  :'  I' 

form  J 
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M     =       SU2l(-^)h=       U'     SiC  (5-40) 

h  o         3  h  .  o  m, 

i  n 

Another  external  force  to  be  considered  is  that  of  thrust. 
This  becomes  an  important  factor  in  most  flight  conditions.  ,  The 
location  of  the  thrust  axis  is  seldom  coincident  with  any  reference 
axis  in  the  body. 


Thru s  t  Line 
(in  XZ  plane) 


FIGURE  5-2 


AXIS  SYSTEM  FOR  LOCATING  THRUST  IN  AN  AIRFRAME 
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Assuming  the  thrust  line  lies  in  the  xz  plane,  then  the  thrust 

forces  can  be  resolved  in  this  generalized  manner. 

XT=Tcofl9T  (5-41) 

ZT=?sineT  (1.44) 

M_  =  T  (5-43) 

l  z 

For  a  steady  state  condition  the  equations  become 

X      =  To  cos  em  (5-44) 

0_.  T  4 

T 

Z      =  To  sin  OT  (5-45) 

°T 

!  1 

M      =  T  (5-46) 

o_        oz 
T 

while  in  a  disturbed  condition  where  the  body  axis  remains  fixed  in 

the  mass, 

XT  =  Tj  cos  0T  (5-47> 

ZT=  T1  sin0T  (i  ^"(5-48) 

MT=  T^  (    (5-49) 

and  T     =  To  +  A  T  '  (5-50) 

i! 

Assuming  the  air  density  remains  constant  so  the. engine  thrust  will 
not  be  affected 
' .'?.<-!£,    *+ts4£     >      **?"  (5-51, 

or 

T=Tu+Tw     .     ,         ^RPM 
u  g  RPM 
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These  are  the  considerations  for  the  external  forces  on  a  body 
in  motion.     While  the  dimensionless  coefficient  has  been  used  exten- 
sively by  the  aeronautical  engineering  field,  it  will  not  be  used  further 
in  this  work.     The  primary  mission  of  this  work  is  to  investigate  the 
control  response  and  paths  for  possible  compensation-.  , 

The  generalized  equations  of  motion  then  appear:  ( 

(5-52) 
m(U  +  QW  -  RV)  +  mg  sin  0     (cos  ©  cos.M*  )  -  mg  (cos  9     sin  0  ) 

>  . 

(cos  0  sin  V)  +  mg  (cos  ©    cos  0  )  (sin  0)  =  X    +  X'   u  +  Xs    w  + 

o  o  o  u  w 

X'  -w  +  X'    q+  X'nn.      +  X'JLO.       +  X'    v -i-  T     cos  ©^  +  T     cos  © _u  + 
H  H  vo  Tu  T 


T  >•  cos  ©  ,  r 

9  rpm  1  O 


o 


w  q 

co_ 
rpm  1  O  rpm 

(5-53) 

m(  V  +  RU  -  PW)  ;+  mg  sin  0     (cos  V  sin  0  sin  0  -  sin  V  cos  0)  -  mg 

(cos  ©     sin  0  )  (cos  ty  cos  0  +  sin  tysin  q  9{n  0)  .  mg  (cos  ©     cos  0  ) 

o  o  o  o 

(cos  ©  sin  0)  =  Y     +  Y'    v  +  Y'  .v  +  Y'   p  +  Y'    r       +  Y'*X      +  Y'  f   +  Y'f? 

I  (5>-54) 

m(W  +  PV  -  QU)  +  mg  sin  ©     (cos  Vsin  0  cos  0  +  sin    y      sin  0) 

-  mg  (cos  0   sin  0  )  (  sin  V  sin  ©  cos  0  t  cos    V    sin  0)  -  mg 

(cos  ©     cos  0  )  (cos  ©  cos  0)  =  Z     +  Z1  u  +  Z'    w  +  Z'.  W  +  Z'  q  + 
o  o  o  u  ww  q^ 

• 

Z'v/l     +  Z'lft        -  T     sin  ©^  -  T     sin  ©^u  T  T  w„„w    sin  ©„,  *»„„w 
■V  "X  o  T  u  T  £RPM  T^RPM 

(5-55) 

Af  +  aK.  +  QR  (C-B)  -  E  (PQ  +  R)  +cK,Q  -  bK,R  =  L    +  L'   v  + 
1  3  Z  o  v 
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(5-56) 
BQ  +  bK    +  PQ  (A-C)  +  E(P2  -  R2)  +  a  K  R  =  cK   P  =  M     +  M'    u 

M  X  ~t  ^  ** 

+  M'    w  +  M'  -w  +  M'    q  +  M'.q  +  M'^H,  +  M'n'X  +  T    z  +  T   u 
w  w  q  qn  1l  l  o  u 

+  T   w  £  rpme 

£  rpm     •    r 


(5-57) 


CR  +  cK,  PQ  (B-A)  +  E  (QR-P)  +  bK,P  =  aK,Q=  N    +  N1    v  + 

' , .     .  3  2  1  o  v 


N 


'    p  +  N'    r  +  N'£ J    +  INT^T   +  N'^r 


Thus,   it  has  been  shown  that  the  equations  of  motion  have  the  above 
external  forces  which  will  affect  the  rigid  mass  in  motion.     The  assump- 
tions made  are: 

1  .     The  earth  is  assumed  to  be  fixed  in  space  and  the  Earth's 
atmosphere  is  assumed  to  be  fixed  with  respect  to  the 
,    ;  Earth. 

2.  The  mass  is  a  rigid  b,ody. 
i      '  ■ 

3.  The  mass  of  the  body  is  assumed  constant  during  the  time 

of  the  problem.  i|  '' 

i    4.      For  an  aircraft  or  missile »  the  xz  plane  is  assumed  to  have  i 

mirror  symmetry. 
These  assumptions  were  made  in  the  first  section.     In  addition  the 
following  assumptions  and  simplifications  are  made  in  considering  the 
aerodynamic  force 8  :  i 
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1 .  The  flow  is  assumed  to  be  quasi- steady . 

2.  The  effects  of  auxiliary  aerodynamic  equipment  such  as 
flaps,   speed  brakes,   slats  and  spoilers  are  neglected.     lt 

3.  The  thrust  can  be  affected  only  by  change  of  speed  and/ or 
revolutions  per  minute.     No  Change  is  considered  for  change 
in  fuel  mixture,  manifold  pressure,   side  slip  or  propeller 
fin  effect.  * 

r 

These  assumptions  were  made  for  several  reasons.     In  the  first  and 
third  cases,  absence  of  adequate  information  on  many  of  the  coefficients 
makes  ignoring  or  setting  them  equal  to  zero  mandatory.     In  all  probability 
these  coefficients  are  very  small.     Assumption  two  covers  special  cases. 
While  these  items  may  have  considerable  effect  when  used,  they  are 
eliminated  from  the  general  case  at  this  time  because  an  aircraft  pro- 
bably uses  them  less  than  two  percent  Of  its   flying  time.    For  the  special 
cases,  these  effects  may  be  quickly  added. 

The  work  in  these  first  two  sections  is  developed  primarily  for 
background  purposes.     There  is  nothing  new  or  startling  in  this.     For 
further  examination  of  this  subject,  the  reader  is  referred  to  the  works 
in  Ref .   7  and  8  which  closely  parallel  the  development  in  this  work. 

The  aerodynamic  force  arid  moment  terms  carry  a  primed  notation 

for  ease  in  publication.     In  later  sections  these  terms  will  be  divided  by 

mass  or  moment  of  inertia.     These  terms  after  division  will  not  carry 

the  primeC  )  notation.     There  are  far  more  of  the  latter  terms  used.     Thus, 

it  is  for  ease  of  publication  only  and  does  not  affect  theory  in  any  way. 
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5.3      The  Airframe  Equations  in  a  Standard  Form 

This  section  will  take  the  general  equations  developed  in  the 
previous  chapters  and  simplify  them.     The  simplification  will  extend 
to  an  aircraft  in  straight  and  level  flight  and  subject  to  small  pertur- 
bations .     This  is  the  standard  method  of  dealing  with  the  problem  as 
found  in  present  day  texts.     However,   these  equations  will  be  used 
in  block  diagram  form  and  in  dete  rmmantal  arrays  as  developed  in 
Chu's  work.   (1)     The  results  will  be  compared  with  present  day 
classical  methods  to  determine  validity  of  this  determinantal  concept. 
The  first  simplification  is  the  concept  of  the  lateral  system  being 
independent  of  the  longitudinal  system.     Therefore,,   the  six  equations 
of  motion  are  divided  into  the  two  groups.     The  longitudinal  group 
includes  the  directional  motion  in  the  X  and  Z  directions  and  rotary 
motion  about  Y  axis.     The  lateral  group  consists  of  directional  motion 
along  the  Y  axis  and  rotary  motion  about  the  X  and  Z  axis.     The  axis 
system  referred  to  here  is  the  body  axis  system  of  the  airplane. 
Secondly,   the  aircraft  is  initially  in  a  steady  state  condition. 
This  is  the  reference  state  as  mentioned  previously.     When  the  dis- 
turbance terms  are  all  zero  Equation  5-52  reduces  to: 

m  (Uo  +  Qo  Wo  -  RoVo)  +  mg  sin  0o=  Xo  +  To  cos  0  (5-58) 

However,   it  was  stated  that  this  situation  was  restricted  to  longi- 
tudinal or  lateral  systems.     Thus  for  motion  along  the  X  body  axis, 
a  yaw  motion  Ro  and  a  side  motion  Vo  would  be  zero.     Further  it  was 
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stated  that  the  aircraft  was  in  a  steady  state  reference  condition. 


Therefore     U  would  be  zero.     Thus  Equation  5-58  reduces  to 

m  (Qo  Wo)  +  mg  sin  0     =  X'o  +  To  cos  0  T  (5-59) 

Similarly,   it  may  be  observed  that  Equation  5-53  for  the  steady 

state  reduces  to 

Y'o  =  0  (5-60) 

Likewise  Equation  5-54  reduces  to 

m  (Qo  Vo)   -  Z*     +  mg  cos  0     +   To  sin  0  T  =  0  (5-6l) 

o  o 

The  moment  equations  all  reduce  to  zero  provided  the  rotor  effects 
are  neglected  or  considered  zero. 

Lo  =  Mo  =  No  =  0  (5-62) 

Recalling  that  only  small  disturbances  were  to  be  considered 
and  these  disturbances  were  to  be  in  only  the  symmetrical  or 
longitudinal  system  for  the  longitudinal  group  the  derivative  due 
to  assymmetric  variables  of  V,   P,    R,   r  and   y  are  zero.     In  a 
similar  manner,   when  the  asymmetric  system  effects  are  considered 
the  symmetric  disturbances  are  ignored.     This  leaves  the  variables 
V,    W,   Q,   and  Tj  out  of  the  asymmetric  equations.     Finally,   because 
of  the  small  disturbances,   the  change  angle  0  will  be  small  so 
the  approximation  cos  ©  ==   1  and  sin  0  *=  ©  is  valid. 

The  result  of  the  previous  discussion  is  subtracting  Equations 
5-59  and  5-52,    eliminating  the  asymmetric  changes  and  supplying 
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t'he  small  disturbance  considerations.     The  result  is: 

m(u  +)  +mg   G  cos  6     =X  »U  +  X   »W  +  X.'W  +X  9Q  #r !    ,'    % 

ou  w  w  q  (5-63) 

+  X   'Y)      +  X-  »fl+Tu°  cos  G^U  +  Tr'    n„wcos  9 J^RPM 
*\   l  H   t  TO      RPM  T  ° 

For  the  z  direction  taking  5-61  from  5-54  and  using  the  same 
assumptions: 

m  (w-)  +  mg  e  sin  0     =Z'U+Z'W+Z.'W+ZIQ+  (5-64) 

o         u  w  w  q  v  ' 

* 

Z  yi  *fl   +   Z  •  !  H      -   TuU  sin  6m  =  T  r  ^ 

I    L  \     L  T  <5RPM   ORPM    sin  0T 

For  the  moment  in  the  y  axis,   the  following  applies: 

BQ  =MUJ  +  MwW  +  MwW  +  MqQ  +  MqQ  +  M^r[    +  M^/j  +      (5    65) 

T    'U    1   +  T  Cddu     6rPMz 
u      zl  5  RPM 

To  complete  the  case  for  the  symmetrical  group,   the  action  of  the 

elevators  must  be  included.     From  appendix  I  the  elevator  hinge 

moment  effects  are  reduced  to  : 

Ie  Y\^   =  H    'J     +  H    '.  T\   +  HewW  +  HewW  +  HeuU  +  HeqQ  +  £Fe 


In  most  cases  where  moment  of  inertia  terms  are  handled  it 
is  usually  expedient  to  use  an  axis  system  which  is  the  principal 
inertia  axis.     This  means  that  the  product  of  inertia  terms  are  zero 
For  aerodynamic  motion  problems,   this  applies  for  some  cases.  (8) 
This  case  is  not  one  of  them.     What  is  gained  in  simplification  of 
the  rotary  motion  equations  by  use  of  principal  axis  is  lost  in 
attempting  to  handle  the  force  equations  with  the  wind  velocity 
considerations  in  sine  and  cosine  terms.     By  orienting  the  X  axis 
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into  the  reference  wind  velocity  the  result  is: 

2  2  2  2 

V       =  U     +  v     +  w  (5-67) 

c 

where  U  =  U     +  u. 
o 

Therefore,   in  the  steady  state  reference  condition  V     =  U     and 

c  o 

V     =  W     =0.     With  the  assumption  that  u,   v  and  w  were  very 
o  o 

small,   their  products  and  squares  can  be  neglected. 


V     ~   \/  U   2  +  2U   u 


and  since  U     '  is  much  greater  than  2  U    u,   then  V     =     U       -  U  and  these 
o  o  c  o 

terms  are  interchangeable  in  this  section. 

As  previously  stated,   the  reference  situation  is  a  steady  state 

motion  and  includes  first  order  effects  only.      The  orientation  of  the 

axis  U     is  a  finite  value  while  V     =  W     =0.     Furthermore,   the  author 
o  o  o 

is  hard  pressed  to  visualize  an  aircraft  or  missile  tumbling  along  with 

a  Q     motion  although  P     exists  in  spinning  missiles  and  bullets       Furthe: 
o  e         o  r  o 

X      =  X     =  Z  •  =  Xf)  =  XY\     =  T  »*  =  0  as  an  assumption . 

w  q  w  't  L  (JRPM  ^ 

These  terms  are  usually  very  small. 

The  "standard"  determinant  form  as  stated  in  Chapter   3  of  Chu     is: 

i  +  pFj    pfh  pf2         pf21  pf3      pf31 

-PTj  1+      PF2  PF22  PF3  PF32 

-  PF  -      PF  1  +    PF  PF 

13  2  3  33 

where  the  PF*  indicates  a  performance  function       The  system  may  have 

fa    nodes.      This  example  contains  only  three  nodes.      The  "standard" 
block  diagram  for  this  form  is: 


*In  sections  2  and  3  the  symbol  G  was  used  to  designate  the  transfer  func- 
tion of  a  component,   and  the  component  itself  was  normally  a  distinct  piece 
of  hardware.     The  performance  function  PF  is  also  a  transfer  function, 
but  as  used  with  regard  to  the  airframe  it  designates  a  relationship  be- 
tween two  variables  which  are  not  visualized  as  the  input  and  output  of  a 
piece  of  hardware.      The  use  of  PF  notation  from  here  on  is  therefore  a 
deliberate  choice    to     emphasize  this  distinction. 
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FIGURE  5.  3    STANDARD  BLOCK  DIAGRAM  DRAWN  FROM  3x3 
DETERMINANT. 
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The  nodes  are  labeled  in  Roman  Numerals.     The  basic  rules 
of  nodal  signals  allow  several  inputs,   but  only  one  output.     This 
output  from  the  node  goes  directly  to  a  performance  function  with  a 
single  subscript  number.     Signals  may  be  picked  off  from  the  node 
output  and  fed  forward  to  another  node.     An  example  of  this  is  con- 
tained in  the  signal  going  through  PF      .     In  this  case  the  signal  is 
also  multiplied  by  the  function  PF        and  then  fed  to  node  III.     Feedback 
takes  place  in  a  similar  manner.     The  signal  output  from  the  first  system 
PF     is  picked  off  and  fed  through  PF       back  to  node  I.     Similarly  the 
other  nodes  function  in  the  same  manner. 

For  feedback  blocks,   the  first  subscript  number  denotes  the  out- 
put signal  source  for  the  feed  back  and  the  second  subscript  number 
the  input  node  to  which  the  feedback  signal  is  going.     The  notable 
points  about  the  determinant  are: 

1 .  The  main  diagonal  always  has  terms  of  one  plus  the 

direct  path  performance  function  times  the  self-loop 
feedback  performance  function. 
(1  +  PF1PFn) 

2.  All  terms  above  and  to  the  right  of  the  main  diagonal 
are  feedback  terms. 

3.  All  terms  below  and  to  the  left  of  the  main  diagonal  of 
the  "standard"  determinant  form  are  of  a  feed  forward 
nature.     All  of  these  feed  forward  terms  are  negative. 
The  term  immediately  below  the  main  diagonal  term  is  the 
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direct  path  performance  function.     Thus,  where  the  second 
column  main  diagonal  spot  has  1  +  PF     PF      ,   the  direct  path 
term  PF     will  appear  in  the  location  immediately  below  as 
-PF 
These  rules  are  strictly  mechanical.     For  a  rigorous  derivation  of  this 
system  the  reader  is  referred  directly  to  Chu's  work. 

In  Chapter  9  of  Etkin  (4),   there  exists  a  two  degree  of  freedom 
situation  where  ©  and  W  vary  and  the  speed  in  the  forward  direction  is 
constant.     This  is  the  short  period  approximation.     The  equations  using 
the  LaPlace  transforms  and  nondimensional  coefficients  are: 

{ZM  s  -  G        )  c*    -  (2^  +  Gz   )  s  0  -  Gz^^l      =  0  (5-68) 

-G         &       +  (i.  s2  =  G         )  ©  -  Gm^TI  =  0  (5-69) 

m<*.  p  mqs  'I  v 

2 

's  e 

This  author  also  grouped  his  acceleration  terms  such  that  G 


-Gh0C5       "V     e+(ies2"Ghfl)n=AC£  (5-70) 


\ 

included  C.      and  C,   .    .     Similar  situations  exist  for  G,  ,    G 

h*  h^  hoc  moC 

and  G        .     Note  that  r  =    ty    and  P  =  0  for  this  consideration.     From 
these  equations  a  block  diagram  is  evolved. 
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FIGURE  5.4    BLOCK  DIAGRAM  OF  TWO  DEGREES  OF  FREEDOM  IN 
THE  LONGITUDINAL  SYSTEM  SHOWN  BY  ETKIN4 
Further,   Etkin  derives: 


_e_ 

n 


U  (2^  s  s  Gzoc>  +  GzYL  mo( 


(5-71) 


G* 


Ao( 


(2/U  •  "  GzdC  )  (V  S      "   Gmq  )   "   Gm(2^+   °Zq)   S 

s 


G7(ifls     -  C  )  +  G      (  2^  +   G      ) 

Zvi     f-'  mq  m  'v  Zq 


jn 


(5-72) 


n,    -\ 


<2X  s  "   GZ«>  «6    *'  ~   °mq  >   "   Gr^2A+   GZq>   s 

s 


Now  it  is  intended  to  use  Chu's  work  involving  the  use  of  "standard" 
block  diagrams  and  "standard"  determinant  form  to  prove  that  the 
"standard"  system  applies  for  the  airframe  as  well  as  a  multiloop  servo 
system . 
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The  Z,  M  and  ^1  equations  of  motion  with  the  applied  assump- 


tions and  G     =  O  are: 
o 


w-Uq-Z     W  -    Z    Q  -    Z^  -   Z^  H.  =   0 
o  w  q  'L  v  <~ 

e-M    W-M.W-M    Q-  Mv\  TL-  Mvj 7l_=0 
w  w  q  'V  <• 

Y)-   lOy^U-   lOvi^-lO    W-10.W-10     q=Z\F< 
*-  I  l  w  w  q 

Recalling  that  UO  =  Uo  tan  /^0(  and  for  small  angles 


^v 


(s  -Z    )U    Acl  -  (sU     +  sZ  )e-n(sZ>i+  z«))1  = 
wo  o  q  L  L       >■ 


(5-73) 
(5-74) 
(5-75) 


IP=  U     ^CCand  0  =  q. 
The  equations  in  LaPlace  notation  then  become: 

(s     -   slO  „+  10  fl)>l  -  (sl0»  +   10    )U    ^<X-  (slOq)  0  =  A  Fe  (5-76) 

'  L  l       x  w  w        o 

(s2  -   sM    )  0  -  (sMyj  +  MY)   )  1\_  -  (sMcX  +  M  OC  )  Uo  Ak      =0  (5-77) 


(5-78) 


These  equations  in  standard  determinant  form  now  become: 


n 

A  OC 

e 

(sion     +  ion) 

2 

s 

(slO<X   + 

2 

s 

loot) 

slOq 

2 

s 

(sZfl       +    Zyj     ) 

Y           ZOC 

s(Uo  +   Zq) 

s 

s 

s 

(  sM  Yi       +  M  y)    ) 

(sMci  + 

MO() 

,         sMq 

2 
s 

2 

s 

1             2 
s 

i j 
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Ordinarily,   it  is  not  necessary  to  divide  by  s  or  s    „     However, 
it  is  done  here  to  demonstrate  more  closely  the  standard  form. 

To  take  the  differential  equations  and  arrange  them  by  rows  and 
columns  to  see  the  actual  physical  system  is  to  no  avail.    The  arrange- 
ment does  not  exist  in  the  above  form.     It  will  be  shown  how  to  achieve 
this  form  later. 

Consider  the  "standard  form"  array  as  set  forth  above.     The 
block  diagram  for  this  in  Chu's  "standard  form"  is: 


2 
s 

AFe 

2 

s       t  /* 

V 

_   -(SZ^+Z^) 

W^- 

-(sM^+Mtf) 

ys 

Y 

i 

e 

f      4>» 

M 

k 

\ 

h 

2 

s 

-(siotfripn) 

s/j 

=  Z 

-sM 

q 

2 

s 

< 

s 

-(SZ^H-Zy,) 

-(sMoc+MoO 

S           l 

c. 

s 

-MOflC+iObd 

-s(UQ+Zq) 

— 1 

s<- 

| 

-(sM^+Mo;) 
s^ 

s 

-slOq 

s2 

FIGURE  5.5    RESULT  OF  BLOCK  DIAGRAMING  EQUATIONS  5-76, 
5-77,   and  5-78. 
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The  above  block  diagram  then  is  a  three  node  system  but  the  outputs 
from  the  first  and  second  main  path  blocks  are  nothing  that  is  readily- 
recognized  in  an  airframe.     Because  of  the  fact  that  the  main  block 
following  the  third  node  is  1  ,   the  output  here  is  ©. 

A  few  simple  block  diagram  manipulations  clears  the  situation 
and  shows  true  signal  outputs  that  are  physically  realized  in  the  air- 
frame.    Block  diagram  manipulations  allow  the  changing  of  a  pick  off 
point  provided  the  transfer  function  is  corrected.     Further,   two  blocks 
can  be  combined  into  one  equivalent  transfer  function. 


-(sM^+MJ 
s2 

_-(sM^-Mr|) 

2 

s 

^ 

r(sZ^+Z^) 
s 

which 

-*0— 

1 

— 4k 

-(sZy^+Zj,) 
s 

i 

changes 
to 

A 

f 

-(slOfj+10/j) 

-(slOjj+lOfj) 

I 

S 

S  ° 

J- 

-(sZ^+Zfl) 

b         l 

Figure  5.6     EQUIVALENT  BLOCK  DIAGRAM  TERMS 
Adding  a  dummy  node  this  now  changes  to: 


Q-i-(sionnon 


-(sMri+M*) 

s^ 

(szAfZ^) 


r>"(sZfz^ 


Figure  5*7    FURTHER  BLOCK  DIAGRAM  COMBINATIONS 
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The  result  of  using  these  two  maneuvers  rearranges  the  block  diagram 
to  the  following: 


At? 


-(sM^+M^) 


(sZqf  Zy^) 
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^ 
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1-   sM, 
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l-sdt+Zq) 


-slO 
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Figure  5.8    REARRANGED  BLOCK  DIAGRAM  OF  FIGURE  5.5, 
USING  MANIPULATIONS  OF  FIGURES  5.6  and   5.7. 

Maneuvering  the  above  diagram  location,   the  block  diagram  takes  the 
form  of  Fig.    5.9,   which  coincides  with  Fig.    5.4. 
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Figure  5.9     TWO  DEGREES  OF  FREEDOM  IN  AN  AIRFRAME 


LONGITUDINAL  SYSTEM 


(rearranged  diagram  of  Fig.    5.8  to  compare  with  Fig     5.4) 
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It  is  noted  that  this  block  diagram  approaches  the  form  set 
down  in  Etkin  except  for  nondimensionalizing  the  terms.     The  dashed 
area  in  the  diagram  is  what  Etkin  labels  G_      and  Got*, 
The  determinant  from  the  above  dashed  portion  of  the  block  diagram 


is: 


1 

0 

0 

0 

0 

(sZrj    +   Zyj) 

s 

0 

-s(Uo  +   Zq) 
s 

0 

0 

-   1 
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-  (sMO(+Mo() 
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0               1 

-1 

1 

0  Aq( 

Evaluating  for  the  transfer  functions -r^-    and     ~     ,   the  results 


will  be  in  the  form 


Gjk  A 


n 


\ 


1Y 


This  form  is  the  type  expressed  by  Chu  (1)  and  Ihaler  (12),     The  /\ 
term,   which  is  the  whole  determinant,   reduces  to: 


z\  = 


i  -  z<* 

s 

-1 
0 


•  (sM°C    +  MOC) 


s(U     +   Zq) 


0 

1   -   sM 


s 


which  is  A  =  (1   -   ZoC)    (1   -   sM      )  -   s  (U     +   Z   )  (SMOC      +  MC*)        (5-79) 

.  q  °       q 


or  A  =   1     (s  -   Zrf  )  (s     -   sM   )   -   s  (Uo  +   Z   )  (sM<X  +  MO(  ) 

~  q  q 

s 
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To  consider  the  — — —  term  the  cofactor  Aiy  of  this  determinant  is 
obtained  by  crossing  through  the  first  row  and  the  second  column. 

-(sZn+zn>  0 


-A 


12  = 


s 
0 


(sM  t|     +  My|    ) 

2 

s 


s  (U     +  Z   ) 

o q 


s 
0 


(sM  <i     +  M  <*    )  1  -   sM 


A 


12  =  +1 
3 


+  (sZ  Y)  +   Z  rj  )  (s     -   sM   )  +  s  (U     +   Z   )  (sM^+M-n) 


The  main  feed  forward  transfer  function  between  the  input  signal  from 
node  2  to  the  output  is  +1  in  this  case.     This  is  also  true  of  node  4 
input  to  output.     Hence  G„   A  ^       =  +   1. 
The  transfer  function 


'2A0( 


'AOCYl 


Adi 


n 


=  c^a 


A~A 


12 


A 


(5-81) 


GA«r\ 


(sZY[        +  Zy\)(s     -sM)+s(Uo+Z)(sMr]_+Mr|) 

2  • 

(s  -   Zct)     (s      -   sM    )  -   s(U     +   Z    )  (sMOC  +  M  CX  ) 

q  o  q 


It  can  be  noted  that  this  equation  5-80  is  the  same  equation  as 
5-72  given  by  Etkin  save  the  non-dimensionalizing  terms. 
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To  evaluate  the  transfer  function  9  1  yi  ,   the  same  method  is  used 


1 


G        =  -   1 
40 


A 


14 


(sZy)+   Z-n)  1   -   ZoL 


■V)  +   ^T| 


-1 


(sM>i  +  Myi) 


-(sMOC  +  MOC)        0 


14  = 


(sZ  yi     +  Zyj  )  1  -  ZQC 


-(sMyi      +  Mfl  )  0  -(sMOl  +  Ma  ) 


^14  =J_       (sMY|_+Mn    )(s-   ZQC)  +  (   sM<k+  M(X)  (sZV|+   Zy|) 


(5-82) 


Therefore: 
r        _     © 


°4e   A4 


G  _(sMr\     +  MjQJ  (s  -   Zgl  )  +  (sM(X+  MqQ  (sZy^+  Z^) 

^l         (s  -   ZCX)  (s2-   sM    )  -  s(U     +   Z   )  (sMOC+  McX) 

q  o         q 


(5-83) 
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Here  again,   excepting  the  non-dimensionalzing,   the  result  is  the 
same.        Equation  5-83  is  equal  to  equation  5-71  as  given  by  Etkin. 

The  system  of  "standard"  forms  then  works  for  the  airframe 
differential  equations  after  some  manipulation.     It  can  be  noted  that 
there  is  a  method  of  making  up  this  standard  form  without  the  necessity 
of  the  manipulation  which  will  yield  physically  realizable  output  signals. 
When  the  block  diagram  of  Fig.   5.8  is  converted  into  a  determinant 
for  solution  of  the  transfer  function,   there  appears  a  series  of  alterna- 
ting rows  where  the  equation  -1  +   1  =  0  appears.     It  also  appears  that 
these  simple  equations  occur  in  a  sequence  where  the  plus  one  term  is 
placed  in  the  2n  !    row  and  column,     n  is  the  number  of  differential 
equations  in  the  problem.     The  minus  one  term  is  placed  in  the  2n  ! 
rows  and  (2n  !-l)  columns.     All  feedforward  terms  (those  below  and 
to  the  left  of  the  main  diagonal)  are  moved  one  column  to  the  right 
and  remain  in  the  same  row.     This  doubles  the  size  of  the  determinant. 
However,   there  are  several  zero  terms  introduced  and  the  actual  work 
involved  in  evaluation  of  the  determinant  is  not  markedly  increased. 

The  next  item  is  to  look  at  the  lateral  motion  problem.     It  is 
slightly  more  complex  because  there  are  five  differential  equations 
of  motion  instead  of  four. 
Equation  5-35  with  simplifications  becomes: 

(v  +  Uor)   -  0  mg  cos  0     -   YvV  -   Yp  -   Y     R  -  YrS        =  0  (5-84) 

o  r  R  S 
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9 

(/£>+  r)  -  0  g  cos©         -   Ygj3  -  Y     p  -  Y^r  -  Y  £.S       =0  (5-85) 


u  u         u         u 

o  o  o  o 


The  moment  equation  5-55  becomes 


P  -_E    R  -  1->fi/8  -Lp-Lr-Lv.^       -   L  ~  9  =0  (5-86) 

1\.  J 


and  equation  5-57  becomes 

EP    -  N^/#   -N   P  -  N   R    -  Nv'£      -  N£,£      _N$^        =0  (5-87) 

The  aileron  hinge  motion  equations  from  Appendix  I  are 
^  -   11    p  -   11    r  -   llv£     -H§&      =    ^Fa  (5-88) 


•  « 

R  +  EP 


And  also  from  Appendix  I,   the  rudder  hinge  motion  is 

••  .* 

?  -\Z g/$-   12   p  -   12    r  -   12q§    -   12^5     =    AFr  (5-89) 

From  these  terms  a  determmantal  array  can  be  made.     It  should 

be  noted  that  this  time  a  set  of  equations  of  - 1  +   1  =  0  will  be  inserted 

in  the  system.     This  permits  direct  drawing  of  the  block  diagram  and 

omitting  the  need  of  manipulation.     Table  5-2  then  presents  this  array. 

From  this  Fig.    5.10  is  drawn  directly. 
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Figure  5-10  Block  Diagram  of  Airframe 
Lateral  Stability  System   as  Drawn 
Directly   From   Table   5-1 
(Roman  Numerals   Indicate  node  number) 


qfun-  .,.,   to- 


FIGURE  6-1       Block  Diagram  of  the   Airframe  ond  control  System 
Using  the   Standard  Form 


The  resulting  block  diagram  very  clearly  parallels  that  of  the 
«:eral  system  shown  by  Etkm     in  Chapter  9  (4).     The  "standard"  form 
:( achieve  this  block  diagram  differs 

There  is,   then,   the  question  of  why  does  this  differ  from  the 
jiginal  "standard"  form  set  forth  by  Chu. 

Chapter  Two  of  Chu  (1)  discusses  at  some  length  the  idea  of 
irangment  of  varying  performance  functions  into  a  "standard"  block 
digram.     The  individual  pieces  of  hardware  are  at  hand  with  known 
[dividual  performance  functions.     The  orderly  arrangement  of  these 
srmit  the  construction  of  a  "standard"  determinant  for  the  whole  system 
From  this,   a  multiloop  servo  system  can  be  handled  in  generalized 
r  suits .     In  the  case  of  the  airframe,   the  determinant  formed  by  the 
ifferential  equations  was  known  and  the  block  diagram  had  to  be 
rmed.     Admittedly  a  block  diagram  (Fig.   5.5)  was  made  up  but  the 
:itput  signals  are  not  the  physically  realizable  type  that  are  generally 
lown  to  pilots  and  engineers.     The  same  performance  function  results 
torn  the  system.     Both  are  mathematically  correct. 

The  difference  lies  in  the  fact  that  Chu  is  taking  individual  pieces 
'  hardware  where  the  functions  are  known.     The  blocks  are  made  up 
id  then  the  "standard"  determinant  form  is  made.     In  the  case  of  the 
trplane,   the  equations  of  motion  are  known  and  a  determinant  is 
nmediately  formed.     Each  column  contains  the  performance  functions 
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multiplied  by  the  direct  function  signal  output  term,     Chu's  standard 
form  dictates  that  all  terms  above  the  main  diagonal  are  multiplied  by 
the  node  output  or  input  to  the  direct  function.     Steps  must  be  taken 
then  to  change  this  in  the  differential  equations.     The  addition  of  the 
alternating  -1  +   1  =  0  equations  solves  this  very  nicely  and  in  no  way 
alters  the  value  of  the  equation, 

A  series  of  evaluations  of  determinant  forms  is  given  in  Appendix 
III.     The  results  show  that  the  addition  of  these  terms  in  no  way  changes 
the  final  value  of  a  determinant  or  a  cofactor  value. 

Referring  to  Fig.    5.8  there  are  other  block  forms  that  may  change 
the  values  within   the  blocks,   but  will  not  alter  the  final  outcome.     Con- 
sider the  first  node  in  Fig.   5.8  as  an  example.     The  four  forms  expressed 
in  Fig.    5.11  are  equivalent.      The  original  form  is  shown  in  Fig.    5.11a. 

The  self- loop  feedback  path  is  combined  with  the  direct  path  in  Fig. 

2 
5.  lib.     The  s     term  is  eliminated  and  the  performance  functions  are 

then  in  the  form  shown  in  Fig.    5.11c.     A  further  change  can  be  made  by 

dividing  the  nodal  input  terms  by  the  denominator  of  the  direct  path 

feed-forward  function.     The  result  is  then  in  the  form  shown  in  Fig. 

5.  lid.      The  performance  functions  in  each  block  differ  somewhat,   but 

the  final  result  remains  the  same.     These  manipulations  allow  the 

user  more  freedom  in  the  forms  of  the  performance  functions.     It 

allows  changes  in  form  to  fit  the  needs  of  the  particular  problem. 
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Figure  5.11  (a) 


Figure  5.11  (b) 
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Figure  5.  11  (c) 


Figure  5,11  (d) 


FIGURE  5..  11     PERMISSIBLE  BLOCK  VALUE  VARIATIONS  AT  A 
NODE  . 
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The  pick-off  point  for  the  feed  forward  terms  does  differ  in 
requirements.     It  is  seen  in  Figs.    5.11a  and  (d)  that  the  pick-off  can 
be  either  off  the  output  signal  of  node  I  or  II.     The  final  result  will 
not  change  .     However,   this  is  not  the  case  in  Figs  .    5.  lib  and  (c).     The 
pick- off  point  must  be  off  the  output  signal  of  node  II.     The  final 
result  will  differ  if  the  output  of  node  I  was  tapped  in  the  latter  case. 

For  converting  a  set  of  linear  differential  equations  into  the 
standard  block  diagram  form,   the  following  rules  can  be  applied. 
1  .     Between  each  differential  equation,   an  equation  -1  +   1  =  0 
shall  be  inserted.     An  additional  equation  -1  +   1  =  0  shall 
be  added  at  the  bottom  of  the  set  so  the  number  of  additional 
equations  shall    equal  the  number  of  differential  equations 
originally  in  the  problem. 

2.  The  arrangement  of  the  values  of  these  added  equations  in 
the  determinant  shall  be  such  that  the  positive  one  value 

is  always  in  the  main  diagonal  block.     The  minus  one  term 
shall  always  be  placed  in  the  block  immediately  to  the  left 
of  the  main  diagonal  term. 

3.  All  other  blocks  in  the  rows  of  these  added  equations  shall 
be  of  zero  value. 

4.  All  finite  performance  functions  of  the  original  differential 

equations  which  are  located  below  the  main  diagonal  shall 

be  moved  one  place  to  the  right.     The  vacated  blocks  shall 

be  replaced  with  a  zero  value,      This  last  rule  is  optional 
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to  the  user.     It  does  save  the  possibility  of  error  which 
could  appear  because  of  cases  shown  in  Fig .    5.11b  and  (d)  . 
With  the  simplified  rules  of  making  a  standard  block  diagram 
from  the  standard  determinant,   a  valuable  tool  is  developed  for  many- 
purposes  .     It  can  provide  a  student  with  a  means  to  draw  and  visualize 
a  complex  feedback  system  such  as  an  airframe  from  a  group  of  differen- 
tial equations.     For  the  research  groups,   it  makes  analog  computer 
programming  a  simplified  matter  with  respect  to  amplifier  usage.     In 
this  latter  case,   pitfalls  do  exist.     The  airframe  is  no  simple  matter. 
The  reader  is  referred  to  Howe  (8)  for  an  excellent  coverage  of  the  air- 
frame problem  for  analog  computers.     Most  important,   though,   is  the 
fact  that  this  standard  form  of  Chu's  is  valid  for  any  multi-loop  linear 
feedback  system  whether  it  is  electrical,   mechanical  or  a  combination 
of  the  two.     Thistaenputs  a  regular  form  to  the  multi-loop  control 
system  from  which  design  or  synthesis  can  proceed  in  an  orderly  manner. 
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6.       Expansion  of  the  Standard  Form  for  Greater  Disturbances 

The  system  developed  by  Chu  was  for  linear  functions.     It  is 
the  purpose  of  this  section  to  explore  the  equations  of  motion  and 
determine  whether  the  equations  of  motion  can  be  used  for  larger  dis- 
turbances than  considered  in  Section  V.     If  so,   then  to  what  degree 
can  this  system  be  used. 

For  all  practical  purposes,   the  small  angle  approximation  holds 
quite  well  up  to  one-tenth  radian.     Therefore,   G^r     sin  0  and 
cos  0^1,     With  this  simplification  in  mind,   Equation  5-52  is  then 
rewritten. 
m  (  U  +  QW  -  RV)  +  mg  sin  0     -    lf>  mg  (cos  0     sin  0  )  + 

0  mg  (cos  0     cos  0   )  =  X     +  X»  U  +  X  !W  +  X1.  W  +  X'    Q  +  X'   Y\ 
°  o  oou  w  w  q  ,^N* 

+  X^.         t  X'vV  +  Tq  cos  ST  +   Tu  cos  QTU  +   T&RpM  £RPM  cos  @T 

(6-1) 

The  thrust  forces  are  combined  into  the  forward  velocity  motion 

u  and  a  steady  state  reference  condition  of 

m  (Q   W     -  R     V   )  +  mg  sin  0     =  X     +  T     cos  ©^  (6-2) 

oooo  ooo  T 

is  assumed. 

Subtracting  6-2  from  6-1  the  following  equation  results 

m(u+Q        +W       +qw-R        -Vr-vr) 
ow  oq  ov  o 

-   Vmg  (cos  0     sin  0   )  +  0mg  (cos  0     cos  0      )   -  X' 
5  v  o  o  o  o  u 


X' 
w 


X'.W  -  X'Q    -  XUY)      -  XfA    -  X'V    =  0  (6-3) 

w  q  ^j  V.  r[  v 
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The  product  term  qw  is  dropped  as  a  small  value  product  situation. 
The  same  applies  to  the  vr  term, 

Converting  to  La  Place  notation  and  dividing  by  the  mass  the 
following  equation  will  result 

+  e  u  a&+  sw  e  -  w   u    A  -  IP  g  cos  e    sin  0 

00  o  ^  o      o  r  o  o 


su 


+  0  g  cos  0     cos  0     -  X   U  -  (sX  .  +  X    )  Zl  &    -   sXq  e 
o  o  u  ww 

-  (sX^+X^Y^-  Xgft     =0  (6-4) 


The  y  and  z  direction  force  equations  will  form  up  in  the  same  manner. 

« 

sU.+^u+U      l?     -  0     U      A<3.  -   W     0  -    tyg  sin  © 
QA  '  o  o  o      o  o  o 

-  g  (cos  ©     sin  0   )  -  g  0  cos  ©     cos  0     +  Yg/S  (6-5) 

o  °  ' 

-sY    0  -   sYr  (p.  (Y>    )£    -  (sY^    +     Y^)1?     =  0 

sU    A&  +  P   U     /?    +  V        0  -  Q        -V       ©  4-  g  sin  ©     © 

O  OOr-^  OS  OUOS  O 

+  g  cos  ©     sin  0       0  -   g  cos  ©     cos  0     -   Z   U 
o  o  o  0  u 

-(sZw+Zw)U     -AO^-sZ     ©  -  (sZ  Yj  +  Zy|)  H  =  0  (6-6) 
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The  moment  equations  can  be  formed  in  much  the  same  manner.     From 
Section  Five,  Equation  5-55  would  convert  to  this  equation: 

s2  o  +  (£i5. )   (so    \j>  +  sr    e  )  -  4-  (sp    e  -  sQ  ty  )  -  -4-  s2  ip 

A  o      '  o  A  o  o  A 

(6-7) 
-  (sL^    +  L^  )^<g-   L   g  0  -  (sLr  7J))  -  (sL  J       +L$     ){+lJ 

In  a  similar  manner  Equations  5-56  and  5-57  can  be  brought  to  a  similar 
maximum  linear  condition.     They  are  shown  in  the  following  equations: 


s    ©    +  {-—-^  )  (sP     ©  +  sQ   0)    +  —=-     (sP     0  -   sRo    lf>  ) 
Boo  ^^  o 


-  M   u  -  ( sMw  +M)U      ^   cC  -   s  ( sM  •    +  M        )      Q      -  ( sM  -y% 
u  wo  q  q      '  v 


M^)7[    =o  (6-8) 


s2    IP    +  (-^=A)     (sP   ©      +  sQ     0  )  +  -4~     (sQ     XV  +  sR     0  -   s20    ) 
*  C  o  o  C  o   ^  o 

-  N^  -sN     0    -   sNR  V    -     Nfc£        -  (sNw       +  N^  )?  (6-9) 

Also  the  control  hinge  moment  Equations  5-66,    5-88,   and  5-89 
are  recalled  in  their  same  form. 
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These  nine  equations  represent  the  linear  equations  of  motion 
of  the  airframe  and  include  all  cross-coupling  effects       These  equations 
still  retain  the  requirement  of  a  linear  approximation       The  equations 
with  the  author's  modification  of  adding  the  equations  -1+1   =  0  make 
up  Table  6-1  .     This  is  an  18  x  18  determinant       From  this  table  the 
block  diagram  of  Fig.    6.  1   can  be  drawn. 
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This  block  diagram  then  represents  the  complete  airframe 
response  maintaining  linear  approximations.     Even  then  this  system 
has  limitations,     As  an  example,    consider  some  airplane  flying  with 
a  steady  state  wing  angle  of  attack  of  15  degrees,     (assume  a  "con- 
vention" wing  and  the  reference  axis  is  the  zero  lift  line)    Any  pilot 
will  quickly  testify  that  an  increase  of  5  degrees  angle  of  attack  at 
that  point  will  probably  result  in  a  rather  wildly  oscillating  maneuver 
popularly  known  as  a  spin       This  condition  is  generally  outside  of  the 
range  of  the  aerodynamic  forces  and  moments  usually  computed. 

Some  aerodynamic  coefficients  such  as  lift  (C    )  are  decidedly  non- 

JLi 

linear  in  this  region,.  Therefore,  the  restrictions  imposed  on  this 
derivation  of  terms  must  include  a  restriction  that  the  airframe  is 
not  operating  at    the  very  edge  of  its  flight  envelope. 

The  system  has  been  diagrammed  directly  from  Table  6-1  .     The 
terms  in  Column  One  are  all  connected  at  node  I  of  Fig.    6,  1  .     Likewise, 
the  terms  in  Column  Two  of  Table  6-1  are  all  connected  to  node  II. 
Extra  nodes  appear  in  the  beginning  and  end  of  the  block  diagram. 
These  are  inserted  as  provisions  for  control  functions.     These  control 
nodes  are  lettered  alphabetically.     The  block  diagram  contains  the 
required  loops  for  aileron  and  elevator  control.     The  rudder  control 
loop  can  be  easily  inserted  when  desired. 

Examination  of  Fig.    6.  1   reveals  many  interesting  points  here- 
tofore known  but  sometimes  only  by  intuitive  reasoning.     The  concept 
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of  separate  systems  for  the  symmetric  axis  and  asymmetric  axis 
systems  is  very  well  shown.     The  level  flight  steady  state  conditions 
used  in  Section  Three  stated  that 

Po  =  Qo  =  Ro  =  Vo  =  Wo  =  0  (6-10) 

Examining  the  paths  between  nodes  VIII  and  IX  reveal  that  all 
of  the  performance  functions  are  zero  under  these  specified  conditions 
of  steady  state  flight.     There  is  no  direct,   feedforward  or  feedback 
paths  between  nodes  VIII  and  IX.     Thus  the  restricted  case  of  steady 
state  level  flight  conditions  with  small  perturbation  is  clearly  divided 
into  the  two  systems.     The  longitudinal  system  includes  the  forward 
velocity,   vertical  velocity,   pitching  motion  and  elevator  position. 
The  block  diagram  for  the  system  can  be  further  restricted  to  two 
degrees  of  freedom  by  holding  the  X  velocity  change,   u,   at  zero. 
This  allows  change  in  vertical  motion  and  pitch.     This  would  then  be 
Figs.   5-8  and  5-9  of  Section  5. 

The  so-called  lateral  airframe  response  system  comprises  the 
nodes  from  IX  to  XVIII.     The  rudder  and  aileron  motion,    sideslip, 
yaw  and  roll  comprise  the  motion  for  this  system.     With  the  steady 
state  level  flight  conditions  previously  mentioned,   the  block  diagram 
is  identical  to  Fig.    5-10  of  Section  5. 

An  interesting  note  of  these  separate  systems  is  the  cross 
coupling  terms  that  begin  to  arise  with  larger  disturbances.     An 
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important  one  is  the  pitch  due  to  sideslip.     The  effects  are  labeled 

M>^  in  the  diagram.     Several  aircraft  are  subject  to  this.     Usually, 

this  is  a  tucking  or  nose  down  motion  that  appears  as  sideslip  angle 

(  /p  )     increases.     In  the  simplified  systems  this  effect  was  neglected, 

It  is  of  negligible  value  for  very  small  disturbances,   but  can  become 

noticeable  by  the  time  the  sideslip  angle  (    £j  )     reaches  five  degrees. 

The  effect  exists  and  is  shown  by  this  expanded  system 

Excellent  methods  of  analysis  are  available  for  these  standard 

forms  brought  forth  in  this  work.     Methods  of  analyzing  outputs  in 

response  to  one  or  more  inputs  are  available.     Work  by  Thaler  (10) 

contains  a  very  useable  system. 

Where         N  =    the  output  signal  at  the  n  th  node 

A   =     characteristic  determinant 

A  =     the  cofactor  of  the  determinant  with  "a" 

n  the  node  location  receiving  the  input 

and  "n"  the  node  location  of  the  output 
signal. 

A 

G       =     the  performance  function  on  the  input 

1  9. 

signal  prior  to  entering  the  input  node. 

I        =     the  input  signal 

G        =     the  direct  path  performance  function 
between  the  output  node  and  the  output 
signal. 

Thus  a  single  input  signal  would  yield  an  output  signal 


NT  =  G  . _ 
I  nN 


A       G,         I 
a  1 

n  a 
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(6-11) 


For  two  or  more  input  signals  additional  subscripts  would  be  used 
and  the  functions  are  additive.     The  rules  of  superposition  of  linear 
outputs  justify  this  method. 


N  =  G 


nN 

A 


A         G,  I  +  A.      G_     II  +  A      G.     III  + 

a  1  b        2n  c        3 

n  a  n        b  n        c 


(6-12) 


This  equation  would  give  the  resulting  ouput  of  one  or  more  inputs 
at  various  nodes  in  the  system. 

The  degrees  of  freedom  can  be  restricted  also.     In  the  block 
diagram,   (Fig.   6-2)  consider  a  cut  made  in  the  path  at  point  "s". 


1K> 


FIGURE  6-2    ILLUSTRATION  OF  RESTRICTING  FREEDOM  OF 
MOVEMENT  IN  A  SYSTEM  BY  CUTTING  A  NODE  OUTPUT  IN 


THE  BLOCK  DIAGRAM. 
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This  opening  in  the  path  now  eliminates  the  signal  to  the 
direct  path,   the  feedback  and  the  feedforward  path.     The  result 
of  this  in  the  determinant  is  the  insertion  of  zeros  in  all  the 
performance  functions  in  the  column.     It  leaves  only  the  constant, 
one,   in  the  main  diagonal  block.     The  reader  will  recall  that 
the  main  diagonal  has  the  term  1  +  PF   PF       .     This  is  shown 
below  as  opening  the  output  path  from  a  third  node. 

x  x  0  x  x  x 

x  x  0  x  x  x 

x  x  1  x  x  x 

x  x  0  x  x  x 

x  x  0  x  x  x 

x  x  0  x  x  x 


The  output  path  from  node  three  is  opened  at  point  s  and  the  deter- 
minant is  changed  in  this  manner. 

A  special  case  occurs  for  the  diagram  illustrated  in  Fig,    6r2. 
The  second  node  has  only  one  input  signal.     Because  the  path  was  cut 
at  point  "s",   the  output  from  the  second  node  will  also  be  zero.     Thus, 
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all  feedforward  values  from  this  node  will  be  zero.     The  direct 
and  feedback  paths  of  the  second  node  are  non-existant  to  being  with 
and  remain  so.     The  determinant  form  would  appear  as  now  shown. 

x      x      0      0      x      x 

x      x      0      0      x      x 

x      x      1       0      x      x 

x      x      0       1      x      x 

x      x      0      0      x      x 

x      x      0      0      x      x 

Simple  manipulation  now  reduces  the  six  by  six  determinant  to  a 
four  by  four  determinant.     Therefore,   when  it  is  desired  to  consider 
a  problem  with  one  or  more  degrees  of  freedom  held  at  zero,   the 
simple  method  of  opening  the  output  signal  path  from  the  node  and 
placing  a  zero  for  all  the  performance  functions  in  the  column  leaving 
only  the  constant  value  one  in  the  main  diagonal  point  accomplishes 
this  . 

Referring  to  Section  5.    the  reader  can  observe  the  general 
equations  of  motion  for  mass  particles  and  examine  one  such  as 
Equation  5-55  which  is  repeated  below. 
AP  +  a   K     +  QR  (C-B)   -  E  (PQ  +  R)  +  cK   Q  -  bK^R  = 

L     +  L'v    +   L'.v  +  L'    p  +  L'    r  +  L'tT       +  L'  >  F  (6-13) 

ov  v  p^  r  rJ  rJ 

+   L'  V      +   L'lV    -  0 
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Note  the  third  term  QR  is  a  product  of  (Q  +  q)  (R     +   r) .     The  pro- 
duct when  multiplied  isQR     +Qr+Rq+qr.     From  the  equations 

o    o  o  o 

of  Section  5: 

Q  =  0  cos  0  +    ^  sin  0  cos  0  (6-14) 

R  =  (/'cos  0  cos  e  -  0  sin  0  (6-15) 

When  the  small  angle  approximation  is  no  longer  valid,   the  handling 
of  this  product  term  becomes  rather  involved.     Furthermore,   the  aero- 
dynamic coefficients  of  the  L1     ,    L!     ,    L'      etc  terms  lose  their  linear 

v  p  r 

approximations  over  large  changes  in  flight  conditions.     Most  aero- 
dynamic coefficients  are  of  a  decidedly  non-linear  nature  over  the 
full  span  of  an  aircraft's  maneuvering  capability.     The  coefficients 
are  definitely  non-linear  with  Mach  number  when  the  velocity  range 
covers  the  subsonic,   transonic  and  supersonic  spectrum.     The  reader 
can  observe  the  general  equations  of  motion  in  Section  5.     Considering 
the  complexities  of  incorporating  equations  6-14  and  6-15  for  large 
movements  and  the  non  linearities  of  the  aerodynamic  coefficients,, 
it  is  not  long  before  one  realizes  that  the  non  linearities  are  numerous 
and  cover  all  the  equations.     The  prospect  of  block  diagramming  for 
the  purpose  of  placing  all  non  linearities  in  one  block  to  be  handled 
by  a  describing  function  scheme  is  hopeless. 

The  possibility  of  doing  the  problem  with  an  analog  computer 
and  using  no  approximations  is  fair.     Undoubtedly  the  task,   if  attempted. 
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will  be  of  considerable  magnitude. 

Thus,   it  is  seen  that  this  system  of  standard  determinant  arrays 
and  standard  block  diagrams  can  be  used  to  handle  the  airf rameproblem 
as  long  as  a  linear  approximation  is  maintained.     The  system  is 
flexible  and  can  carry  varied  conditions  of  multiple  inputs.     It  can 
also  be  noted  that  for  the  six  degrees  of  freedom  problem  over  a  large 
variation,   it  is  of  limited  value  except  for  diagramming  a  computer 
problem . 
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7.     APPLICATION  OF  THE  DETERMINANTAL  APPROACH  TO 
A  SIMPLIFIED  AIRFRAME  CONTROL  PROBLEMS. 

The  system  to  be  considered  is  a  surface  to  air  homing  missile 

using  proportional  navigation,     A  sketch  of  the  missile  shape  is 

shown  in  Fig.    7.1,   which  also  defines  the  coordinate  system.     A 

guidance  computer  provides  command  signals  for  the  yaw  and 

pitch  autopilots,  which  have  the  same  form  because  of  the  symmetry 

of  the  missile.     The  missile  is  to  be  considered  a  rigid  body,   and 

gravity  forces  are  considered  to  be  trimmed  out  in  steady  flight. 

The  missile  is  also  to  be  considered  in  steady  state  operation 

such  that  the  speed  (Mach)  is  constant  and  the  pitch  angle  of  attack 


is  const 


tant  at  C(    .     Assuming  that  </?>  —  5;  cos  /zp  - 1.0,    sin  /?  -   jCs  , 
assuming  that  the  mass  and  moments  of  inertia  are  constant,   the 
linearized  equations  of  motion  of  the  missile  can  be  derived  and 


are: 


Yaw  force  equation 

*  '  (7-1) 

3     -  0  sin    C<    +      lp  cos    CX         -  Y^#  -  Y^    _  Yt)    =0 

Yaw  moment  equation 

•* 

^  -N.^  ~N^5  -Ny£  =0  (7-2) 


N^/^  -N^5  ~N>X 


Roll  moment  equation 


0  -   Lp  0  -   L^/3  ~L>t  "L%5^    =  °  (7"3) 

Where  the  coefficients  are  as  defined  in  Section  4„ 


121 


Fig.    7-1      SKETCH    OF  MISSILE  SHAPE 
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Since  pitch  is  constant,   the  equations  for  the  roll  and  yaw 

control  systems  are  needed.     The  measuring  devices  used  are 

accelerometers  which  measure  acceleration  in  the  y     and  z 

direction  or,   within  a  scale  factor 

Y^/4  +  Y£$      +VJ  Y^£       andY^OC+     Y  *£     +      Ytj 

Other  measuring  devices  are  rate  gyros  which  measure  0,    4>    and  ©. 

The  relationships  existing  in  the  controlled  missile  are  defined 

by  the  functional  block  diagram  of  Fig.    7.2. 

If  the  sensors  are  considered  perfect  G,   =G„  =G  =  G  =G  =   1.0, 

12         3       4       5 

and  the  control  equations  may  be  written: 

SJ      =     [K1Y/+K1     YJ$+K1Y55]       +K2^+K3^     (7"4) 
sj      =  (-)        I  K4  0  +  K5  0+  K6  0  1  (7-5) 

The  coefficients  of  these  equations  may  be  arranged  in  a  deter- 
minantal  array  as  shorn  on  Table  7-1. 

A  block  diagram  may  be  constructed  directly  from  this     jL\  , 
then  manipulated  using  the  methods  indicated  in  Section  5,   or  some 
preliminary  algebraic  manipulations  may  be  used  to  define  transfer 
functions.     Using  the  latter  scheme,   each  row  is  divided  by  the 
entry  in  the  main  diagonal  of  that  row,   providing  a  new  determinantal 
array    (Equation  7-7)  shown  as  Table  7-2  and  7 -2a. 

The  standard  block  diagram  may  be  drawn  from  Equation  7-7, 
and  inspection  of  the  diagram  indicates  that  a  rearrangement  of 
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FIGURE    7-2 
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Where  the  performance  functions  of  Section  5  are  defined  as: 


A 


TABLE    7 -2a 


$ 

(K2    s   +KJ 
(a  -  KjY*,    ) 

Ki  V 

* 

5 

(s-Y^) 

-cosof        r 
(s-Y^)          <2> 

sin    o( 
o 

0 

L? 

5 

X 

(■ 

(s-Lp) 
-K6sZ+K5sK4) 

fi 
$ 

(s-KjY„     ) 

KiYi 

2 

s 

] 
V 

(s-K,Yh     ) 
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(s-Y^  ) 

Yi 

$ 

(s-Y^) 

£_         _^_  _f_    ±s_ 


Ji  (s-Lp) 
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nodes  and  blocks  can  produce  a  separation  of  the  roll  and  yaw 
subsystems.     If  this  is  done  the  block  diagram  takes  the  standard 
form  shown  in  Fig.    7.3.     Note  that  the  yaw  and  roll  subsystems 
are  clearly  separated,   and  that  specific  blocks  are  obviously 
representative  of  cross-coupling  effects.     A  load  disturbance 

r       is  indicated,   and  a  command  input  is  also  shown  as 
^fift  +  YV^)  c.     The  input  block  PFa  is  defined  as 


5 


Y^+Y^j  )c  s-KlY 


1  A 

=    PFa 


It  can  be  noted  for  a  simple  missile,   various  cross  coupling 
terms  are  zero  at  7,ero  angle  of  attack   in     the  asymmetric  system 


K-  ys  *r       >       L? 

=      £-       =    L_       =    1 =o  (7-8) 


1  7 

and  since      CL     -  Q,   then  sin    o( 


°-..=  0 


Examining  Fig.    7-3,   and  noting  the  zero  blocks  resulting  from 
Equation  7-8,   the  standard  determinant  can  be  formed  as  shown 
in  Equation  7-9.    (Table  7-3) 
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Yaw  Subsystem 
-K3s+K2 
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TABLE  7-3 


Equation  7-9 


130 


Therefore,   the  selection  of  input  nodes  in  the  order  of 
~*  )   r     /&  &   and   '       has  allowed  the  formation  of  a  yaw 
and  roll  subsystem.      The  array  clearly  shows  the  two  sub- 
systems and  the  two  attendant  cross  coupling  areas 

Specifications  for  dynamic  response  of  simple  missiles  is 
often  most  rigid  at  the  zero  angle  of  attack  condition.     It  can  be 
noted  from  Equation  7-9  that  compensation  can  be  made  separately 
in  the  tvo  subsystems.     After  these  specifications  for  zero  angle 
of  attack  are  met,   the  cross  coupling  consideration  may  be 
brought  into  the  problem  to  check  design  requirement  at  other 
angles  of  attack. 

It  is  evident  that  this  resultant  block  diagram  is  in  a  form 
which  conveys  a  maximum  amount  of  information  to  the  designer. 
With  proper  manipulations  the  design  of  a  well  controlled  aircraft 
or  missile  becomes  a  much  simpler    problem  . 

No  attempt  is  made  here  to  apply  the  results  just  obtained. 
The  purpose  of  this  section  was  to  indicate  how  the  block  diagram, 
determinant  technique  can  be  applied  to  the  simplification  of 
complex  airframe  control  problems.     It  is  felt  that  the  preceding 
considerations  have  accomplished  this. 
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8.     THE  GENERATION  OF  COMPLEX  ZEROS 
8.1     Introduction 

Any  closed  loop  system  which  can  be  represented  by  a  block  dia- 
gram has  a  characteristic  equation  which  may  always  be  manipulated 
into  the  form 

n 
K    JT   (s  +  Z„) 

G        =  -  1  (8-1) 


n  c 

m 


s       TT    (s  +  p.) 

J=l  J 

Where,   04K  <Oo,    -  oo  <  m  <  <£> ,   and  the  z    and  p    may  be  real  and/ or 

1  J 

complex.       K  is  assumed  to  represent  a  variable  gain  element  of  the 

system.     G     represents  the  transfer  function  of  a  compensator  which  may 

be  in  any  single  feed  forward  or  feedback  path,    or  may  be  cascaded  in  the 

main  transmission  channel.     G     is  of  the  form 

c 

n 

k      TT   (s  +  z  ) 

C  1 

G        = = (8-2) 

c  n 

m  .  . 

s     ji.i  ('  +  '»J> 

The  definitions  of  the  terms  in  Equation  8-2  are  the  same  as  Equation 

8-  1  .     The  total  gain  of  Equation  8  .  1  is  then 

KDT    =  KK  (8-3) 

RL  c 

The  selection  of  a  transfer  function  for  G    ,   besides  being  based  upon 

c 

physical  realizability  is  guided  by  the  desire  to  force  the  roots  of  the 
characteristic  equation  to  lie  in  such  locations  on  the  s-plane  that  the 
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transient  response  of  the  system  to  be  compensated  will  meet  specifica- 
tions without  affecting  steady  state  accuracy.     Steady  state  accuracy 
specifications  are  normally  met  before  compensation  by  setting  system 
gain  to  obtain  an  acceptable  steady  state  error  coefficient       However, 
it  is  conceivable  that  compensation  could  be  required  to  meet  steady  state 
specifications  . 

It  is  common  practice  to  select  two  complex  conjugate  root  loca- 
tions and  then  to  choose  G     so  that  the  two  selected  locations  will  be  roots 

c 

of  the  characteristic  equation.     Generally,   it  is  hoped  that  the  selected 
roots  will  prove  to  be  dominant,   but  a  certain  amount  of  experience  and 
foresight  is  required  in  order  to  ensure  that  this  end  will  be  attained. 
Even  with  experience,    some  trial  and  error  is  frequently  involved.     In 
any  case,    only  calculation  of  the  transient  response  will  determine  whether 
or  not  a  particular  solution  is  acceptable.     The  concept  of  dominance 
states  that  if  the  residues  at  certain  roots  of  the  characteristic  equation 
are  large  compared  to  the  residues  at  all  other  roots,   and  if  the  real 
parts  of  those  certain  roots  are  several  times  smaller  than  the  real  parts 
of  any  other  roots,    then  those  certain  roots  are  dominant       A  transient 
response  calculated  using  only  the  dominant  roots  will  be  sufficiently 
accurate  for  engineering  purposes       All  roots  must  be  included  however 
when  calculating  the  phase  angle  of  the  transient  response.     Note  that 
this  concept  may  significantly  reduce  the  labor  involved  in  calculating 
the  transient  response  of  higher  order  systems 
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The  common  use  of  passive  networks  to  effect  compensation  is 
with    G     containing  only  negative  real  poles  and  zeros,   and  usually  with 
K    ^L  1  .     It  is  possible  to  generate  complex  poles  or  zeros  with  passive 
networks,   but  it  is  difficult  to  use  these  networks  practically  due  to  large 
attenuation  in  the  compensator  and  extreme  sensitivity  of  the  complex 
pole  or  zero  location  to  variations  of  components  within  the  network.   In 
some  cases  an  additional  amplifier  is  required  to  correct  for  the  attenua- 
tion of  a  passive  compensator,    so  the  compensator  could  have  been  de- 
signed as  an  active  compensator  from  the  start, 

The  terms  "active  network"  and  "active  compensator"  are  inter- 
changeable and  refer  to  networks  which  contain  passive  components,  a 
summer,   and  one  or  two  amplifiers. 

Frequently,    complex  poles  and/or  complex  zeros  are  present  in  a 
system  before  compensation  is  attempted.     Having  networks  available  to 
use  as  compensators  which  will  produce  complex  poles  or  zeros  will  in 
some  instance  make  the  compensation  problem  easier.     The  active  net- 
works investigated  in  this  report  are  flexible  in  that  the  complex  poles 
or  zeros  resulting  from  them  may  be  relocated  on  the  s-plane  by  varying 
either  an  amplifier  gain  or  the  magnitude  of  one  or  more  of  the  passive 
elements  in  the  network. 

The  passive  networks,   which  are  essential  to  active  compensators, 
may  be  either  RC  or  RL  networks.     However,    RL  networks  pose  some 
problems  due  to  difficulties  in  determining  and  allowing  for  the  resis- 
tance of  the  inductive  elements. 
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Active  compensation  is  not  envisioned  as  replacing  passive  com- 
pensation,  but  rather  as  supplementing  it  by  providing  solutions  to 
compensation  problems  not  easily  attainable  with  passive  compensators 
alone.     In  some  instances  it  may  prove  desirable  to  combine  active  and 
passive  compensators  to  achieve  desired  results. 

Two  major  disadvantages  of  active    compensators  are  cost  and  com- 
plexity.    Both  of  these  factors  arise  from  the  fact  that  amplifiers  and 
summers  are  required  as  parts  of  the  compensator, 

The  passive  networks,   their  transfer  functions,   and  s-plane  pole- 
zero  arrays  which  will  be  used  as  basic  components  in  the  active  networks 
are  depicted  in  Fig.   8.1.     They  are  held  to    simple  configurations  pri- 
marily to  provide  attainable  goals  for  this  investigation.     In  addition, 
all  compensation  problems  attempted  could  be  solved  using  compensators 
which  required  only  these  simple  networks.     The  transfer  function  of 
Fig.    8.  1  are  defined  on  Table  8-1  for  algebraic  manipulation. 
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Fig.   8,1     Networks,   transfer  functions,   and  s-planes  of  passive  networks 
used.     The  network  of  Fig.    8.1  (c)  cannot  ordinarily  be  used  in  the  main 
transmission  path  unless  there  exists  a  parallel  forward  path  through  which 
dc  components  may  pass. 
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TABLE  8-1 


s  +  R    C  £  s  +   Z 
i =  d 


s       + 


(s+-^-   ) 


R2  ^ 

S+Ri+Rz         L 
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Rl  +  R2  1 

■  +  -4—^-  -^  s  +  pd 

R2  R    C  d 


RZ       ^  1         ^ 

(R1  +  R2)  C  8 


1  +  TC  s  +  Zd 


R2 

s  +      —  A  s  +  zd 


R.  +  R9 
,12  s    +  p  , 

s    + ^d 


R  R  A  OC(.+  Z) 
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R    +R_                  L 
1         2 (s+p) 


g 


s    +  R/L  s      +  p^, 

d 

Definitions  of  the  transfer  functions  of  Fig.   8-1  for  algebraic  manipulation 
The  subscript  d  refers  to  poles  or  zeros  resulting  from  a  lead  network. 
The  subscript  g  refers  to  poles  or  zeros  resulting  from  a  lag  network. 
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8.2    Investigation  of  Some  Active  Networks 

Not  all  active  networks  will  yield  complex  poles  or  zeros.     It  is 
therefore  necessary  to  investigate  several  combinations  of  amplifiers, 
summers  ,   and  passive  networks  to  determine  some  of  the  configurations 
which  will  produce  the  desired  complex  poles  and  zeros. 

Two  basic  network  configurations,    one  a  parallel  feed  forward  and 
the  other  a  feedback,   were  selected  for  investigation.     Solutions  were 
obtained  by  placing  a  lag  or  a  lead  transfer  function  and  amplifier  in 
each  path  available.     The  overall  transfer  function  of  the  network  was 
then  determined  and  investigated  for  the  existence  of  complex  poles  or 
zeros . 

Fifteen  different  networks  were  investigated  using  this  plan.     Some 
were  unsuccessful  in  that  they  did  not  yield  the  desired  complex  poles 
and  zeros.     Only  those  producing  the  desired  result  are  treated  here. 
There  are,    undoubtedly  other  simple  active  networks  which  will  produce 
complex  poles  and  zeros. 

Lead-lead  feed  forward  difference  compensator 

Fig.   8.2  is  the  block  diagram  of  this  compensator.     The  transfer 
function  of  Fig.    8.2  is 


V 
o 


A    ■(.+  » d)(   s+pd2)-  A2(s+   ZdZ)(s+    pdl) 


3  -.  (s  +  p      )  (s  +  p       ) 

dl  d2  (8-4) 


For  factoring  by  root  locus,   the  numerator  of  Equation  8t4  may  be 
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manipulated  into  the  form 

,)  (s  +  PH2) 

(8-5) 


VS  +  *di)(s  +  Pd2> 


A2    (s  +  ^d2)(s  +  pdl) 


=     1 


Note  that  Equation  8-5  is  a  zero  degree  root  locus.     The  possible  root 
loci  arising  from  Equation  8-5  are  Figs  .   8.3,   8.4,   8.5,   and  8.6.     It 

can  be  seen  from  Figs.   8.5  and  8.6  that  complex    roots    of  Equation  8-5 

Al 
will  obtain  for  certain  values  of     ™ —    .     These  complex  roots  will  become 

A2 

complex  zeros  in  Equation  8-4  and  result  in  Equation  8-4  having  the  form 

Vo  (ArA2)    (s+<r+jU))(s+<r-jW)  _ 


V.  (s  +  p  .   )     (   s  +  p    „  ) 

i  *di '    *  *d2 


This  transfer  function  is  of  further  interest  because  p..    and  p  n_  ,   the 

dl  d<£ 

real  poles,   may  be  placed  far  out  in  the  left  half  plane  if  desired. 

There  can  be  no  complex    root    resulting  from  the  root  loci  of  Figs. 
8.5  or  8.6.     Therefore  they  are  of  no  further  interest  in  this  investigation, 
Lead-lead  negative  feedback  compensator. 

Fig     8.7  is  the  block  diagram  of  this  compensator.     The  transfer 

function  of  Fig.   8.7  is 

(S+    *dl)(S+Pd2> 

V 

__?_      -       _ —  (8-7) 

For  factoring  by  root  locus,   the  denominator  of  Equation  8-7  may  be 
manipulated  into  the  form 
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Fig.   8.2    Block  diagram  of  the  lead-lead  feed  forward  difference 

compensator. 
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Fig.   8.4      A  possible  root  locus  of  Equation  8-5 
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Fig.   8.5      A  possible  root  locus  of  Equation  8-5. 
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Fig.   8,6      A  possible  root  locus  of  Equation  8-5 
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Fig.   8.7      Block  diagram  of  the  lead-lead  negative  feedback  compensator 


Fig.   8.8      A  possible  root  locus  of  Equation  8-8 
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V  +  'V(B  +  V 


(s  +  Pdl)(s  +  pd2) 


(8-8) 


Fig.    8    8  is  the  root  locus  of  Equation  8-8.     Any  locations  of  p.,,   p,_, 

dl         ad 


'dl 


'dl 


or  z        will  produce  the  locus  of  Fig.   8.8  unless     z 
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d2 


>    P 


'dl 


or 
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'62 


The  complex  roots  which  obtain  for  certain  values  of  A 


1 


will  become  complex  poles  in  Equation  8-7  and  result  in  that  equation 
having  the  form 


o 


(s+Zd1)(S  +  Pd2) 


(8-9) 


i  (Al  +  1)  (  s  +  o-+  jW)  (  s  +cr-j">) 

This  transfer  function  is  of  further  interest  because  one  of  the  real 
zeros  of  Equation  8-9  may  be  placed  close  to  the  origin  while  the  other 
may  be  placed  as  far  out  in  the  left  half  plane  as  desired.     Note  from 
Fig.    8.8  that  the  complex  poles  will  always  be  in  the  left  half  s -plane. 
Lag-lag  feed  forward  difference  compensator. 

Fig.    8.9  is  the  block  diagram  of  this  compensator.     The  transfer 
function  of  Fig.   8.9  is 


A     cCj  (  s  +   Z      )  (s  +  p      )-  A20(2(s  +   Z      )(b  +  P     ] 


V 


(  s  +  p  )  (  s  +  Pg2) 


(8-10) 


For  factoring  by  root  locus,   the  numerator  of  Equation  8-10  may  be 

manipulated  into  the  form 
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V; 


Fig.   8.9      Block  diagram  of  the  lag- lag  feed  forward  difference 
compensator . 
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Fig.   8.10      A  possible  root  locus  of  Equation  8-11  . 


oO  A,  j 


ncreasing    J 


P 


uO 


r 


a£Hfc  increasing 


<r 


Fig  .8.11      A  possible  root  locus  of  Equation  8-  1 1 


144 


IK) 


•jt — e— ■ 


% 


*- 


-<r 


Fig.   8.12      A  possible  root  locus  of   Equation  8- 11 
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Fig.   8.13      A  possible  root  locus  of  Equation  8- 1 1 
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Ai^('  +  V(*  +  V     ,  (811, 

Az<^(.  +  *g2)(-  +  Pgl)      "  l    " 

Note  that  Equation  8-11  is  a  zero  degree  root  locus.     The  possible  root 
loci  arising  from  Equation  8-11  are  Figs  .   8.10,   8.11,   8.12,   and  8.13. 

It  can  be  seen  from  Figs.    8.10  and  8.11  that  complex  roots  of  Equation 

Al  o(i 

8-11  will  obtain  for  certain  values  of     - — — —    .     These  complex  roots  will 

A2^2 
become  complex  zeros  in  Equation  8-10  and  will  result  in  the  equation 

having  the  form 

V  (A    dC,    -  A      C<J  (s  +o-+ju>)  (s  +cr-  j  to) 

-2-      =      -J— 1 1 1 (8-12) 

Vi  (s+Pgl)  (s+Pg2) 

This  transfer  function  is  of  further  interest  because  p    ,  ,   and  p    _,   the 

gl  g2 

real  poles  of  Equation  8-12,  may  be  placed  fairly  close  to  the  origin  in 
the  left  hand  s -plane. 

There  can  be  no  complex  roots  resulting  from  the  root  loci  of  Figs. 
8.  12  and  8.13.     Therefore  they  are  of  no  further  interest  in  this  investi- 
gation. 
Lag- lag  negative  feedback  compensator. 

Fig.   8.  14  is  the  block  diagram  of  this  compensator.     The  transfer 
function  of  Fig .    8.  14  is 

V_o_    _      <*,    <s+egl"s+Pe2) (8.13) 

Vi  (s  +  pgl)  (  s  +  pg2)  +  AjCijOCj,  (s  t  Zgl)  (s+Zg2) 

For  factoring  by  root  locus,    the  denominator  of  Equation  8-13  may  be 

manipulated  into  the  form 
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Fig.   8.  14      Block  diagram  of  the  lag-lag  negative  feedback  compensator 
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Fig.   8.15      A  possible  root  locus  of  Equation  8- 14, 
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a  cCTOC2(s  +  s  ^(8+*    ) 
— ai i -  _  i  (8-14) 

The  root  locus  of  Equation  8-14  is  Figr    8.  15.     It  can  be  seen  from  Fig. 

8.15  that  complex  roots  of  Equation  8-14  will  obtain  for  certain  values 

of  A  dC  cC   •     These  complex  roots  will  become  complex  poles  in  Equation 
i     1     ^ 

8-13  and  will  result  in  that  equation  having  the  form 

Vo  <*1  (  S  +  V  (  S  +  Pg2> 

_2_     =      i I£ if (8-15) 

i  Ai°Ci^2+   1)  (  s  +CT+  jU))  (   s  +CT-  JOJ) 

This  transfer  function  is  of  further  interest  because  one  of  the  real  zeros 

of  Equation  8-15  may  be  placed  near  the  origin  while  the  other  may  be 

as  far  out  in  the  left  hand  s-plane  as  desired.     The  complex  poles  will 

always  be  in  the  left  hand  s-plane. 

Note  that  Fig.   8.15  will  be  the  root   locus  of  Equation  8-14  unless 

|p_|->|zl       orlp      I  >       I  z    „  I      .     In  neither  of  these  cases  could  com- 
I    g2K  I    gl|  I    gl|  I     g2' 

plex  roots  be  obtained  from  Equation  8-14.     They  are  therefore  of  no 

interest  in  this  investigation. 

Lag-lead  feed  forward  summing  compensator. 

The  block  diagram  of  this  compensator  is  Fig.   8.  16.     The  transfer 

function  of  Fig.   8.1 6  is 

V  Al(<(s+Z     )  (   s  +  p)  +  A.  (   s  +  a   )(S+  p   ) 

-^    =      -J g d 2 d^ g_       (g_l6) 

1  (s  +  p   )  (   s  +  p   ) 

For  factoring  by  root  locus,   the  numerator  of  Equation  8-1 6  may  be 

manipulated  into  the  form 
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Fig.   8.16      Block  diagram  of  lag-lead  feed  forward  summing  compensator 
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Fig.   8.17      A  possible  root  locus  of  Equation  8- 17 . 
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A1^(s+Zg)(s  +  pd)  ^  (g_17) 

A     (   s  +  Z   )  (   s  +  p     ) 
L  d  g 

Fig .   8  .  17  is  the  root  locus  of  Equation  8-  17  .     It  can  be  seen  by  inspec- 
tion of  Fig.   8.17  that  complex  roots  will  obtain  for  certain  values  of 

— so  long  as    I   z     !>*     z  ,        and      p  .    >  |p  .     These  complex 

a2  e      i  g    i  d  pd  •*  rg  i  F 

roots  will  become  complex  zeros  in  Equation  8-1 6  resulting  in  that  equa- 
tion having  the  form 

V  A1oC+  A.)  (  s  +cr+  jW)(   s  +  cr  -  j  t»» 

-1 1 (8-18) 


V 

i  (s  +  p  )  (s  +  p  ) 

g  d 

This  transfer  function  is  of  further  interest  because  one  of  the  real 
poles  of  Equation  8-18  may  be  placed  close  to  the  origin  in  the  left  hand 
s-plane  while  the  other  may  be  placed  as  far  out  on  the  negative  real  axis 
as  desired.     The  complex  zeros  will  always  be  in  the  left  hand  s-plane. 
Lag-lead  positive  feedback  compensator;  lag  in  the  feedback  path. 

The  block  diagram  of  this  compensator  is  Fig.   8.  18.     The  transfer 
function  of  Fig.    8.18  is 

V  (   s  +  Z  ,)(   s  +  p   ) 

-    i i (8-19) 


i         ( s  +  p   )  (   s  +  p   )  -  A  dC(   s  -   Z   )( s  +  Z   ) 
g  d  1  g  d 

For  factoring  by  root  locus,   the  denominator  of  Equation  8-19  can  be 

manipulated  into  the  form 

A    oC    (s  +   Z   )  (   s  +  Z  ) 
— i 1 ?■ =   1  (8-20) 

(   s  +  p     )  (   s  +  p     ) 
g  d 
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Note  that  Equation  8-20  is  a  zero  degree  root  locus.     Possible  root  loci 
resulting  from  Equation  8-20  are  Figs.   8.19     8.20>   8.21,   and  8.22.     It 
can  be  seen  by  inspection  of  Figs.   8.21  and  8.22  that  complex  roots  of 
Equation  8-20  will  obtain  for  certain  values  of     A    oC\     These  complex 
roots  will  become  complex  poles  in  Equation  8-19  and  result  in  that  equa- 
tion having  the  form 

V  (  s  +  2  ,)  (  s  +■  p  •  ) 

T°-      =      ■ * -1 (8-21) 

i  (  1  -  A^)  (  s  +(T+  j  U))  (  s  +<T  -  jui) 

This  transfer  function  is  of  further  interest  because    z     and  p    ,   the  real 

d  g 

zeros  of  Equation  8-21  may  be  placed  as  close  to  the  origin  of  the  s-plane 
as  desired.     There  can  be  no  complex  roots  resulting  from  the  root  loci 
of  Figs.   8.19  and  8.20.     Therefore,   they  are  of  no  further  interest  in 
this  investigation. 
Lag- lead  positive  feedback  compensator;  lead  in  the  feedback  path. 

Fig.   8,23  is  the  block  diagram  of  this  compensator.     The  transfer 
function  of  Fig.   8,23  is 

V  (   s  +  2     )  (   s  +  Pj  ) 

-2.        = i d (8_22) 

V-  (  s  +  p  )  (  s  +  p  )  -  A.<*(s  +   &,)  (s  +  2   ) 

g  d  1  d  g 

For  factoring  by  root  locus,   the  denominator  of  Equation  8-22  may  be 

manipulated  into  the  form 

A,C((   s  +  2 J  (   s  +   2     ) 
-J * i =  1  (8-23) 


(  s  +  pd)  (   s  +  p   ) 
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Fig.   8.  18      Block  diagram  of  the  lag-lead  positive  feedback  compensator 
Lag  in  the  feedback  path. 
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Fig.   8.19      A  possible  root  locus  of  Equation  8-20 
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Fig.   8.20      A  possible  root  locus  of  Equation  8-20 


152 


^i_ 


3      t 


/  c<  increasing 


Fig.   8.21      A  possible  root  locus  of  Equation  8-20 
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Fig.   8.22      A  possible  root  locus  of  Equation  8-20 


153 


',-*6 


s  +   Z. 


TT 


Pd 


C((s+Zg) 


s  +  p. 


Fig.   8.23      Block  diagram  of  the  lag-lead  positive  feedback  compensator 
Lead  in  the  feedback  path. 
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Equation  8-23   is  identical  to  Equation  8-20,     Therefore  the  root  loci 
of  Equation  8-23  will  be  those  of  Figs.   8.  19,   8.20,   8.21,   and  8.22. 
Thus;    Equation  8-22  will  be  of  the  form 

V  (s  +  ^  )  (  s  +  p    ) 

^        =        ^— ^ (8-24) 

1  (  l  -  Ap£)  (  s  +cr+  jw)  (  s  +cr-  ju>) 

The  same  comments  may  be  made  for  Equation  8-24  as  for  Equa- 
tion 8-21  except  that  z     and  pia   the  real  zeros  of  Equation  8-24  be  placed 

g  d 

far  out  in  the  left  hand  s-plane  if  desired. 
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8.3    Algebraic  Equations  for  the  Loci 

The  complex  zeros  generated  by  the  active  compensators  all 
lie  on  circles.     For  some  purposes  the  equations  of   these  circles 
are  of  interest.     The  derivations  are  simple  and  are  not  given 
here.     The  results  fall  naturally  into  two  sets  of  algebraically 
similar  equations,  which  are  listed  in  Tables  8-2  and  8-3. 
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9.     MANIPULATIONS  AND  DESIGN 

9.  1  Introduction 

The  equations  of  the  root  locus  circles,  while  of  definite  interest,   are 

not  suitable  to  work  with  when  designing  active  compensators.     A  convenient 

form  may  be  obtained  by  manipulation  of  the  equations  of  Table  8-2  and  8-3. 

For  complex  zeros  the  relationship  of  Table  8-2  may  be  manipulated  into 

the  form: 

(9-1) 

(Zl"Z2+P2  "Pi*  {x2+y2)  +2x  <Zi  P2_Z2Pl)  +  P1P2  (Zl"Z2)  +  Z1Z2(P2"P1)  =  ° 
for  use  in  design.  Table  9.  1  gives  the  particular  form  of  Equation  9-1 
applicable  to  each  of  the  complex  zero  producing  compensators.  It  will  be 
shown  that  four  of  the  six  variables  in  Equation  9-1  may  be  chosen  to  meet 
the  needs  of  a  specific  problem.  Equation  9-1  will  then  reduce  to  an  equa- 
tion in  two  unknowns  which  may  be  chosen  so  as  to  satisfy  that  equation  and 
the  requirement  of  physical  realizability . 

A  significant  simplification  of  Equation  9-1  occurs  if: 
Z    /  p     =   Z    /  p 

1.1*6  (92) 

Z1P2=  Z2P1  =  K1 
Substituting  Equation  9-2  into  Equation  9-1  and  simplifying  yields 

x2+y2  =  K1  (9-3) 

Note  that  due  to  restrictions  imposed  by  Equation  9-2,   Equation  9-3  can  be 
used  only  in  those  active  compensators  which  contain  two  lag  or  two  lead 
networks. 
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For  complex  poles,   the  equations  of  Table  8-2  may  be  manipulated  into 

the  form 

(■Bj  +  a2  -Pj  -  p2)     (7C2+y2)    +  2X(Pj  p2  -  2j  Z^)  + 

flj  22  (  Pj  +  p2)  -  Pl  p2  (  Bj  +  a2) .    0 

for  use  in  design.     Table  9-2  gives  the  particular  form  of  Equation  9-4 
applicable  to  each  of  the  complex  pole  producing  compensators.     It  will 
be  shown  that  four  of  the  six  variables  in  Equation  9-4  may  be  chosen  to 
meet  the  needs  of  a  specific  problem.     Equation  9-4  will  then  reduce  to  an 
equation  in  two  unknowns  which  may  be  chosen  so  as  to  satisfy  that  equa- 
tion and  the  requirement  of  physical  realizability . 

A  significant  simplification  of  Equation  9-4  occurs  if 

Zl   a2    =  Pl  P2    =      K2  (9"5) 

Substituting  Equation  9-5  into  Equation  9-4  and  simplifying  yields 

x2  +  y2    =  k2  (9-6) 

Note  that  due  to  restrictions  imposed  by  Equation  9-5,  Equation  9-6  can 
be  used  only  in  those  active  compensators  which  contain  a  lag  and  a  lead 
network. 
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9.  2    Construction  of  the  Locus    of  Complex  Poles  or  Complex 
Zeros  Producing  a  Constant  Phase  Angle  at  a  Point  in  the 
S-plane. 

It  has  been  shown  by  Carpenter  (19)  that  if  the  angle  at  a 
point  in  the  s-plane  due  to  an  array  of  poles  and  zeros  be  known, 
it  is  possible  to  construct  the  locus  of  complex  conjugate  poles  or 
complex  conjugate  zeros  which  will  cause  this  point  to  be  on  a 
root  locus o     This  locus  of  complex  poles  or  zeros  proves  to  be  an 
arc  of  a  circle.     Carpenter's  derivation  is  as  follows: 

Figure  9.  1   shows  the  necessary  construction  in  proving  the 
locus  of  complex  conjugate  poles  or  zeros,   R  and  R',   producing  a 
constant  phase  at  a  point  ps   in  the  s-plane,   is  an  arc  of  a  circle 
passing  through  that  point.     The  center  of  the  circle  lies  on  the 
real  axis . 

From  the  figure: 

a  =  90°+C*+/^ 
b  =  90°  +QC 
a  +  b  =^+  20C+   180  =M=angle  contributed  by         (9-6) 

complex  pair  at  point  p. 

But/^+  2  CC+  y>=   180°  (9_7) 

Subtracting  9.7  from  9.6 

360     -  M  =  ^  =  a  constant 

The  vertex  of   y7,  therefore,   lies  on  an  arc  pp' ,   of  a  circle 
whose  center  is  on  the  real  axis,   by  a  fundamental  theorem  of 

plane  geometry. 
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Since  the  arc  always  intersects  the  real  axis,   the  arc 
may  be  quickly  constructed  by  first  locating  that  intersection. 
The  intersection  is  found  by  constructing  the  angle, 

0    =  M/2  (9-8) 

at  the  point  P. 

It  can  also  be  shewn  that  the  construction  angle  at  the  point 
P  which  locates  the  center  of  the  circular  arc  on  the  real  axis  is 

r=  m  (9-9) 

The  necessary  construction  to  show  this  is  Fig.  9-2.     F  is  the 
intersection  of  the  circular  arc  with  the  real  axis  defined  by  the 
angle  0.     HG  is  the  perpendicular  bisector  of  PF .     Since  the  center 
of  the  circular  arc  is  on  the  real  axis,   the  intersection  G  must  be 
the  center  of  the  circular  arc,   by  plane  geometry.     From  Fig.  9-2; 
0+g=9O°  (9-10) 

20+  2g  =  180°  (9-11) 

2g  =  h  (9-12) 

Subtracting  Equation    r;.5    from  0.4  obtains 

20    =       180°  -  h  (9-13) 

h+f=  180°  (9-14) 

f=  180°  -  h  (9-15) 

Subtracting  Equation  9-15  from  9-13 

20  =  r  =  M 
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f  rom  Equation  $.$,   thus  proving  the  statement. 

With  the  center  of  the  circular  arc  and  its  point  of  inter- 
section with  the  real  axis  known  from  0  and  K,   the  arc  may- 
be easily  constructed. 

This  construction  provides  a  simple  means  of  finding  where 
on  the  s-plane  complex  poles  or  zeros  may  be  placed  in  order  to 
force  a  selected  point  to  be  on  a  root  locus.     This  technique  will 
not  always  be  needed.     For  example  when  it  is  desired  to  cancel 
existing  complex  poles  with  complex  zeros,   the  desired  location 
of  the  complex  zeros  is  specified  by  the  location  of  the  complex 
poles.      However,   in  the  more  general  cases,   this  construction 
becomes  a  powerful  tool  in  problem  solution. 

Note  that  the  circular  arcs  resulting  from  this  construction 
will  tell  the  designer  where  complex  compensating  poles  or  zeros 
should  be  placed  in  order  to  cause  a  desired  point  to  be  on  a  root 
locus.     The  circular  equations  for  the  active  compensators  of 
Section  9-  1  then  allows  one  to  place  complex  zeros  at  the  required 
locations . 
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9.3    A  Technique  for  Compensation  Using  Active  Networks 

Relocation  of  the  root  loci  of  the  characteristic  equation  so  that 
they  pass  through  a  specified  point  is  not  in  itself  any  great  problem. 
Simultaneously  meeting  the  requirements  of  dominance,    steady  state 
accuracy,   and  transient  response  is.    however,    quite  another  matter.     A 
considerable  number  of  techniques  exist  for  the  simultaneous  solution  of 
the  steady  state  accuracy  and  root  location  problems.     See  for  example 
20,   21    .     These  techniques  apply  only  to  compensation  with  negative  real 
poles  and  zeros  and  do  not  easily  extend  to  the  use  of  complex  poles  and 
zeros.     The  technique  presented  here  requires  some  intuition  and  trial 
and  error  as  thus  far  no  more  straight  forward  approach  has  been  found 
For  the  trial  and  error  phases,   the  Electro-measurements,   Inc.  ,   ESIAC 
algebraic  Computer,   Model  10  was  a  valuable  tool. 

The  techniques  of  solution,   whether  for  complex  poles  or  complex 
zeros     are  practically  identical.     Therefore  they  will  be  presented 
separately  to  avoid  any  possibility  of  confusion. 

(A).     Solutions  requiring  complex  zeros. 

a.     Select  the  compensation  path.     Then  determine  the  open  loop 
transfer  function  and  manipulate  the  characteristic  equation  into  the 
form: 


n 

K:T1        (S+Z.)  „  ,  (9-17) 


m  TT      /  > 

s      l  =  i      (  s  +  V 


G    =  -  1 

c 


168 


b.     Plot  the  z     ,   p.,   and  s       on  the  s -plane, 
i         J 

c       Select  the  real  poles  of  G     so  as  to  be  most  convenient. 

c 

Convenience  will  be  determined  by  the  particular  problem.     Plot  these 
poles  on  the  s-plane.     This  will  usually  determine  which  of  the  complex 
zero  producing  active  networks  is  to  be  used. 

d.  Select  the  desired  locations,   P  and  P' ,    of  two  complex  con- 
jugate closed  loop  roots.     Plot  P  and  P'  on  the  s-plane.     As  mentioned 
previously,   it  is  hoped  that  these  roots  will  prove  dominant,   but  there  is 
no  way  of  ensuring  this. 

e.  Calculate  the  phase  angle  at  P  or  P1  due  to  z.    ,   p. ,    s      ,   and 
the  real  poles  of  G    .     Then  construct  the  locus  of  complex  zeros  in 

the  manner  described  in  Section  9.2,   which  will  result  in  P  and  P'  being 
on  a  180    degree  locus. 

f.  Choose  two  complex  conjugate  points,   Q  and  Q' ,    on  this  con- 
structed locus.     These  will  be  the  complex  zeros  of  G    . 

g.  Determine  the  root  locus  gain,   K  which  will  cause  roots 

to  be  at  P  and  P'  with  the  z ,    ,   p.  ,    s      ,   and  the  chosen  real  poles  and 

i  J 

complex  zeros  of  G    . 

c 

h.     Substitute  K         and  the  real  poles  and  complex  zeros  of  G 

RL  c 

into  the  open  loop  transfer  function  to  determine  if  steady  state  accuracy 

specifications  are  met.     If  not,    choose  two  new  locations  of  the  complex 

zeros  of  G     and  repeat  steps  f  and  g, 
c 

i.     Determine  x  and  y  from  the  chosen  locations  of  the  complex 


zeros  of  G 
c 
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j  .     Substitute  x  and  y  and  the  selected  real  poles  of  G     into  the 
appropriate  equation  of  Table  9-1  to  determine  the  algebraic  relation- 
ship between  the  real  zeros  of  the  compensator  passive  networks. 

k.     Construct  a  graph  of  this  equation  or  evaluate  it  for  several 
values  of  one  variable  and  select  any  convenient,   physically  realizable 
values  for  the  zeros  of  the  passive  networks.     The  z/p  ratios  of  lead  net- 
works should  be    J    =0.1. 

1.     Substitute  the  real  poles  and  zeros  of  the  passive  networks  into 
the  numerator  root  locus  equation  of  Table  9-1   corresponding  to  the 
active  network  chosen.     Plot  this  equation  and  the  selected  complex  com- 
pensator zeros  on  another  s-plane. 

m.     Determine  the  value  of  A     /A      which  will  cause  the  roots  of 

1         " 

the  numerator  of  G     to  lie  at  the  selected  locations  for  the  complex 

c 

zeros . 

n,     Using  the  results  of  steps  g  and  m     A     and  A     may  be  calcu- 
lated  and  the  compensator  design  is  completed. 

o.     The  transient  response  must  be  computed  and  compliance  with 
specifications  determined. 

(B).     Solutions  requiring  complex  poles. 

a.     Select  the  compensation  path.     Then  determine  the  open  loop 
transfer  function  and  manipulate  the  characteristic  equation  into  the 

form: 

n 

K   JT      (s  +  2.  .) 

G     =   -   1  (9-18) 


c 


m     JJ- 


s  }\     (    s  +  p.) 
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b.  Plot  the  z.  ,   p     ,   and  s       on  the  s -plane. 

1        J 

c.  Select  the  real  zeros  of  G     so  as  to  be  most  convenient. 

c 

Convenience  will  be  determined  by  the  particular  problem.     Plot  these 
zeros  on  the  s-plane.     This  selection  will  determine  the  particular  active 
network  to  be  used. 

d.  Select  the  desired  locations,   P  and  P1,    of  two  complex  conju- 
gate closed  loop  roots.     Plot  P  and  P'  on  the  s-plane       As  mentioned 
previously,   it  is  hoped  that  these  roots  will  be  dominant;   but  there  is  no 
way  of  ensuring  this. 

c.     Calculate  the  phase  angle  at  P  or  P'  due  to  z.  ,  p.,    s      ,   and 

the  real  zeros  of  G    .     Then  construct  the  locus  of  complex  poles  in  the 

c 

manner  described  in  Section  9.2,   which  will  result  in  P  and  P'  being  on 
a  180  degree  locus. 

f.  Choose  two  complex  conjugate  points  on  this  constructed  locus. 
These  will  be  the  complex  poles  of  G    . 

g.  Determine  the  root  locus  gain,   K        ,   which  will  cause  roots 

R  J-/ 

to  be  at  P  and  P'  with  the  z„ ,   p. ,    s      ,   and  the  chosen  complex  poles  and 

i        J 

real  zeros  of  G    . 
c 

h.     Substitute  K         and  the  complex  poles  and  real  zeros  of  G 

RL  c 

into  the  open  loop  transfer  function  of  the  system  to  determine  if  steady 

state  accuracy  specifications  are  met.     If  not,    choose  new  locations  for 

the  complex  poles  of  G     and  repeat  steps  f  and  g. 

c 

l.     Determine  x  and  y  from  the  chosen  locations  of  the  complex 

poles  of  G 
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j  .     Substitute  x  and  y  and  the  selected  real  zeros  of  G     into  the 
appropriate  equation  of  Table  9-2  to  find  the  algebraic  relationship  be- 
tween the  real  poles  of  the  passive  networks. 

k.     Construct  a  graph  of  this  equation  or  evaluate  it  for  several 
values  of  one  variable  and  select  any  convenient     physically  realizable 
values  of  the  passive  network  poles.     The  z/p  ratio  for  lead  networks 
should  be     o    =0.1. 

1.    Substitute  the  values  of  the  poles  and  zeros  of  the  passive  net- 
works into  the  denominator  root  locus  equation  of  Table  9-2  correspond- 
ing to  the  active  network  chosen.     Plot  this  equation  and  the  selected 
complex  compensating  poles  on  another  s-plane. 

m.     Determine  the  value  of  A     which  will  cause  the  roots  of  the 
denominator  of  G     to  lie  at  the  selected  locations  of  the  complex  compen- 
sating poles.     This  completes  the  design  of  the  compensator. 

n.     The  transient  response  must  now  be  computed  and  compliance 
with  specifications  determined. 


172 


9  .  4    Numerical  Examples 

This  section  applies  the  techniques  of  solution  to  further  clarify 
their  use.     The  examples  presented  herein  are  straightforward  so  as  to 
illustrate  the  use  of  active  networks  as  compensators  without  becoming 
unnecessarily  involved  in  other  considerations.     It  is  assumed  that  active 
compensation  was  necessary  in  these  examples  due  to  unspecified  condi- 
tions .     The  subsections  of  each  numerical  example  correspond  to  the   steps 
in  the  applicable  technique  of  solution. 

Example  1 

Given: 


The  feedback  control  system  of  Fig.    9.3    K      =5.0. 

Requirements: 
Use  cascade  compensation.     Desired  root  location  to  be  such  that 
0.5^  j^—  0.7;    settling  time  =   1  .0  seconds.     K     not  to  be  reduced. 

Solution: 


a.     The  compensated  open  loop  transfer  function  of  Fig.    9.3  is 

F     =     5 - =   G  (9-19) 

°  s  (   s       +  5s  +   100)  C 

The  characteristic  equation  of  the  compensated  system  is 

1  +       2-* =  0  (9-20) 


s 


(  s     +  5s  +   100) 


b.     The  z   ,   p  ,   and  sm  of  Equation  9-19  are  plotted  on  Fig.   9    4. 
i       J 

c       Select  p         =  -  25,   p         =  -   30,    so  as  not  to  add  any  more  poles 
dl  d2 

in  the  region  near  the  origin,     p        and  pd2    are  also  plotted   on  Fig.    9-4. 
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d.  Select  the  points     s=  -  5  +     j     7  to  be  P  and  P' .   the  roots  of  the 
characteristic  equation  which  will  meet  the  requirements. 

e.  From  Fig.   9-4,   the  lead  angle  required  at  P  is  307 .  3  degrees . 
The  construction  angles  at  P  are 

0    =   153.65  degrees 
P  =   307. 3  degrees 
The  arc  of  the  circle  on  which  complex  compensating  zeros  can  be  located 
was  constructed  on  Fig.    9-4  using  these  angles. 

f .  On  Fig .   9-4,    choose     s  =  +    j  8.8  to  be  Q     and  Q        the  locations 
of  the  complex  compensator  zeros. 

g.  The  root  locus  gain  at  P  using  all  poles  and  zeros  is 

KRL    =  3'25° 
h.     Substituting  this  value  of  gain  and  the  chosen  compensator  poles 

and  zeros  into  the  open  loop  transfer  function  obtains  K     =   3.51  .     This 
does  not  meet  requirements.     Therefore,    choose  two  new  complex  compen- 
sator zeros  a  little  closer  to  the  desired  root  locations.     Let  these  new 
locations  be     s  =  -2.5  +    j  8.4,   Q     and  Q'-,-     With  these  new  zeros,   the 
root  locus  gain  at  P  is 

KRL    =6,578 

Substituting  this  root  locus  gain  and  the  new  complex  zeros  into  the  open 

loop  transfer  function  yields  K     =6.72.     This  meets  the  requirements. 

i.     From  the  chosen  locations  of  the  complex  zeros, 

x  =   -2.5 

y  =  +  8.4 

176 


j  .     Substituting  these  values  of  x  and  y,   p,,   and  p        into  the  working 

1      *dl  *d2  6 

equation  of  the  lead-lead  feed  forward  difference  compensator  of  Table    9-1 
yields 


ZJ1   =   140.34 
dl 


Z  ^  +     5446 
d2 


135. 34  -   2 


d2 


(9-21) 


The  algebraic  relationship  between  zni   and  z  ,„     . 

dl  d2 

k.     Several  values  of  z  ,_    were  substituted  into  this  equation  and 

d2 

z  __   =   -  2.5,    z._  =   -   3.0  were  selected  as  the  zeros  of  the  lead  networks  . 
dl  d2 

1.     Substituting  the  chosen  z       ,    z       ,   p       ,   and  p         into  the  numerator 

dl        d2        dl  d2 

root  locus  equation  of  the  lead-lead  feed  forward  difference  compensator  of 
Table  9-1  yields 


A     (   s  +  2.5)  (   s  +  30) 
A2  (   s  +  3  )  (   s  4-  25) 


=   1 


(9-22) 


Equation  9-22  and  the  chosen  complex  compensator  zeros  are  plotted  on 
Fig.    9.5. 

m.     From  Fig.   9-5,   A   /A     =  0,8368  will  cause  the  numerato 

J.  L-t 

to  be  at  the  chosen  complex  conjugate  locations. 

n.      The  values  of  A     and  A     may  now  be  computed  from  the  equations 

i  *-> 


r  zeros 


A   /A     =  0.8368 

1  £ 


K<ArV  =  KRLr6578 


(9-23) 


Since  K  is  a  variable  gain  element,   any  values  of  K,   A    ,   and  A       which 
satisfy  Equations  9-23  will  yield  satisfactory  results. 
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o.     The  characteristic  equation  of  the  compensated  system  may  now 
be  written  as 


6578  (s+2.5  +  j8.4)(s+2.5-j8.4) 
s  (s     +   10s  +   100  )  (   s  +  25)  (   s  +  30  ) 


1  (9-24) 


The  roots  of  this  equation  were  obtained  by  calculation  on  the  Control  Data 

Corporation  1604  digital  computer  as 

s=  -    5.11  +  j  7.26 

s=  -    5. 11  -  j  7.26 

s=  -    4.69  +  j  11 .69 

s=  -    4. 69  -  j  11 .69 
8=  -40.4 

The  roots  as  s  =~5  .11  +  j  7.26  and  s  =   -  5.11   -j  7.26  correspond  to  those 

selected  early  in  the  solution.     There  is  some  error  due  to  the  graphical 

techniques  required  in  the  solution,   however  the  agreement  is  close  enough 

for  engineering  work.     The  transient  response  may  now  be  calculated  and 

final  agreement  with  requirements  determined. 

.Example  2 . 

Given- 

The  feedback  control  system  of  Fig.    9.6.     K     =100 

7  °  v 

Requirements 
Use  cascade  compensation.     Desired  root  location  such  that  0.5  £■  j 
Jk     0.7,    settling  time  =1.0  seconds.     K     not  to  be  reduced. 

Solution: 


a.     The  compensated  open  loop  transfer  function  of  Fig.   9-6  is 

F    = G  (9-25) 

°        s(s+l)2 
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The  characteristic  equation  of  the  compensated  system  is 

1  + r-     G     =  0  (9-26) 

s  (  s  +  1)  ° 

b.  The  z   ,   p. ,   and  sm  of  Equation  9-26  are  plotted  on  Fig.   9-7. 

c.  Select  p.,     =  -40.   p  ,_  =  -  50  so  any  root  loci  from  these  poles  will 

dl  d2 

not  affect  the  region  where  the  desired  root  will  be  located,     p       and  p       are 

dl  ad 

also  plotted  on  Fig.   9.7. 

d.  Select  the  points  s  =  -4+    j  5  to  be  P  and  P1  ,   the  roots  of  the 
characteristic  equation  which  will  satisfy  the  requirements. 

e.  From  Fig.   9.7,   the  lead  angle  required  at  P  is  203.8  degrees. 
The  construction  angles  at  P  are 

0    =101.9  degrees 

P  =  201. 8  degrees 
The  arc  of  the  circle  on  which  complex  compensating  zeros  can  be  located 
was  constructed  on  Fig.   9-7  from  these  angles. 

f.  Choose  the  points  s  =  -3.2  +    j  2  to  be  Q  and  Q' ,   the  locations 
of  the  complex  compensator  zeros. 

g.  The  root  locus  gain  at  P  and  P1  using  all  poles  and  zeros  is 

K         =   15,316 
R  i-j 

h.     Substituting  this  value  of  gain  and  the  chosen  poles  and  complex 

zeros  of  the  compensator  into  the  open  loop  transfer  function  obtains 

K     =   108.8. 
v 

This  meets  requirements. 
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i.     From  the  chosen  locations  of  the  complex  compensating  zeros- 

x  =  -   3.2 

y  =  +  2.0 

j  .     Substituting  these  values  of  x  and  y,   p  ln   and  p  ,„  into  the  working 

dl  d2 

equation  of  the  lead-lead  feed  forward  difference  compensator  from  Table 
9-1  yields 


ZJO  =  169.424 
d2 


Z      -    .08405 
dl 

ZJn  +   175.824 
dl 


the  algebraic  relationship  between  z,,   and  z,_. 

dl  d2 

k.     Several  values  of     z..   were  substituted  into  this  equation  and 

dl 

z,,   =   -30,    z,„  =  -  34. 95  were  selected  as  the  zeros  of  the  lead  networks  . 
dl  d2 

1.     Substituting  the  chosen  z  zn>   p,,,   and  p         into  the  numerator 

dl        d2        dl  d2 

root  locus  equation  of  the  lead-lead  feed  forward  difference  compensator 
on  Table  9-1  yields 


A     (s  +  30)  (s  +  50) 

A     (  s  +  34.95  )  (  s  +  40) 


=  1 


(9-27) 


Equation  9-27  and  the  chosen  complex  compensator  zeros  are  plotted  on 

Fig.   9.8. 

m .     From  Fig.   9.8  A   /A     =  0 .  9314  will  cause  the  numerator  zeros 

to  be  at  the  chosen  complex  conjugate  locations. 

n.     The  values  of  A     and  A     may  now  be  computed  from  the  equations 

12 


A   /A     =  0.9314 
K(A1-A2)  =  KRL=15'316 


(9-28) 
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Since  K  is  a  variable  gain  element,   any  values  of  K,   A    .   and  A    ,   which 
satisfy  Equations  9-28  will  yield  satisfactory  results. 

o.     The  characteristic  equation  of  the  compensated  system  may   now 
be  written  as 


15.316  (  s  +  3.2  +  j  2)  (  s+  3.2  -  j  2)  =  _   ^  (9_2<?) 

s  (  s  +  l)2  (  s  +  40  )  (  s  +  50) 


The  roots  of  this  equation  were  obtained  by  calculation  on  the  Control  Data 

Corporation  1604  digital  computer  as 

s=  -    3.39  +  j  5.4 
s=  -    3.39  -  j  5.4 
s=  -    4.33 
s=  -19.73 
s=  -61.06 

The  roots  as  s  =  -3.39  +  j  5.4  and  s=-3.39-j5.4    correspond  to  those 
selected  early  in  the  solution.     There  is  some  error  due  to  the  graphical 
techniques  required  in  the  solution,   however  the  agreement  is  close  enough 
for  engineering  work.     The  transient  response  may  now  be  calculated  and 
final  agreement  with  the  requirements  determined. 
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9 .  5    Summary 

For  both  the  examples  chosen,   the  lead-lead  feed  forward  difference 
compensator  proved  to  be  the  most  practical  for  solution.     However,-   the 
use  of  the  other  active  networks  follows  exactly  the  same  procedure  as 
demonstrated. 

It  was  found  that  there  are  some  regions  on  the  circular  loci  of  the 
active  networks  at  which  the  gain  ratio  A   /A     is  extremely  critical;  i.e.  , 
a  small  variation  in  the  gain  ratio  results  in  a  large  movement  of  the  com- 
plex zeros  of  the  active  network.     This  sensitivity  may  be  overcome  to  some 
extent  by  a  trial  and  error  method  of  varying  the  real  poles  and  zeros  of 
the  passive  networks  so  that  the  critical  area  of  the  circular  locus  is  moved 
away  from  the  desired  location  of  the  complex  zeros. 

The  active  networks  presented  afford  a  simple  method  of  generating 
complex  poles  and  zeros.     The  locations  of  these  complex  poles  and  zeros 
may  be  accurately  predicted  and  are  easily  varied  by  adjusting  an  ampli- 
fier gain  or  one  or  more  of  the  passive  elements  of  the  networks.     The 
varying  of  an  amplifier  gain  is  the  more  attractive  of  these  two  methods  as 
it  results  in  the  complex  poles  or  zeros  moving  along  a  circle  which  is  pre- 
cisely defined  by  the  passive  circuit  elements. 

Since  the  location  of  the  complex  poles  and  zeros  is  a  function  of 
amplifier  gain     a  requirement  for  a  very  stable  amplifier  is  generated. 
This  requirement  may  be  relaxed  somewhat  if  the  compensated  system  is 
not  too  sensitive  to  the  compensator  pole  and  zero  locations. 
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The  lag-lead  feed  forward  summing  compensator  has  characteristics 

which  merit  special  attention.     Since  the  K     of  this  compensator  is  A     +  A 

c  \  c. 

the  possibility  presents  itself  of  using  this  compensator  as  the  main  trans- 
mission path  amplifier.     The  complex  zeros  resulting  from  the  network 
could  be  adjusted  to  best  suit  the  needs  of  a  specific  system  while  the  over- 
all gain  remained  independent. 
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10.     APPLICATION  TO  A  JET  AIRPLANE 

In  any  discussion,  the  argument  is  strengthened  by  an  analysis 
of  a  practical  problem,   feasibility  permitting .     The  agreement 

of  calculated  values  and  test  responses  of  the  actual  vehicle 
is  proof  for  the  skeptic.     This  section  will  use  the  arrays  presented 
in  previous  sections  and  substitute  the  aerodynamic  stability  deri- 
vatives of  the  North  American  Aviation    In      Jet  Trainer  (Navy 
designation  T2J-1).     This  particular  vehicle  was  selected  as  a 
model  because  of  the  wealth  of  flight  test  and  aerodynamic     data 

immediately  available  to  the  authors. 

The  discussion  will  be  limited  to  the  longitudinal  system 

where  the  aircraf  .  is  originally  in  equilibrium  at  a  level  flight 

attitude.     An  additional  comparison  will  show  the  differences 

between  a  two  degree  and  three  degree  of  freedom  calculated 

estimate. 

Considering  first,   the  following  condition  of  flight: 

h  =  30,000  feet,    J0  -  8892.7  x  10-7  -^^-     ,      M=.6,  Vs     9?ft/sec 

r  ft3 

2 
W  =  9660  pounds,   m  =   300  slugs,     B  =   13157.7  slug-ft 

The  aerodynamic  stability  effects  were  calculated  and  the  following 

values  found! 

Xu  =     -   .01264  Xw  =   .0384 

zu  =  -    .0924  Mw  =  -    .0288 
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Mu  =   .0004  Mw  =  -   .00134 

Mq=   -    1.603 
zw  =  -   1.092  Mf^  =  -  24.771 

zq  =-  4.00  10)^  =  -   1.400 

z^  =  -   142.0  10Y1=  -  172.00 

Xw  =  Xq=zw=M  Y[  =10u  =  10w=  10w=  10q=0 

Placing  these  numerical  values  in  a  determinant,   the  following 
form  results  as  shown  in  Table  10-1. 

Evaluating  the   Zl  of  this  determinant,   the  result  is  c\  =  (s   +172 s+  1.4) 
(s4+3, 5056s  +18.99122s2+,2535s+.0998)  (10-1) 

This  results  in  the  following  set  of  roots, 

s  =  -    .01       ~\ 

/elevator  motion 
s=  -   171.99y 

s  =  -   .006205+  j.  07232  (  phugoid  motion) 

s  =  -   1  .7466  +  j  3.9864  (short  period) 

If  only  a  short  period  estimate  is  required,   then  the  determinant 
of  Table   10- 2  is  used.     This  assumed  no  speed  change  (u)  during 
the  short  period  and  the  X  force  equation  is  eliminated. 

Evaluating  the  Z3of  the  determinant,   the  results  are: 

s  =  -   .01 

(elevator  response) 
s  =  -   171.99 

s  =  -  1.745  +  j   3.985  (short  period) 
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The  short  period  response  using  the  results  of    Zjt  in  Equation  10-1 


is 


2     IT 

^  =  q«A a ~   1.575  seconds  per  cycle 


and  the  damping  to  1  /  2  amplitude  is  . 


.69 

=   . 395  seconds 


1/2     '   1.7466 

From  the  North  American  Aviation  flight  report  NA-594-282,   the 
flight  test  results  for  the  T2J  at  conditions  similar  to  the  above 
calculated  flight  conditions  were  T=  1 .  6  seconds  cycle. 

t    .       =  .4725  seconds 

There  is  very  good  agreement  between  calculated  and  flight  test 
results. 

There  was  insufficient,  information  available  in  the  North  American 
report  with  respect  to  the  long  period  longitudinal  oscillation. 
Therefore,   a  second  set  of  level  flight  condition  derivativeswas 
omputed  using  the  same  aircraft  model. 

The  conditions  specified  are: 
h=  20,000  feet,  /?=  12664  x  10-7  slugs/ft3       M=.463       U   =480ft/sec 

2 

W  =9660  pounds  m=300  slugs  B=  13157  .  7  slug-ft 
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For  these  conditions  the  following  stability  derivatives  were 
found: 

xu  =  -   .01382  Zq  =  -  4.60 

zw  =  -    1 .425  Mq  =  -   1  .  87 

-4 
Mw  =  8.59  x  10  z  n=  -   130.0 

Xw  =   ,0222  MY\  =  -24.45 

zw  =  -   1 .265  10h=  -  2.02 

Mw  =  -   .03365  10)1=  -   179.0 

zw=*».0  10u=10w=10wl0q=Xw=Xcpb 

Mw  =  -    .00193  M^=zM=X»i  =  Xv^=0 

Thus  a  determinant  can  be  formed  as  shown  in  Table  10  =  3. 

Evaluating  the  resultant    A  ,   it  is  found  that  the  following  roots 

exist 

-  _      ni?R        ~A 

elevator  motion 

s=  -   .006698  +j.  10129  (phugoid  motion) 
s  =  -  2.0257  +  j   3.7756  ( short  period  motion) 
Using  the  results  from  Table  10-3,   the  period  of  phugoid 

=62.0  seconds  per    cycle. 


motions  is     T  =     2 


.10129 

.69 
The  time  to  dampen  to  1/2  amplitude  is  t  =  ,_       =   10  3.0 

seconds.     If  a  two  degree  of  freedom  estimate  is  desired;    it  can  be 

determined  from  the  above  table.     For  the  phugoid  motion,   the  esti' 

mate  can  be  made  by  considering  angle  of  attack  constant 
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and  the  moment  equation  equal  to  zero  in  all  nodes.     Thus 

Table  10-4  can  be  made  from  Table  10-3>     The  roots  from 

this  determinant  are: 

(10-2) 

^\=(s+178.9575)  (s  +   .0125)  (s  +.0069  +  j*0979)  (s+ . 0069-j . 0979) 
The  flight  test  information  for  the  phugoid  was  determined 
by  the  authors.     On  5  January  1961  s   a  T2J-1  (Navy  148204), 
piloted  by  McCamey  and  Borthwick  yielded  the  following  results 
for  h=20,000  feet,  M-.463: 

T  ~  64  seconds  per    cycle 
t    .      =   103  seconds 

These  results  show  excellent  agreement  with  the  calculated 
estimates  and  prove  the  system  to  be  very  valid  one  for  both  long 
and  short  period.     It  reaffirms  the  previous  knowledge  that. 

1,  Given  reasonably  accurate  aerodynamic  stability  derivative, 
an  accurate  estimate  of  dynamic  response  can  be  made. 

2.  For  conventional  subsonic  aircraft,   a  two  degree  of  freedom 
approximation  is  quite  close  to  the  three  degrees  of  freedom 
result  and  should  be  sufficiently  accurate  for  preliminary 
design  estimates.     Furthermore,   factoring  a  second  order 
polynomial  is  far  simpler    than  a  fourth  order  polynomial,, 

As  mentioned  in  previous  sections  the  performance  function  is 

stated  thus 

PF  =  °j  K  Aij  (10-3) 

A 
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It  is  noted  that  the  /\  values  have  already  been  found  in  solving 

the  natural  periods  of  the  airframe.     A  simple  additional  step 

yields  the  performance  function  desired. 

For  example,   if  the  output    tt= —    1S  desired  the    A.,  determinant 

ZiFe  ij 

is  formed  by  cancelling  the  first  row  and  5th  column  of  Table  6-1 
the  results  form  in  Table   10-5  is  shown. 

The  result  of  this  determinant  in  algebraic  form  where  the  zero 
value  terms  are  omitted  is: 


15     S 


M  Y]  (Zw  -  1)  -Zy\     Mw 


1 


*l 


+  s 


Xw 


Zw  +Z         (Mw  Xu-Mw) 


M  Y)  (Zw-Xu 


T 


Zw-  1 


(10-4) 


My\     (XwZu-XuZw)  +   ZV)     (MuXw-MwXu) 


A 


X 


Substituting  in  the  appropriate  values  for  the  T2J-1  at  h=20,000 
feet,   the  following  quadratic  equation  results 


■^=  s2    (24.20)  +s    (26.900)  +  (  .535038) 
1  5 


(10-5) 


Recalling  that  the  direct  path  function  for  all  cases  in  this  system 

is  one  ( 1)  ,   then  G   ,   =  1 . 
56 


Therefore 


Q 
AFe 


G56^5 


or 
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(10-6) 

©_     _    24.20  (s+   .0203)  (s+  1.0907) 

AFe      (s+.0125)(s+178.9875)(s+.006698+j  0    10129M s+2 . 0257+J3  .  7756) 

Thus  the  performance  function  for  the  pitch  output  with  relation 
to  elevator  force  is  known.     The  block  diagram  of  the  internal 
sensing  and  activating  equipment  can  be  added  as  additional  parts 
to  control  this  multiloop  system.     Such  an  example  is  shown  in 
Fig.    10-1. 

The  functions  Gl  constitutes  an  error  sensing  and  signal 
amplification  device  such  as  an  autopilot.  G2  is  the  internal  control 
power  such  as  an  electric  or  hydraulic  servo  mechansim  G3  is  the 
performance  function  0/  /AF& 

Assume  for  the  moment  that  two  input  signals  are  desired 

in  the  system.     For  example,   the  aircraft  may  also  encounter 

turbulence  in  the  form  of  vertical  gusts  (w)  as  well  as  changes  in 

elevator  force  (  /\Fe) .     The  equation  for  consideration  of  double 

inputs  in  general  form  is 

qout=/lX      -i—       G.  k      (A  q  +A       q  J,  (10-7) 

H  (k)       /\  J  J         1  j         ^lm         2j       ^2m  J 

This  last  equation  representing  multiple  inputs  is  shown  only  to 

point  out  the  ability  to  use  this  system  for  such  purposes.     It  will 

not  be  used  with  numerical  coefficients  at  this  time,   however 
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FIGURE     10-1     SIMPLE  BLOCK  DIAGRAM  OF  PITCH 


REFERENCE  IN  AN  AIRPLANE 
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This  section  has  demonstrated  the  practical  use  of  the  dis- 
cussion of  previous  sections.     It  has  shown  the  agreement  of 
calculated  values  and  actual  test. 

This  particular  airframe  is  stable,   but  more  often  supersonic 
vehicles  are  not.     While  these  calculations  were  limited  to  the 
longitudinal  system,   the  use  for  the  lateral  system  involving 
the  yaw  and  roll  channel  is  equally  valid  for  the  appropriate 
equations  of  motion.     The  system  as  developed  in  prior  sections 
shows  an  orderly   process  whereby  electrical  and  for  mechanical 
processes  can  be  combined  into  one  multiloop  block  diagram. 
From  this  system,   with  the  various  criteria  set  forth,   further 
research  in  multiloop  compensation  theory  can  be  attempted. 
The  airframe  problem  shows  one  practical  use  for  this.     There 
are  many  others. 
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APPENDIX  I 
BASIC  EQUATIONS  OF  MOTION 

The  basic  concepts  of  mass  in  motion  and  the  equations  of 
motion  can  be  derived  in  vector  notation.     This  is  done  using  mass 
particles  in  motion  and  determining  their  effects  in  two  separate 
axis  systems.     The  basic  reference  system  shown  in  Fig.   AI-I  is 
the  starred  axis.     This  is  the  non-  rotating  reference  system  and 
is  also  known  as  the  Galilean  Axis.      The  axis  system  is  fixed  in  the 
rigid  mass  but  may  rotate  about  the  fixed  space  axis  system.     It  is 
the  relative  system.     A  special  set  of  notation  is  used  here  in  this 
appendix  which  differs  from  the  body  of  this  work. 

The  definition  of  the  special  terms  in  this  appendix  is  as  follows: 
i      =         the  vertorial  unit  along  the  X  axis 

->    A 

j      =        the  vertorial  unit  along  the  Y  axis 


t   * 


the  vectorial  unit  along  the  Z  axis 


■**    A  * 

i       =  the  vectorial  unit  along  the  X     axis 

-**  A  ,  * 

j       =  the  vectorial  unit  along  the  Y     axis 

->*  £  # 

k     =  the  vectorial  unit  along  the  Z     axis 

0     =  the  origin  of  the  relative  axis  system 

0     =  the  origin  of  the  space  axis  system 

-*    A 

R     =  the  vectorial  distance  from  the  space  axis  origin  to  the 


relative  axis  origin 

the  vectorial  distan< 
mass  particle  in  space 


->    A 

r     =        the  vectorial  distance  from  the  space  axis  origin  to  a 
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7M 

fJ  =      the  vectorial  distance  from  the  relative  axis  origin  to 
the  mass  particle 

A 

dm  =     the  mass  particle 


P     =     the  vectorial  distance  from  the  mass  center  to  the  origin 
of  the  relative  axis  system 

n  A 

J  t-  —     rotational  veloctiy  vector  of  the  relative  axis  system  about 
the  space  axis;  an  absolute  velocity  vector 

L&   =     absolute  rotational  velocity  vector  of  a  rotor  whose  axis  is 
attached  to  XYZ  axis  system 

-*    A 

L     =    momentum  vector  for  the  entire  mass  body 

~>    A 

H    =    moment  of  momentum  vector  with  respect  to  the  space  axis 

->  A 
G    =    moment  vector  with  respect  to  the  relative  axis 

^   A 

F     =     the  summation  of  the  external  forces  on  a  body 

c  A 

j    =     the  summation  of  the  internal  forces  on  a  body 

This  appendix  will  derive  the  basic  equations  of  motion  for  any 
mass  moving  in  space.     After  the  basic  equations  have  been 
developed,   the  effects  of  control  actions  will  be  added       The  final 
equations  will  then  be  converted  to  standard  aeronautical  terminology 
which  is  used  in  the  body  of  the  work  for  aircraft  stability  considera- 
tions 

These  equations  have  no  doubt  been  derived  several  times  in  the 
past  one  hundred  years       However,    the  author  was  unable  to  find  any 
clear  derivations  in  any  of  the  reference  material  covered       There- 
fore,  the  derivations  in  this  appendix  are  submitted  for  the  reader's 
interest  and  for  background  information 
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Fig.   AI-I    MASS  BODY  IN  MOTION 


otor 


mass  cenjei, 
-Rigid  body- 


Consider 


P-*       *       =* 


ix  +  jy  +  kz 


UO  =  11  =  iw    +  j  QJL    +  k   t0_ 

12  3 


Momentum  vector  for  the  whole  body     L  =     /    dm  dr 

J  "dT 

~tr         ~^  -*-     "^         /^ 

—  =  (Vo  +  (  60x/0)  +  /^of  rotor) 

L=    /    dmdr    =jdm  (  Vo  +(  U>x  P)  +  /^  rotor) 

'  dT    / 

^> 

-7—   =   /dm  d    r      = -f  +  F  but  summation  of  internal  forces  equals  zero 

dt 


Hence 


dL 
dt 


dm    d    r  =  F 


dt 
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For  an  airframe  with  a  large  rotor  (i.e.   helicopter,   STOL  air- 
craft or  GEM)  consider  addition  of  the  rotor  to  the  airframe 
using  the  mass  center  of  the  aircraft  as  the  origin  of  the  axis  system 

dt         (dt 


Tp> 


transport  terms 


relative, 
acceleration 
term 


coriolis 
effect 


Vo  =  dR 
dt 

Consider 

->  >        ~~ ~>    ~>     -5>         -^>      — >       "^      ~~^>     -2> 

d   r    =  dVo      +_fLx  A-0-x(_a.x/?)  +   f>+  ZlLxf 


dt 


dt 


then 
2-> 


r.     =  i  "Vx  +  j*Vy  +  k*Vz  +  i*  (iLVz  -  ^Vy)  +  j*  (AVx-  -O-V 


dt 


+  k*  (J^-  Vy  -  _fl -Vx)  +  0+0+0    +    0  where  /^=  0 


m 


Hence  F  =  ma^r 

dt2 

i*  Fx  =  i*  (Vx  +  O-  Vz  ~0-3Vy) 

j*  Fy  =  j*  (Vy  +  -ft-    Vx  -  -A-  V  )  m       >    neglecting  rotor  effect 

k*Fz  =  k*  (Vz  +A    Vy  -A;  Vx)  m 


Therefore  in  scalar  form 

Fx  =  (Vx  +jf)-2  Vz  -H3  Vy)m 

Fy  =  (Vy  +.il3Vx  -  -0^  V   )  m 

Fz  =  (Vz  +-0:   Vy^-Qr  Vx)  m 
l  £ 


Concurs  with  Etkin 

4.4.3 

Also  Duncan  Page  64 
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//I* 


rotor 


Figure  AI-2  VECTOR  DIAGRAM  OF 
ROTOR  MOTION 


Y* 


h*      is  fixed  in  the  rigid  body.     The  rotor  spins  at  an  absolute 

angular  velocity  of  but  the  mass  center  of  the  rotor  system  is  at 

point  P.     Therefore,    M    and   P   =  0.     It  is  also  assumed  that 

c  c 

_ i)-=  constant.      Therefore,   the  additional  term-O-x  (  _Ti.x   0  ) 

'  c 

enters  the  problem. 


(Pi        6^        c^3 

#  fc  /*c 

x  l      y 
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A  x  (B  x  C)  =  (A.C)  B  -  (A  .  B)  C  is  a  basic  identity,   therefore 


jix  ^£X   p    =  (  J?^)_0_  (^O-.n.)  '/$*. 


Note:     coordinates  i,  j,  k  are  those  in  the  body  and  not  the  inertial 

or  earth  system.     The  values  of  Jlare  primed  for  the  same 

reason.     It  is  quite  possible  that  fs        will  be  0  considering 

•       cy 

alignment  of  the  rotor  on  the  xz  plane. 

F      =  F.  +  F_  =  m     a  +  m_a  Fn   =  airframe      F„  =  rotor 
T  1  2  1  2         1  2 

In  resolving  this,   it  must  be  noted  that  the  mass  centers  differ 
from  those  of  the  combined  masses. 

It  should  be  noted  that  the  equations  for  the  acceleration  of  a 
particle  also  contain  rotational  and  relative  velocity  terms,   however, 

I    and    /  =  0  as  these  particles  are  in  a  rigid  body.     Further,   the 
summation  of  particles  about  the  mass  center  makes 

dm         (  -0-x    P  )    +  .A  x      (  _H-  x    P)  equal  to  zero. 

Moment  terms  may  be  considered  using  Moment  of  Momentum.   (H*) 

From  this  the  differential  may  be  taken  and  the  moment  is  then 
derived.     Consider  H       to  be  the  moment  of  momentum  of  the  system 
with  respect  to  the  inertial  (non- rotating  axis).      The  derivation  will 
include  a  rotor  spinning  on  the  rigid  body. 
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Figure  AI- 3    MASS  DIAGRAM 
FOR  MOMENT  ANALYSIS 


M  =  mass  of  the  rigid 
body  less  rotor 

m  =  mass  of  the  rotor 

The  absolute  angular  velocity  of  the  axis  OXYZ  is  A=i/L+  j-^*-+ki\. 

The  absolute  angular  veloctiy  of  the  body  is  l/)=i  iP   +j  [/)  +  k  10 

The  relative  angular  velocity  of  the  body  is    [p       -  -Clf=  n( scalar) 

The  angular  motion  of  a  rotor  in  motion,   but  attached  to  the  body  will  be 

IaJ-  7L=  i   (u)j  -jcij)  +7  (u>2  -A,)  +  k  tu>3  -  .0^)  . 

Henceforth  this  relative  term  will  be  noted  as  K  (KAPPA)  . 

# 

The  equation  for  H     is: 

H     =  J    r  x  dr      dm  where  r  =  R  + /O 

y     at  /    -^ 

Recalling  that  the  momentum  equation  was  L  =     /  dm  dr 

y  "dt 
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then 

H*::RxL+|^xdr        dm 

dt 


The  term  under  the  integral  will  be  known  as  H 

— >    ~^>  s>     -£>   —r> 

Recalling  thatjir  =  Vo  +_Q.  x  P  +  f    (rotor)  then 

dt 


H  = J  I  f>     x  (Vo  +  (I±x/°)  +    {>  (rotor) 


dm 


This  can  further  be  resolved  into  the  moment  of  momentum  about  its 
own  Oxyz  axis;and  the  velocity  moment  with  respect  to  the  inertial 


axis  . 


^ 


— 5»       -—>>  .  ~~£>    -^>    ~-^>  x 

H  =  (M  +  m)    f  c  x  Vo  +  C~\  P  -X.H-K  P)     dm    +j. P    x  (Kx   f    )  dm 


Ho     4  (/&    x  (JLx    P)  dm 


Ho 


i  (A 


rel     /  ' 


x  (  K  x    /°    )  dm 


To  derive  moments  consider  recalling  H     = 


dH*  = 
dt       J 


dr  dr     )    dm  -1 

dT  dt     / 


(r  x  dr )  dm 

"dt 

(R  +  P)  x  d    r 

dt2- 


dm 


dr 

dt 


dr      =    -0  and  recalling  that    F  =  dL,  = 


dt 


dt 


2~> 

d    r        dm 

1? 


Then  substituting  these  in, 

-^>    ^    _>  ^>  ~^>  zf 

dH*  =RxdL+  G  where  G  = 


dt 


dt 


h 


d    r 


I-        dt      J 


dm 


However,   we  also  noted  that 

'   px  d7 
dt 
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Therefore- 

dH*  =  dR 

dt  dt 


->      ->      ~^>  ^> 

x  L     +  R  x    dL       +  dH 


dt 


dt 


Equating  the  two  results  of    dH*  it  is  seen  that 

dt 


\ 


/ 


-\      -^>        -^>    ~>         ^=>\  ^$>         ^> 

R  >vaL     +  G  =  dR      x  L  +  R\x  xdL     +  dH 

dt~~ 


s  dt 

Recall  that 

dt 

dR  =  Vo 
dt 

/\dt 

and  that 

-^  ^>  ~^  *>   "^ 

L  =  (M  +  m)  (Vo  +ilx  P     +   P     rotor) 


Now_dH  =  d_ 
dt        dt 

dH  =  (m  +  M) 
dt 


combining 


->    ^  ->   _>. 

(M  +  m)     U     x  Vo  +  Ho  +  Ho      , 
/     c  rel 


0  ■*     "^       -^ 

d/^  x  Vo  +  /?    x  dVo 

'    C  /    C        -r— 

"dt"  dt 


+  _d_ 
dt 


r 


Ho      +  Ho 


rel 


-^> 


m 


G  =  Vo  x  ( Vo  +/lx  /^  )  M  +  Vo  x  (Vo  +  ILx   Z9     +    /*    rotor) 

-^,  c  _ f      c        '    c 

+    (   ^     x  Vo  +  (  SLx  f>    )    x  Vo)    m  +  (  f>      xVo+  (Stxf  )  x  Vo)  M 


+  (M  +  m)  /*    x  dVo       +  d 


dt 


dt 


— > 

i>       " 

Ho 

+ 

Ho 
rel 

A  further  assumption  is  stated  at  this  point.     This  is  that  the  point 
of  rotation  of  the  rotor  is  fixed  on  the  body  and  remains  so.     The  shaft 
does  not  have  nutation  about  some  axis  but  remains  fixed  on  a  mount 
of  the  rigid  body.     Therefore,   the     r^    of  the  rotating  body  is  zero. 
Also,   logically,   in  a  rigid  body  the  mass  does  not  shift  hence      /? 
is  also  zero. 
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In  crossing  a  vector  with  itself  Vo  x  Vo  the  result  is  zero. 


Vo  x  Vo  = 


i  J         k 

V        V         V 
1         2  3 


V        V         V 
1         2  3 


=  Mv2v3-V2V3)+T(vlv3-vlv3) 

+  k(Vj  V2  -  V1V3   )  =  0 


Vo  x  Vo  =  0 

Furthermore,    crossing  two  dissimilar  vectors  in  opposite  order  yields 
opposite  signs  and  equal  magnitude 

J  k 


VoxT  = 
•  c 


v,  V„       V„ 

1  2  3 


P     Pi    A 


-=>> 


ft  x  Vo  =  f*  (  ^2  v3  -  ^3v2  )  if(   f  v,  -  ^  v3)  t  M  <«v2-/»  V  ) 


Thus  Vox     f>  +    /?   x  Vo  =  0 

Incorporating  these  two  facts,   the  equation  for  6  can  be  simplified  to 
the  following: 

G  =JL  Ho  +  Ho  +  (M    +  m)  (  0    x  dVo  ) 


dt        L 


r^> 


dHo      =  Ho      +  Si-  x    Ho  in  a  general  equation 
dt~~ 
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~> 


Thus  it  is  necessary  to  solve  dHo    to  find  G. 

dt~ 


Considering  Ho  = 


-5> 


-S>    ~^>. 


/O  — •>>    ^^ 


which  accounts  for  the 


absolute  velocity. 

Recalling  an  identity  of  vector  algebra: 

Ax  (B  x  C)  =  (A*C)  B  -  (A.  B)  C  then 


Ho  = 


2         2         2     ->         -^»  -> 

(x    +  y    +  z  )  (i  (X^  +  j  \j)2  +  k  t^3)  -  (x uJj •+  y  (^  +  z  V)  j 


-5*    -~>    -> 

(ix  +  jy  +  kz) 


dm 


-5> 
i 


(y     +  z  )  (>i)  -  xy  {jJl  -  xz  (,0, 


+  j  |(x     +  z   )   u)2  -  yzU^  -  xy  LOj 

f-       2        2  "I 

+  k        (x     +  y   )  (^>3  -  xz   u^  yz  U^ 


dm 


The  above  terms  can  be  recognized  in  part  as  moment  and  product 

of  inertia  terms. 

Ix    =      |  (y2  +  z2)  dm 

Iy    S      [  (x     +  z   )  dm 


Iz 


A     f ,    2  2, 

=      I    (x     +    y   )  dm 


Ixy  =     J      (xy)  dm 
Iyz  =     I       (yz)  dm 


Ixz  = 


(xz)  dm 


21<" 


Therefore: 

Ho  =  i  (  Ix  OOj   -  Ixy  CO    -  Ixz  uJ)  +  j  (-IxyU^  +  Iy  Ul  -  Iyz  u») 

+  k  (Iz  U53  -  Ixz  UD    -  Iyz    uJ  ) 


Considering  the  effects  of  the  aircraft  as  a  rigid  body  first 


dHo  =  Ho    +jTLx  Ho  and  for  a  rigid  body     I/O  =  -fl- 
dt 

dHo  =  i  (Ix  JX    -  Ixy-^-     -  Ixz  SL  )  +  j  (-  Ixy -A-    +  ly  O-  -Iyz  A) 
dt  12  3  12 


-£> 


+  k  (-IxyiX    -  Iyz  P-   +  Iz  -Hr)  + 

1  £  J 


.0-1 


J 

&2 


k 
-0-3 


(Ix-Ixy-Ixz)     (Iy-Ixy-Iyz)    ( -Ixz-Iyz+Iz) 


For  any  ordinary  aircraft  or  missile,   the  produce  of  inertia  terms 
Ixy  and  Iyz  can  be  considered  zero  because  of  the  mirror  symmetry 
aspect.    Ixz  may  be  zero  if  the  principal  axis  is  selected.     Thus  with 
no  loss  in  generality,    consider  Ixy  =  Ixz  =  Iyz  =  0 

Mb  ="?  [ixi\  +  iz  -0-2  SL3  -  iyA,  n.3  1 

+  j    [iyJL,  +  IxJlj   _0-3   -  Izf^    n^    i 

+  k  [izAj  +  iy^  -/L  -  ix/ij  a2  i 


dt 


21$ 


However,   including  the  effects  of  Ixz,   the  product  of  inertia  term, 
which  is  often  no'i    zero: 
dHo    =  i 

ly  A  ^  +  ixA,    A     -  Ixz  A, 2      +  Ixz /I,  2  -  Izfl    I*- 


dt 


ixA  r  ixz  A  3-ixz/li  A2  +  izA2/l3-iyA2  /T 


+  J 


1 


1 


1 


+t 


-   Ixz  fi    +  Iyrt    A     -   ly  fl       fL        +  Ixz  A      -A 

J  X  Xw  1  £*  t*  ~) 


/I 


Using  standard  aerodynamic  terminology,   the  following  items  are 

defined: 

A^Ix 

B#Iy 

C  =  Iz 

E  =  Ixz 

Thus 

->     — >  /? 

G  =     dHo       where  system  is  chosen  such  that    f        =  0 

dt  C 


G  =  i 


+  J 


+-t 


AP  -  ER  -  EPQ  +  QR  (C-B) 
("BQ  +  PR  (A-C)  +  E  (P2  -  R2) 
("CR  -  EP  +  PQ  (B-C)  +  EQR 


These  equations  concur  with  Duncan,    Page  66,   and  Etkin,    Chapter  4 
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To  consider  the  rotor  effects  it  should  be  remembered  that 

W  =  -**-+  K.      where|s   is  the  relative  motion  of  a  rotor  on  the  body. 

Thus  for  a  rotor  there  will  be  a  combination  of  terms  considering 

the  absolute  angular  velocity  of  the  OXYZ  system  plus    the  angular 

velocity  of  the  rotor  with  respect  to  the  body. 

Thus 

dHo  =  Ho  +  tO  x  Ho    whereU— -fr-+/C 
dt 

It  should  be  noted,   that  for  a  symmetrical  rotor  or  disk  the 
product  of  inertia  terms  are  zero  where  the  disk  is  at  least  a 
three  bladed  propeller.     The  two  bladed  propeller  has  a  product 
of  inertia  dependent  on  its  angular  position.     This  case  will  not  be 
considered,   but  the  three  or  more  bladed  propeller  situation  will 
be  used  in  all  considerations.     Hence  with  a  dynamically  and 
statically  balanced  rotor    r   =  0  and  the  product  of  inertia  is  zero. 

The  moment  of  inertia  of  the  disk  must  be  considered  from 

the  axis  of  the  rigid  body.        Therefore,   the  moment  of  inertia 

2 

must  include  the  md     consideration. 
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I„_  =  IA  A   +  md 
BB        AA 


The  moment  of  inertia  terms  for  the    rotors  will  be  thus  defined 
in  the  most  general  sense 

a  =  Ix  rotor  =  Ix  +  md       and  aligned  with  the  reference  axis 

A  A  2 

b  =  Iy  rotor  =  Iy  +  md        and  aligned  with  the  reference  axis 

c  =  Iz  rotor  =  Iz  +  md        and  aligned  with  the  reference  axis 


G  =  dH 
dt 


dHo    +  dH 
dt 


rel 


dt 


Recalling  the  derivation  of  moments  for  the  rigid  body  without  the 
rotor,   the  same  method  applies. 
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5> 


dHo       =  i  a  A-  +  jb  A,  +  kc  XI    +  -fix  H 

1  2  3 

dt 

^>      •         __,     »        _$>  —5J,      -*>     _>>  for  a  rotor  only 

+  iaKi+jbl<2+kcK3+<xHo 


Combining  this  to  a  system  consisting  of  a  rigid  body  plus  a  rotor, 
the  following  results  occur 

G  =  dH    = 


dt 


+  J 


+  k 


^5> 


ix  i^    -  ixzi^  +  ixz  JX  Jl  +i     Jl2  iT  -  Iy  -O    A. 

Iyi\+  Ix  rijHj  -  Ixz/l2     +  IxzfL2  -Izf^    J\ 
Izjfl3  -  Ixzfl    4-  IyXl     ^-+  Ixil  -^-+I.xz  A2  -H- 


+  i  a^  +  jb  )<2+"kc  ^3  + 


-2>  -> 

1  J 


ill  _fL2  Jl3 


^1  b<2  4 


Combining  and  using  aerodynamic  notation, 


t=T 


+1 


+  k 


AP  -  ER  +  a  \L  -  EPQ  +  (C-B)  QR  +  ck  Q  -   b/^R 
BQ  +  b  K^  +  (A-C)    PR  +  E(P2  -  R2)  +  a  £  R  -  c  ^P 

CR  -  EP  +  c   K     +  (B-A)  PQ  +  EQR  +  b)C.  P  -  a£  Q 


This  result  concurs  with  Etkin,   page   116. 

Therefore,   using  a  standard  aerodynamic  notation  the  moment: 

M  +  kN 
and  the  scalar  quantities  are  as  follows 


G =  i  L+  j 


2 


L  =  AP  +  a  K     -     E  (PQ  +R)  +  QR  (C-B)  +  cK,Q  -  b  K_R 
1  3  2 


M 


+  BQ  +  b<2  +  PR  (  A-C)  +  E  (P2  -  R2)  +  aft    R  -  cK_ 


R 


N  =  CR  +  c  K3  +  E  (QR  -  P)  +  PQ  (B-A)  +   b^P-a^Q 

Thus,   the  force  and  moment  equations  derived  in  this  appendix 
are  applicable  to  a  rigid  body  of  constant  mass  where  the  body  is 
in  motion.     It  also  considers  the  effects  of  a  three  bladed  or  more 
rotor  system  which  is  attached  on  some  fixed  point  of  the  rigid  body, 
The  rotor  effects  may  be  included  or  discarded  depending  on  the 
magnitude  of  rotor  inertia  and  the   relative  velocity  of  the  rotor. 

It  should  be  noted  that  the  derivations  are  done  in  a  general 
form,   but  in  the  concluding  equations,   the  usual  aerodynamic  nota- 
tion is   substituted  to  bring  the  equations  to  this  specific  problem. 
Furthermore,   missile  or  aircraft  shape  and  density  has  been 
considered  symmetrical  about  the  xz  plane.      Thus,   the  produce  of 
inertia  terms  Ixy  and  Iyz  have  been  considered  zero.     No  generality 
is  lost  by  this  and  if,   for  some  reason,   the  body  is  sufficiently 
unsymmetrical ,   the  terms  can  be  quickly  added  in. 

The  following  is  taken  from  Chapter  4  of  Etkin  (4)  .     The  concept 
of  motion  on  the  control  surfaces  is  shown  in  the  following 

n         ,*" 


Contr 
Lever 


» 


Hydraulic 
Cylinder 


A  SIMPLE  ELEVATOR  CONTROL  SYSTEM 
2°,0 


The  LaGrange  equation  of  motion  in  a  moving  reference  frame  is 


dt      V 


q 


1 


^T    =.7 


<?qic 


T  =  kinetic  energy  of  the  system  relative  to  the  chosen  frame  of 
reference 


J^ 


^w     =  generalized  force 


W  =  work  done  on  the  system  by  the  external  forces  which  act  upon 
it 

q-|/=  generalized  coordinate 

H  =  aerodynamic  moments 


f- 


rotation  of  control  surface  in  radians 


"  2  *r 

Kinetic  energy  T  =  1/2    l£     and      ^T    =  l£ 


^=     ^W 


For  elevator  motion 

^W  =  He/^e)  +    FeoT(J"e)  +  ^Wi 
Fe  =  generalized  elevator  control  force 
A  Wi  =  work  done  by  inertia  forces 
Work  done  by  the  control  forces  is 

Fed  (cfe)  =  Pcf  sp  +  J(f  sj 

Fe  =  P  dsp        +  J  dsJ 


dc/  e  <Ve 

=  gear  rate  of  pilot's  stick  to  control  surface 


dsp   /i 

-JL    =  "gear  ratio"  or  advantage  of  hydraulic  cylinder  to  control  surface 

22  1 


Evaluating     £  Wi 


2 
consider  dFi  =  d    r    dm 

"dt~ 


recalling  from  previous  derivation  that 

~f>  ->      -=>      ->  ->     -^    ^>     ^P      ^>    ~P 

jr     =dVo    +  fix  /^+H-xilx/?  +  r   +  2ilx/<? 


dt  dt 

Then 
dFxi  =  - 


dVo        +  2  Qz  -  2  Ry  -  x  (Q     +  R2)  +  y  (PQ  -  R)+z(PR-Q) 


dt 

dFyi  and  dFz  are  similar  in  form. 

Looking  at  an  elevator  assembly,   it  can  be  seen  that  assuming  a 
lamina  in  the  xy-plane,   the  displacement  in  the  C      direction  is 
amount    <$      (  £  e)  . 
The  work  done  by  the  inertia  forces  is    <3  W     =    (^    A  (  o  e)  dFzi 

qT Wi         =  f     dVoz      A:  dm  -  (PR-Q)   j   x    /^  edm  -  RQ  +  P  \    y/om 

Jwr.)  /    dt  y 

The  coriolis  terms  2  Py  and  -2  Qx  are  zero  as  there  is  no  work  done 
in  those  directions.     Furthermore  z  =  0  because  of  the  assumed  lamina 

acz  =  dVoz    and  acz   /    r      e  dm  =  m    1    a 


e    e    c 


dt  '  z 


Q)  f   x   fe 


(PR  ■-  Q)     /       x   /*e  dm  =  P       (PR  -Q) 

x 

P         product  of  inertia  of  the  elevator 
e 
x 
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cT 


If 


By  symmetry      I    f       Y    dm  =  0 


e ' 


Wi         =  -  m    1    a        -    P        (PR-Q) 
e    e    c  ex 

z 


1    cf       =   &*   m    1    a         -    P        (PR-Q)    +H        +    F 
e        e  eec  e  e  e 

z  x 

7 

If/  +         m    1    a        +  P        (PR    -Q)    =H     +  F  (Etkin  4   8,    13) 

e^e  eec  ex  ee 

z 


Similarly  for  the  rudder 


*» 


I     0       -ami  -  (R  +  PQ)  P        -  (RQ  -  P)  P       =  H     +  F 

r  &    r  cv     r    r  rx  rz  r  a 


and  for  the  aileron 


I     _(  +  2  P        (RQ  +  P)  =  2  H     +F 
a  ^  ay  a  a 
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APPENDIX  II 
ROLL,   YAW  AND  PITCH 

The  derivation  of  the  yaw,    roll  and  pitch  rates  is  given  here  to 
show  the  actual  motion  taking  place       Frequently  the  motion  is 
considered  small  and  small  angle  approximations  are  used. 
These  will  be  noted  at  the  conclusion. 

First,    consider  the  Eulerian  Axis  used  in  the  aeronautical 
field  shown  below.     The  rotation  takes  place  in  the  order  of  yaw 
(  ty   ) ,  pitch  (0  )    ,   and  roll  {0  )  . 


x?.X 


£Lo,Z-{ 


Figure  All- 1     EULERIAN  AXIS  SYSTEM  SHOWING  ROTATION 
ABOUT  THE  Z,    Y,  AND  X  AXIS 

224 


The  initial  conditions  at  X    ,    Y    ,   and  Z     and  the  relations  of 

o         o  o 

Of  ,  ©,   and  0  can  be  related  to  P,   Q  and  R  of  the  body  axis  at 

x,   y,    z.     In  the  above  diagram  consider  the  angular  velocity 

P  which  is  about  the  OX  axis.     It  must  be  the  sum  of  the  velocities 

about  the  OZ    ,   OY,  ,   and  OX_,  axis  respectively, 
o  1  2 

rr 

OZ     makes  an  angle    — - —      +  -0  with  OX.     The  rate  thus 

^  cos  (  4|L     +  e)  =  -  ipsin  e.     OY     is  perpendicular  with  OX 
TT 

so  ©  cos      — - —     =  0    and  OX  is  parallel  with  OX     or  the  motion 

0  cos  ©  =  0  . 

Thus  the  summation  P  =  0  -     ty  sin  0. 

To  evaluate  Q  the  angle  YOZ     must  be  determined.     Specifically 

y\ 

the  value    W  cos  YOZ     is  desired. 
'  o 

/\  /\  /<>X 

cos  YOZ     =  cos  YOZ,  cos  Z   OZ 
o  2  2        o 

IT 

-    Cos    (  — y-      -    0)    COS    0 

=   sin  0  cos  © 


Thus  the  contribution  of  OZ     axis  motion  is    Y  sin  0  cos  ©. 

o 

The  angle  between  OY  and  OY     is  0  so  the  motion  here  results  as 

©  cos  0.     OY  is  perpendicular  to  OX     hence  the  angular  motion 

*  *7T 

0  cos  -4f—     -  0. 

Thus  Q  =  ©  cos  0  +      ^     sin   0  cos  © 

In  the  third  angular  velocity  motion  R  the  yaw  motion  appears  as 

*  »  "TT  • 

(f)  cos  0  cos  0  while  the  motion  of  pitch  is  0  cos  (-y-  +  O)  =  -0  sin 

0.     The  rolling  term  vanishes  0  cos  -y—    =  0.     Thus  R  =  -  ©  sin  0 

+   CPcos  0  cos  ©. 
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Hence  the  relationship  of  the  body  yaw,   pitch  and  roll  rates 
(P,   Q,  and  R)  to  the  designated  reference  system  angular  velocities 

QJ ,   0  and  0  is: 

*         • 
P  =  0  -    Cp  sin  0 

Q  =  0  cos  0  +    QJsin  0  cos  0 

R  =  tp  cos  0  cos  e  -  0  sin  0 

These  results  concur  with  Duncan  page  81  and  Etkin  page  116. 

It  should  be  noted  that  for  small  angles  the  approximation  below 

is  often  used: 

P^0 

A  second  consideration  which  must  be  made  regarding  movement 

about  the  reference  axis  is  the  velocity  or  displacement  with  respect 

to  time.     It  is  desired  to  obtain  the  velocity  with  respect  to  the 

fixed  frame  or  inertial  axis  I*  x,  j*y,   and  k*z. 

Consider  first  that  in  scalar  values 

dx*    =U 
dt 

dy*    =V 

dt 

dz*    =W 
dt 

It  should  be  noted,   however,   that  in  the  three  position  system  shown 

on  the  previous  page  that 
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u3  =  u 

V     =  V  cos  0  -   W  sin  0 

W     =  V  sin  0  +   W  cos  0 

and  W     =  -  U     sin  -0  +  W     cos  0 

v„  =  V 
2  3 

U     =  U     cos  ©  +  W     sin  0 
and  U   =  U     cos  [^  -  \     sin  (V 
Vj  =  U2  Sinf  +      V2  cos  ^ 

Wl  =  W2 
Combining  the   1  ,    2  and  3  position  terms 

dx*  =  U  cos  0  cos  >-p+  V  (sin  0  sin  0  cos  ty-   cos  0  sin  ty) 
dt 

+  W  (cos  0  sin  ©  cos  (fS+   sin  0  sin  (f> ) 

dy*  =  U  cos  0  sin  if +   V  (sin  0  sin  ©  sin^+  cos  0  cos  V^ 
dt 

+  W  (cos  0  sin  ©  sin  (f~   sin  0  cos  Lp  ) 

dz*    =   -U  sin  ©+  V  sin  0  cos  ©  +   W  cos  0  cos  © 
dt 

Thus,   the  velocity  or  rate  of  change  with  respect  to  the  inertial 
or  earth's  reference  axis  in  this  case  is  expressed  in  these  terms 
as  derived.     These  velocity  terms  are  derived  in  Etkin  Page  102 
and  conincide  with  the  derivation  shown. 

Therefore,   this  appendix  accounts  for  the  changes  in  angular  and 

directional  displacement  of  a  rigid  body  with  respect  to  a  fixed 

reference  axis  system. 
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APPENDIX  III 
NOTES  ON  DETERMINANT 
MANIPULATION 
The  body  of  this  project  contains  several  determinant  expan- 
sions where  the  order  is  doubled.     These  manipulations  do  not 
alter  the  values  of  the  determinant  and  the  resulting  values  of  the 
performance  functions  are  the  same. 
The  simple  determinant 

Al  Bl  CI 

A2  B2  C2 


A3 


B3 


C3 


A- 


evaluated  will  yield  the  following  result 
A1B2C3  +  B1C2A3+  C1A2B3 
-C1B2A3  -  A1C2B3  -  B1A2C3 


=  A 


A3-1 


This  expansion  by  the  equations  -1  +   1  =  0  will  result  in  the  following 
determinant. 


Al 

0 

Bl 

0 

CI 

0 

-1 

1 

0 

0 

0 

0 

0 

A2 

B2 

0 

C2 

0 

0 

0 

-1 

1 

0 

0 

0 

A3 

0 

B3 

C3 

0 

0 

0 

0 

0 

-1 

1 

The  resulting  evaluation  will  yield  the  same  value  as  result  (A3-1). 
Furthermore,   if  the  terms  below  the  main  diagonal  are  allowed  to 
remain  in  their  original  columns  in  lieu  of  shifting  one  place  to  the 
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right  as  is  done  in  the  main  body,   the  resulting  Za.is  still  the  same 
The  results  of  evaluating  cofactors  for  particular  performance 
functions  remain  unchanged.     For  example  it  is  desired  to  evaluate 
the  cofactor    ^—^-,7  of  the  original  determinant.     Thus: 

1  w 


-    -Al Bfi    —    -eh 


12 


A2 
A3 


B|2 
b'3 


C2 
C3 


=  (-1)      ^A2C3  -  A3C2/  A3-2 


In  the  expanded  determinant  the  equivalent  cofactor  i 


is  A 


13 


-  -f  A-i &    —      -B'l  — 


1 

A2 

0 

A3 

0 


b;2 


bi3 


-Cl- 
0 

C2 
0 

C3 
-1 


-  0 
0 
0 
0 
0 
1 


=  A 


13 


1 

A2 
0 
Ae 


0 

C2 
0 
C3 


A2        0 


A3        0 


C2 


C3 


A^  =  (-1)      ^2C3  -  AeC2j 


A3- 3 


2^9' 


The  result  of  Equation  A3- 3  equals  Equation  A3=2.     With  the 
values  below  the  main  diagonal  left  in  their  original  columns, 
the  following  occurs: 

Al  0  Bl  0  CI  0 


A 


13 


-1 

1 

0 

0 

0 

0 

A2 

0 

B2 

0 

C2 

0 

0 

0 

-1 

1 

0 

0 

A3 

0 

B3 

0 

C3 

0 

0 

0 

0 

0 

-1 

1 

By  inspection  of  columns  two  and  six,    two  reductions  are  made 
to  the  following: 

A2  0  C2 

0  1  0 

A3  0  C3 


13 


/\13  =  (-D     £A2C3  -A3C2^ 


A3-4 


The  resulting  equation  (A3-4)  is  the  same  as  equation  (A3-2). 

A  similar  examination  using  a  four  by  four  determinant 
system  was  made.     This  was  done  to  insure  that  the  simplified 
effects  of  a  three  by  three  or  lower  order  system  did  not  cover 
up  an  error  in  a  four  by  four  or  higher  order  determinant.     The 
same  manipulations  can  be  performed  and  the  resultant  value  of 
the  determinant  is  not  altered. 
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Thus,   the  individual  has  the  option  of  placing  his  values 
m  either  column.     The  location  of  the  values  is  dependent,  on 
desired  location  of  the  pick-off  point  on  the  block  diagram. 
The  final  result  is  not  altered 
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